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ON THE QUALITY OF COLOR-TELEVISION IMAGES 
AND THE PERCEPTION OF COLOR DETAIL*t 



By 

0. H. Schade, Sr. 

RCA Election Tube Division, 
Harrison. N. J. 

Summary— A theoretical, and experimental study of the NTSC color 
system, supported by color photographs, shows that contrast range and 
color saturation obtained with commercial tricolor kinescopes provide a 
larger color space than provided by color motion pictures. 

In fine detail more than 60 per cent of full color information is trans- 
mitted and reproduced by the NTSC system because the bandwidth restric- 
tions of the electrical color signals (I,Q) do not affect definition in the 
vertical dimension and have a smaller effect on the reproduction of horizontal 
color detail than indicated by earlier evaluations which disregarded the 
two-dimensional nature of the image. 

The. detail color reproduction appears adequate to the eye because the 
color errors remaining, although perceptible, are small. This fact is 
significant because the spatial sine-wave response functions of _ the color 
discriminators of the visual system are found to be substantially inde- 
pendent of the color of light and similar to the spatial sine-wave luminance 
response function of the eye. 

PART I— ELECTRICAL AND OPTICAL CHARACTERISTICS OF 
NTSC COLOR TELEVISION SYSTEM 



General Characteristics of Color Systems 

XT 



-SJ HE reproduction of color in a television or photographic system 
based upon the trichromatic theory. The analysis of color 
a television or photographic camera requires a discriminator 
mechanism having three different spectral sensitivities, resulting in 
three integrals or "primary signals" R, G, B"* (or any linear trans- 
formation thereof, such as the NTSC values Y, 1, Q) which are 
independent two-dimensional intensity functions of the object-point 
coordinates. In the synthesis of a colored image, luminance and color 
of all image points is restored by letting the primary signals control 



* Manuscript received November 10, 1958. 

t This paper was presented at the Convention of the Society of Motion 
Picture and Television Engineers of Detroit, Michigan, on October 22, 1958. 
It is also being- nublished in the December, 1958 issue of the Journal of 
that Society. 

** The spectral sensitivities or wavelength functions /(X) should be such 
that the primary signals are approximately proportional to the tristimulus 
values of the object, the tristimulus values being those of the receiver 
primaries. 
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three properly chosen "primary lights'": the reproducing primaries 
red, green, and blue. The reproduction of color by a linear or a non- 
linear system must be independent of illumination intensity to conform 
with the requirements of vision. A stable color balance necessitates 
a constant ratio of the three signal functions, which requires: 

(a) Matched transfer characteristics, 

(b) Matched spatial frequency spectra of the system, and 

(c) Matched noise-levels. 

These specifications express a registry requirement. The registry of 
three transfer characteristics is an old problem in photography and a 
new problem in television, demanding high precision in maintaining 
relative and absolute gain and black-level stability in linear or non- 
linear amplifiers as well as in camera tubes and kinescopes, where 
sensitivity or response within the image area must be highly uniform. 
Registry of the spatial frequency spectra of three color images re- 
quires accurate time-delay and phase equalization in electrical channels 
and a closely matched image geometry in three scanning rasters. 
Matched noise levels become important when the signal-to-noise ratios 
are low, because unmatched levels result in a mismatch of transfer 
characteristics. 

For the reason of "compatibility" in television, the color informa- 
tion, together with a luminance signal suitable for monochrome 
receivers, must be transmitted without increase of video bandwidth. 
The two additional color-transmission functions (/, Q) are, therefore, 
restricted to smaller passbands than the luminance function (Y). 

The present analysis does not concern itself with faults of a tem- 
porary nature in the operation of a color system, but with the range 
and purity of colors obtainable, in particular with the unequal pass- 
bands of NTSC television signals and with commercially available 
color kinescopes. The computed colorimetric performance is checked 
by measurements, supported by visual observations (color photo- 
graphs), and interpreted with reference to the performance of color 
photography and the visual system. 

Color Television System for Measurements and Visual Tests 

The color signals are generated by a light-spot slide scanner (see 
Figure 1) which is followed by three gamma-correction amplifiers and 
band-limiting filters feeding the cross mixing network or Matrix I 
(RCA "colorplexer") for translation of the input signals E R , E a , 
and E B into E Y , E h and E Q signals. The matrix coefficients for the 
system are given in Table I. After band limitation by / and Q filters 
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PHOSPHOR 

CORRECTION FILTERS ■ 

PICKUP UNIT GAMMA (4 Mc) 



LIGHT-SPOT 
SCANNER 



OBJECT 
TRANSPARENCY 




FREQUENCY TRANSFER 
SPECTRUM 



Fig. 1 — Block diagram of color system. (Modulator and demodulator system 
(not shewn) can be inserted ahead of Matrix II at point M.) 



in the colorplexer, the signals are translated back to E' R E' a E' B in 
the inverse Matrix II. The primes indicate that each signal is now 
a frequency-dependent mixture of components from three unequal 
passbands as indicated in Table I and discussed later. These signals 
are then modified by the transfer characteristics and frequency char- 
acteristics of the electron guns of the color kinescope and superimposed 
on the kinescope screen, where they are converted to light energy. 
Modulation and demodulation circuits (not shown) can be switched 
in at point M to observe effects introduced by transmission of the 
multiplexed NTSC signal over a single channel. 

For measurements of the dynamic transfer characteristic of the 
system, the light intensity in a fixed small area on the kinescope screen 
is observed by a multiplier phototube while a calibrated step tablet 
image is slowly drifting over this area in a vertical direction. The 
step signal is generated by a color transparency having the charac- 
teristics indicated in Figure 2. The slow vertical drift is generated 
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Fig. 2 — Test pattern for measuring transfer characteristics. 
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by synchronizing the vertical oscillator of the light-spot scanner with 
a stable oscillator differing very slightly in frequency from the field 
frequency. The screen of the kinescope is thus illuminated continu- 
ously in its entire area by the drifting colored test pattern picture, 
while the monochrome step intensities are traced on a recorder. The 
optical step signal is generated in two sections. The main step tablet 
(a) covers a 70-to-l range, and a 10-to-l neutral filter strip (b) extends 
this range to 700. Measurement of the electrical signal steps at inter- 
mediate points of the system furnishes the various transfer charac- 
teristics shown in Figure 3. 




LOG RELATIVE INPUT SIGNAL (B|/B| MAX..E|/E| MAX.. OR E 2 /e 2 max} 

Pig. 3 — -Transfer characteristics of color system. 

Picture-Tube Transfer Characteristics 

The dynamic transfer characteristic shown in Figure 3 was meas- 
ured on a 21AXP22 metal color kinescope. It follows a power law 
(close to a square law) with an additive constant B 0 , which expresses 
the black-level illumination or light-bias on the screen resulting in a 
"toe" in the characteristic. The light bias is caused by diffuse electron 
excitation (secondaries) , optical diffusion in the screen, and ambient 
light. It is variable, therefore, and has its lowest value for a dark 
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viewing room and for pictures having a high ratio of peak-to-average 
luminance. For normal picture material, the range of the transfer 
characteristic is similar to that measured with the colored test pattern 
(Figure 2), i.e., approximately 600 to 1 in a dark room." This range 

is much higher than obtained in motion picture theaters (60 to 100 
to 1) where the light bias, B 0 , is caused by projection lens flare and 
ambient light; in fact it is even higher than the contrast range of 
the color film positive itself* which is in the order of 470 to 1.** The 
large range of the color kinescope cannot be photographed in its 
entirety on any Kodachrome or Ektachrome film and is further reduced 
by the color printing process in the various illustrations. 



x 

Fig. 4 — CIE diagram showing' location of television receiver primaries 
(R,G,B), I and Q axes, and color planes. 



The degree of color saturation obtainable in a color image is deter- 
mined by the reproducing primaries and the contrast range of the 
reproducer. The red, green, and blue primaries of the color kinescope 
are located in the CIE diagram as shown in Figure 4. A straight line 
drawn through the white point (illuminant C) intersects the spectrum 
locus at two points. The excitation purity, S, of a color located on 
such a line is, by definition, 1 the ratio of the distance from the color 
point to the white point, to the distance from the spectrum locus to 
the white point. The excitation purity is thus specified by a mixture 

* The reader may compare these values with those given in Reference 
(2), showing- transfer characteristics of early (1954) color kinescopes 
covering- a range of hardly more than 20 to 1. 

"'"Today's color kinescopes (21CYP22) have an even higher contrast 
because of reduced electron diffusion by secondary emission. 

1 Dean 15. Judd, "Colorimetry," National Bureau of Standards Circular 
U78, March, 1950. 




Excitation Purity (Color Saturation) 
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of a spectral light with white light. The fixed light bias, B 0 , caused 
by ambient or scattered light will be assumed to have the neutral color 
of illuminant C. The addition of this white light moves all color points 
towards the white point and decreases the excitation purity of the 
color from the value S e , computed for £„ = 0 (see Table II), to a 
smaller value, 5. The purity reduction by a fixed white light energy 
is obviously dependent on the excitation level of the color and its 
relative stimulus energy, defined as follows : The maximum luminance 
Y max in a given color system is that of the white light (F max = Y w max ) 
obtained when the three primary lights (phosphors) are excited to 
selected maximum luminance values (F rmas , Y g max , Y b raax ) giving the 
specified white light (Illuminant C for color television). The maximum 

Table II — Tristimulus Energy Factors, w„ Luminance Factors, U, and 
Excitation Purity, St, NTSC Standards 

Dominant X Trichromatic Coefficients 



Color 


We 


I. 


(millimicrons) 


S, 


X 


V 


z 


W 


1.00 


1.0 




0 


0.31 


0.316 


0.374 


R 


0.286 


0.299 


eii 


1.0 


0.67 


0.33 


0.00 


G 


0.261 


0.587 


535 


0.85 


0.21 


0.71 


0.08 


B 


0.453 


0.114 


470 


0.88 


0.14 


0.08 


0.78 


Y 


0.547 


0.886 


573 


0.88 


0.45 


0.512 


0.038 


Mo 


0.739 


0.413 


535 compl. 


0.625 


0.345 


0.177 


0.48 


Cv 


0.714 


0.701 


490.5 


0.52 


0.165 


0.31 


0.525 




Munsell Lightness Value, V, 


and Luminance, 


Y 






V 


10 


9 8 


7 


6 


5 








102.56 


78.66 59.1 


43.06 


30.05 


19.77 






V 


4 


3 2 


1 


0.5 


0.2 






Y(%) 


12. 


6.56 3.126 


1.21 


0.581 


0.237 





luminance of a color, c, is therefore obtained when at least one of the 
primary lights reaches the excitation limit established for maximum 
white light excitation. 

The excitation level of a color is hence specified by the excitation 
factor 

fc c =F c /y omal gl. (1) 

The maximum luminance of a color relative to the maximum luminance 
is expressed by the luminance factor 



(2) 



and the relative tristimulus energy of a color is given by the product 
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k,w c = k c (X+Y + Z) c max / (X + Y + Z) w max . (3) 

The tristimulus energy factor, w r , can be expressed in terms of tri- 
chromatic coefficients, y, and the luminance factor, l r , by the simple 
relation 

wjl c = y w /y 0 . (4) 

The tristimulus energy factor, w„, of the light bias is specified by the 
maximum contrast ratio C, 

W Q = wJC = 1/C. (5) 

The mixture ratio k c w c /w 0 of the relative tristimulus energy of the 
color to the fixed relative tristimulus energy of the light bias is equal 
to the distance ratio S/(S C — S) , which gives the desired purity relation 

S/S c = k a wj (k c iv c + tv 0 ) , (6a) 

and, with Equation (5), 

S/S e = k c wj (k c w c + 1/C) . (6b) 

The total relative luminance of the color and light bias energy is the 
sum 

Y/Y max = k c l c + 1/C. (7) 

Inspection of Equation (6a) shows that only an ideal color reproducer 
having an absolute black level (tv 0 = 0) provides a constant excitation 
purity for all colors, independent of excitation, k c . A plot of excitation 
purity, S, as a function of relative luminance (F/F ma J for all possible 
colors furnishes a color space in which arctan y/x is the vectorial 
direction of the color from the white point (see Figure 4), excitation 
purity is the vector length, and relative luminance is the elevation. 2 
The axis of this space, erected over the white point, represents the 
gray scale, and vertical sections through this axis furnish color planes 
such as shown in Figures 5-7. The boundaries of these color planes 
include all possible colors and permit a comparison of different 
reproducers. 

Given the NTSC primaries and constants (Table II) and the con- 
trast range, C, of the picture tube, the lower boundary of television 

2 w - L - Bl ^ wer . J - H. Ladd and J. E. Pinney, "Brightness Modification 
1954° S Televl smg Color Film," Proc. I.R.E., Vol. 42, p. 174, January, 
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EXCITATION PURITY S 
1.0 0.5 0 0.5 1.0 




Pig. 5 — Red-cyan color plane of additive television process (shadow-mask 
kinescope C = 500), solid lines, and subtractive M.P. process {C = 71) 

broken lines. 2 

color planes is constructed by selecting a color on the color triangle 
(Figure 4) and computing the coordinates S, (Y/Y max ) of points on 
the boundary curve with Equations (6) and (7) by assigning values 
between 0 and 1 to the parameter k c . 

The computed boundary extends up to the relative luminance l c + 
1/C at full excitation (k c =l) of the primary color (or two-component 
mixture) . Higher luminance values in the particular color plane can 
therefore be obtained only by adding light of complementary color 
which combines with a certain fraction k' of the color energy to a 
partial white excitation, having the relative stimulus energy fc'ttv 

EXCITATION PURITY S 



1.0 0.5 0 0.5 1.0 




Pig. 6 — Green magenta color plane of additive television process (shadow- 
mask kinescope C — 500), and subtractive M.P. process (C = 71) broken 

lines. 2 
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The remaining relative stimulus energy of the color is (1 — k')w c . 
The excitation purity is hence determined by the mixture of the re- 
maining color amount (l — k')w c with the amounts of white k'w w + 
w 0 = k' + 1/C, which leads to the expression for the upper boundary 
of the color space; 

S/S c - (1 - k') wjl (1 - k') to, + k' + 1/C], (8) 

where k' Si 1 (white excitation factor) . 

The relative luminance is given by 

Y/Y max = (1 - fc') l c + k' + 1/C. (9) 

EXCITATION PURITY S 
1.0 0.5 0 0.5 1.0 




Fig. 7 — Yellow-blue color plane of additive television process (shadow-mask 
kinescope C = 500), and subtractive M.P. process (C = 71) broken lines. 2 

The relative luminance scale in Figures 5-7 is distorted, because a 
better appreciation of the contrast range and color purity of a repro- 
ducer is obtained by using the psychophysical Munsell liyhtness scale, 
V, which divides the dynamic luminance range of the visual system into 
uniform lightness steps. (The conversion from per cent luminance to 
lightness, V, is given in Table II). 

A "perfect" color reproduction requires an absolute black level 
(C = °o) and reproduces all spectrum colors with 100 per cent purity; 
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i.e., the color triangle is replaced by the spectrum locus. The lower 
boundaries of "perfect" color planes are therefore rectangles, indicated 
by C = «) and S = 1 up to the relative luminance for the spectral color 
point as exemplified by the red television primary. The upper boundary 
of the "perfect" color plane is determined by Equations (8) and (9) 
with l/C = 0 and w c equal to the tristimulus energy factor for the 
spectral color point. Because their form is not strongly dependent on 
the contrast, C, the upper boundaries of the television color planes 
R, G, B, and Y approach those of the "perfect" color planes for the 
corresponding spectral color points. It is seen from Figure 4 that the 
cyan and magenta regions of the color triangle lie approximately half 
way between the spectrum locus and the white point. The correspond- 
ing color planes in Figures 5 and 6 have, therefore, approximately one 
half the width of the corresponding "perfect" color planes. It is fur- 
ther evident that the 600-to-l contrast range obtained with a color 
kinescope comes much closer to a "perfect" color reproduction than 
does a motion picture. The color space boundaries of a color motion 
picture (broken lines taken from Reference (2) ) are, by comparison, 
much more restricted and illustrate the basic fact that a subtractive 
process rapidly loses saturation at higher luminance values, because 
an increase in film transmission reduces the dye concentration and its 
color filter action. 

To approach the chroma obtainable by the addititive color television 
reproduction, color film could be modified to have a neutral density 
range of approximately 3.7, to be used with a minimum neutral high- 
light density near unity to retain a sufficiently high dye concentration. 
Unfortunately this modification leads to the requirement for ten times 
more light from the projector which results in heating and other film 
problems. 

The difference in performance between normal color film and color 
television can be demonstrated strikingly by color photographs of tele- 
vision images having constant luminance (constant Y-signal) and 
various degrees of color saturation, obtained by a progressive increase 
of / and Q signals (chroma control). Photographs were taken with 
chroma signals 1, 2, and 3 times normal. The 2 times chroma value 
produced very high color saturation on the kinescope. The correspond- 
ing color photographs however taken with normal exposures show only 
minor increases in chroma as illustrated by plates I and II and expected 
from the color diagrams Figures 5-7, while the denser photographs 
reproduced in Plates III and IV taken with y 2 normal exposure give 
a better reproduction of the kinescope chroma at the expense of a 
shorter distorted contrast range. Relative chroma values are fairly 
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well reproduced in the prints, although the color purity is lower than 
in the transparencies and considerably lower than in the original 
kinescope image. 

Sine-Wave Spectra and Optical Passbands 

The transmission of color requires three independent video signals 
as compared to one for a monochrome image. Equal definition in a 
color image requires thus, in theory, a transmission system having 
3 times the information capacity of a monochrome system. An ap- 
praisal of the total information capacity of the NTSC color system 
can be obtained by comparing its information capacity with that of a 
color system having 3 equal independent channels, taking into account 
a number of non-ideal conditions arising in practical systems and 
deserving particular attention. 

The electrical frequency spectrum of a stationary monochrome 
television picture signal is a line spectrum of discrete frequency com- 
ponents which are harmonics of the frame frequency (30 cycles). Be- 
cause of this fact, a second line spectrum, the color signals (I, Q), 
can be added to the monochrome or luminance signal (Y) by inter- 
leaving its frequency components with those of the Y spectrum. Inter- 
leaving is accomplished by modulation of a color carrier frequency 
(3.579545 megacycles) which is made an odd multiple of the half-line 
and half-frame frequencies. To permit separation of color signals by 
synchronous demodulators in the receiver, one color signal (Q) is 
transmitted with double sidebands. It is, therefore, limited to a 600- 
kilocycle bandwidth, and a 600-kilocycle filter is required after de- 
modulation to eliminate all higher cross-product frequencies. The other 
signal, /, (also double-sideband up to 600 kilocycles) can have its band- 
width extended by single-sideband transmission. It is restricted to x / 2 
the color carrier frequency, i.e., a bandwidth of 1.8 megacycles, and a 
1.8-megacycle filter is required after demodulation to eliminate cross- 
talk. 

A perfect separation of the interleaved luminance and chrominance 
signals can be achieved with interleaved comb filters when the image is 
stationary. The inexpensive continuous passband filters used in prac- 
tical receivers, however, give rise to spurious color signals upon syn- 
chronous demodulation in the / and Q channels, caused by high-fre- 
quency F-signal components, and the chrominance signal (modulated 
color carrier) produces periodic errors in the luminance signal. These 
errors change polarity in successive frames and would cancel out in 
a linear system when integrated (by the eye) over two frame periods 
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(1/15 second). Practical systems, however, are not linear (all kines- 
copes are rectifiers), and integration by the eye is incomplete, par- 
ticularly for bright pictures. The errors, therefore, do not cancel 
completely even in stationary pictures. They become larger for moving 
objects, and attain full magnitude for random signals such as noise. 

The high-frequency "crosstalk" from the common 2-4 megacycle 
region of a "flat" Y spectrum into the color demodulator band causes 
orange and blue color tinting of horizontal monochrome detail (color- 
ing of resolution wedges in a monochrome test pattern) , and normally 
fine-grained camera noise in this region of the Y channel is heterodyned 




FREQUENCY- MEGACYCLES 

Pig. 8 — Filter response and "crosstalk" in a color system. 

into objectionable coarse color-noise (streaks) by the demolulators, as 
indicated by the large values of the broken-line cross product curves 
in Figure 8b. 

These undesirable effects can be substantially eliminated by intro- 
ducing a bandwidth limitation of 3.6 megacycles and a gradual "roll- 
off" into the Y channel before synthesis of the NTSC signal in the 
colorplexer of the transmitting station as shown by the solid-line curves 
in Figures 8a and 8b. Subsequent reemphasis (aperture correction) 
of the Y channel in the receiver (see Figure 8c) after removal of the 
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color carrier restores a good Y-signal response. The various degrees 
of color crosstalk are illustrated in Plates V and VI. (The bluish back- 
ground color of the "white" test pattern photographs resulted from the 
accidental omission of the ultraviolet filter on the camera lens.) The 
striking reduction of noise obtained by the bandwidth limitation and 
"roll off" is not reproduced because of integration by the time exposure. 
Considerable improvement in noise crosstalk is obtained by the band- 
width limitation alone. 

The complete elimination of the chrominance signal (modulated 
color carrier) from the Y signal (in the receiver) is not possible with 
continuous filters. It is therefore common practice to suppress the 
color carrier and its lowest sideband frequencies in the receiver by 
insertion of a filter (trap) which limits the Y-channel response to 3.6 
megacycles (see Figure 8c). The carrier interference is thus com- 
pletely eliminated in large areas, leaving only beat patterns of reduced 
amplitude (occurring near sharp vertical edges) from the remaining 
sideband components (7 amplitude in Y channel Figure 8c) which 
generally contain little energy and permit aperture correction of the 
remaining Y channel at the receiver.* 

There are three electrical passbands available for the transmission 
of color signals; these correspond to three two-dimensional optical 
passbands for the color image. Because of the rectilinear scanning 
process, the three optical passbands in the "vertical" (v) coordi- 
nate are alike. They are determined by the raster line number, and 
have the theoretical equivalent passbands: 3 N e(Y)v = N eU)v = N e(Q)v = 
490 lines. The equivalent passbands in the horizontal coordinate (h) 
are unequal and have the theoretical values N e(Y)h = 320 for A/ = 4 
megacycles, 2V e(m = 144 for A/ = 1.8 megacycles, and N e(Q)h = 48 lines 
for A/ = 0.6 megacycle. 

The equivalent symmetrical passbands 3 

_ / 4 \V2 

N e = l—N elv) N eih) \ (10) 

are hence N elT) = 446 lines, 

AT e(7) = 300 lines, 
iV e(0) = 173 lines. 



*The aperture correction approximately doubles the chrominance signal 
amplitude in the Y channel. 

3 Otto H. Schade, "Image Gradation, Graininess and Sharpness in 
television and Motion-Picture Systems, Part IV— Image Analysis in Photo- 
graphic and Television Systems," Jour. S.M.P.T.E., Vol. 64, p. 593, Novem- 
ber, 1955. 
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FREQUENCY- MEGACYCLES 



£t 0.463 1.25 3.25 




Fig. 9 — Sine-wave spectra of horizontal (F, I, and Q) transmission 
passbands used in the analysis. 

Their sum S N" e = l ines is hence 68 per cent of the sum of three 
equal passbands of 446 lines each. This simple appraisal of total 
information capacity in such a color system assumes theoretical rec- 
tangular frequency spectra having abrupt cutoff which are neither 
practical nor desirable for image transmission (monochrome or color) 
because of strong edge transients and spurious signals generated in 
the signal-separating process. 

The sine-wave response characteristics of practical electrical pass- 
bands used in the following analysis are shown in Figure 9. Their 
equivalent spatial horizontal passbands are N e{V)ll = 260," N F(tih = 100 
and N e{(/)ll = 37. The spatial vertical passband is determined by the 
cascaded value of the camera and kinescope sine-wave spectra and is 
in the order of N elr> = 200 for each of the three signals. The equiva- 
lent spectrum spaces for the three practical channels are the rectangu- 
lar solids illustrated In Figure 10. The equivalent symmetrical pass- 




CHANNEL Y I Q 

N e (v) 200 200 200 

N. (h) 260 100 37 



Fig. 10 — Two dimensional sine-wave spectra (Y, I, and Q passbands). 



* This Y channel is somewhat wider than for receiver use (Figure 8) 
because it does not contain a trap circuit for the carrier frequency. 
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bands obtained with Equation (7) form cylindrical spectrum spaces 
having the radii 

7t HY) - ■ 257 lines, 
W cU) - 160 lines, 
W C{Q) — 97 lines. 

Because of the coordinate transformation or cross-mixing processes 
and unequal passband limitations, the electrical sine-wave spectra 
("frequency response") for different colors are not equal in this system 
and can be determined as follows. Referring to Figure 1, it is seen 
that the electrical frequency spectra of the color signals, E R , E a , E B 
entering Matrix I are alike. The linear cross-mixing process to Y, I, 
and Q signals and back to E' R E' G E' n in the inverse Matrix II does not 
disturb the signal ratios in the range of 600 kilocycles where the 
frequency response in the Y, I, and Q channels is alike, i.e., E' R = E R , 
E' G = E G , E' B = E B . This range includes the complete electrical fre- 
quency spectrum required for transmission of the vertical spatial pass- 
bands of the three color functions which remain, therefore, unaffected 
by the cross-mixing processes. The horizontal passbands are normal in 
the 600-kilocycle range, i.e., a green signal E a , for example, will result 
in a green signal E' G on the green kinescope gun and in zero signals 
on the other two guns as indicated in Table I under the column A/ = 
0 — » 0.6. In the range from 600 kilocycles to 1.8 megacycles, however, 
all Q coefficients in Matrix II are zero because of the Q-filter cutoff. The 
matrix is no longer the inverse of Matrix I, and causes signal voltages 
to appear on all three kinescope grids as shown in Table I. Beyond 
1.8 megacycles, both Q and / coefficients are zero, with the result that 
all three kinescope guns receive equal signals. The last three columns 
of Table I show that the total luminance would remain constant in all 
sections of the total passband for a hypothetical linear kinescope, al- 
though the color or color mixture does not remain constant because of 
"spurious" signals. The complete video frequency spectra (at the 
control grids of the kinescope guns) have been plotted in Figures 11 
and 12 for the six colors listed in the Table I. Note that all of them 
are different, and some of the spurious responses have negative lobes 
(negative signs of coefficients indicate a phase reversal). 

Effect of Unequal Passbands on Color Detail 

It is a rather widely accepted opinion that the NTSC color trans- 
mission provides a three-color presentation for large areas, a two-color 
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presentation for medium sized areas and a presentation of fine detail 
without color information. This cannot be true because full color 
information is transmitted in the vertical dimension up to the finest 
detail which is not affected by the inequality of the electrical passbands. 
Detail color information is also transmitted in the horizontal dimen- 
sion, because only the true or "fundamental" color signal components 
have a normal frequency spectrum including a d-c component, while 
the frequency spectra of the "spurious" color components contain only 
horizontal a-c components and represent one-dimensional high-pass 
filters. It is improper to disregard the difference in d-c components and 
conclude from Figure 11, for example, that a vertical line group in 
a green sine-wave test pattern (for which the blue and red camera 



° 0 
o 

UJ 

m 

i? 0.5 



e(b') green (eg'J 




100 


200 


300 




RED 


(ErO 




3 




\ log/ 


200 


300 




LINE-NUMBER— N 4 



Pig. 11 — Horizontal sine-wave spectrum components at kinescope grids for 
green, red, and blue camera signals. 



signals are zero) would appear in black and white in the fine detail 
"mixed-high-frequency" region of line numbers N ]t > 200. 

Actually these lines appear green, because the green a-c signal only 
is raised completely above the black level by a d-c component, while the 
spurious color signals (E\ R) and E\ R) ) have an electrical zero-level 
axis, are rectified by the kinescope, and can only produce light of much 
lower intensity during positive half cycles. In computations of color 
mixtures for small areas it must be remembered that the normal Fourier 
relations between impulse forms, unit functions, and their spectra do 
not hold in systems containing nonlinear elements such as a kinescope. 
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To obtain the correct answer, it is necessary to convert frequencies 
to intensity functions (waveforms) at the input terminals of nonlinear 
elements such as the kinescope guns, project the waveforms over the 
nonlinear transfer characteristics, pass them through the low-pass 
filters representing the gun and optical performance, and then combine 
the light signals to color mixtures as indicated in Figure 1. The gen- 
eral effect of the kinescope "spot" size on fundamental color signals is 
a reduction of high-frequency components as indicated by broken line 
curves in Figure 11 for small signals. The spurious color signals, 




O 100 200 300 

LINE-NUMBER — N 4 



Fig. 12 — Horizontal sine-wave spectrum components at kinescope grids for 
yellow, magenta, and cyan camera signals. 



however, are rectified by the kinescope, so that each sine-wave fre- 
quency is replaced by a series of even-order harmonics which is atten- 
uated more rapidly by the kinescope low-pass spectrum than indicated 
by the dotted lines, which lose their meaning as a sine-wave response. 

The observed transfer characteristic of the actual kinescope is 
substantially a square-law characteristic. The intensity function of 
the fundamental green sine-wave pattern on the kinescope screen 
(neglecting kinescope spot size) can hence be calculated by squaring 
the instantaneous sine-wave grid-signal values measured from the 
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current cutoff point, as illustrated in Figure 13 (top) for two fre- 
quencies (points Pj and P 2 in Figure 11) . Because of the missing d-c 
component, the squared spurious red and blue a-c signals become very 
small and only their positive half-waves produce light, which adds to 
the green sine-wave light, on the kinescope screen, diluting the chroma 
as indicated. Similar conditions prevail for test patterns of other 
saturated colors (see Figure 13, bottom) as summarized in Table III. 

The fine-line groups in the "mixed-high-frequency" region, there- 
fore, retain substantially the color of the fundamental signal, although 
they acquire a spurious tint. Their relative amplitude and luminance 
is determined largely by the amplitude of the fundamental color com- 
ponent which is proportional to the Y component of the particular 
color and therefore lowest for blue (11.4 per cent). These findings are 
confirmed visually by turning off one or two of the three equal color 
signals generated by a "white" SMPTE test pattern ahead of Matrix I 
(see Figure 1) and comparing the color in the vertical line-wedge 
reproduction with the pure color reproduction in the horizontal line- 
wedges of the kinescope image.* 

A similar result is obtaind by an analysis of line images. The 
images of fine horizontal lines are transmitted without color change 
(because of the matched vertical sine-wave response in all channels) . 
The color of fine vertical lines, however, is affected by the unequal 
horizontal-frequency spectra and is determined by the sums of the 
three squared impulse components. The impulse forms computed for 
unit signals in the Y, I, and Q channels are shown in Figure 14, cor- 
responding to the frequency spectra shown in Figure 9. For a given 
color-line signal the impulse forms at each of the three kinescope con- 
trol grids are sums of these three impulses modified in amplitude by 
the amplitude factors for the particular color signal mixture (Table I, 
column group 2). The squared positive amplitudes of these impulse 
sums are transduced by the kinescope into light, i.e., into line-images, 
which are then added and analyzed for color mixture. 

The impulse forms at the control grids of a color kinescope are 
shown in the top row of Figure 15a for a green line object. Note that 
only the fundamental color (green) has a positive average value (the 
spurious signal impulses red and blue contain no d-c term) . A low- 
intensity green line object superimposed on a white background would 
represent a small signal on a d-c level. In this case the spurious as 
well as the fundamental signal components are transduced linearly and 



" It is apparent that overemphasis of the electrical high-frequency 
response by excessive aperture correction increases the spurious color signals 
and chromatic errors. 
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Fig. 13 — Intensity function and components (broken lines) of green sine- 
wave images on the kinescope screen (top), and of red sine-wave images 

(bottom) . 

without rectification by the kinescope, resulting in the impulse sum 
shown at left in Figure 15(c). A high-intensity line object against 
a black background is a large signal on a zero d-c level, requiring sepa- 
rate squaring of the gun signals to obtain the intensity functions of 
the three colored line images on the kinescope screen shown in Figure 
15b, which have no negative amplitudes (clipped by the kinescope) . The 
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Fig. 14 — Unit-impulse forms of Y, I, and Q channel from Matrix I 
and of kinescope. 
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sum of these optical intensity functions is shown at the right side in 
Figure 15c. Since the transfer gain is adjusted so that equal electrical 
amplitudes result in "white," the amount of white contained in the 
impulse sum is given by multiplying the smallest one of three common 
components (R in this case) by three at any one distance from the line- 
image center. It is seen that the white dilution is lower in a squared 
large-signal line image than in a small signal, which by itself is not a 
pure color. This condition is also obtained for other colors as shown 
in Figures 1G and 17. It follows again that vertical line images retain 
substantially their original hue somewhat tinted, desaturated, or re- 
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Fig. 15 — Relative intensity of impulse components and sums forming the 

green line image. 

duced in amplitude as summarizel in Table IV and substantiated visu- 
ally in the color photographs Plate VII by comparison with the normal 
color reproduction of horizontal line objects. The illustration shows 
small sections of 5*4 x 7-inch picture reproductions. The faint blue 
tint in the vertical green lines is barely visible. 

Color transitions at sharp edges can be synthesized in the same 
manner from the component step functions shown in Figure 18. The 
results are plotted in Figures 19-23. Inspection shows that both color 
and waveform distortions vary with color and are again relatively 
smaller in high-contrast transitions than in low-contrast transitions 
because of the square-law kinescope characteristics. 

It should be pointed out that the effect of camera tube and kinescope 
beam sizes as well as the %-power gamma correction for the kinescope 
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have not been included in the line image and step-function evaluations. 
The gamma-correction amplifier itself would not change an ideal unit 
impulse or a step function, but it does round off the low-level corner 
of the actual S-shaped step function from the camera relative to its 
high-level corner, thus making the large-signal step functions more 
symmetrical and all step functions less steep with reduced transient 
ripples. To include the effect of the kinescope beam size, it is necessary 
to perform a convolution of its line image or its step-function response 
with the computed response forms. The finite kinescope beam size 

SMALL SIGNALS LARGE SIGNALS 

(linear) (rectified AND SQUARED) 




Pig. 16 — Relative intensity of red and blue line images. 

broadens the line images and transitions somewhat more, and further 
reduces high-frequency transients. The visibility of the luminance 
errors caused by the unequal passbands of the NTSC system can be 
assessed by computing the visual lightness of horizontal step functions 
and line images for red, green, and blue as shown in Figure 24. Con- 
stant-bandwidth vertical step functions and line images are shown for 
comparison, and indicate that the errors in the blue are readily visible 
as confirmed by observation (color plate VII). 

It can thus be concluded that considerable color information is 
transmitted and reproduced even in fine horizontal detail with NTSC 
color signals and that a high gamma (square-law) kinescope transfer 
characteristic is a definite asset because it reduces the color distortion 
introduced by unequal color-transmission channels. Furthermore, the 
fact that relatively small color changes can be seen in the fine detail 
of color television images raises doubts that the eye's color mechanism 
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does not function when observing fine detail, and suggests a reappraisal 
of the capabilities of the visual system at normal luminance values. 
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Fig. 17 — Relative intensity of yellow, magenta, and cyan line images. 



PART II— SOME CHARACTERISTICS OF VISUAL SYSTEM AND 
EFFECT ON OVER-ALL SYSTEM CONCEPT 



Transfer Functions, Wavelength Functions, and Incremental 
Sensitivity for Luminance and Color 

The visual system is a very complex system which terminates in a 
"computer" (the brain) having random connections, and very little is 
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Fig. 18 Unit step-function components from Matrix I and of kinescope. 

known about its memory and interpreting processes.* There is ample 
evidence that the system is nonlinear in many of its sections. It is 
well known that both hue and lightness of colored objects are functions 
of the surrounding background, although there is no spectral or colori- 
metric change. It is also known that the apparent "lightness" of a 
color mixture does not necessarily agree with the luminance value 
computed as the linear sum of component luminances. (A bluish-white 




Fig. 19 — Relative intensity of step-function components and sums forming 
the green step function image. 



* An interesting theory on this subject is discussed in Reference (4). 
4 Frank Rosenblatt, "The Design of an Intelligent Automaton," Re- 
search Trends, Cornell Aeronautical Lab., Inc., Vol. VI, No. 2, Summer, 1958. 
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Pig. 20 — Relative intensity of step-function components and sums forming 
the red step function image. 



screen, for example, appears brighter than a screen of equal luminance 
illuminated by a low-color temperature light source.) Observations 
made by the author on equiluminous sine-wave mixtures indicate, like- 
wise, that the fixed relation of lightness and luminance given in Table 
II does not necessarily apply to incremental amounts of the components 
in a mixture of primary colors, but that it is strongly dependent on the 
relative total amounts of the primaries. 

With regards to chromatic sensitivity it is known that the just- 




Fig. 21 — Relative intensity of step-function components and sums forming 
the blue step function image. 
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perceptible amounts of pure color added to white is of the order of 
1 per cent for red and blue and about 2 per cent for yellow and green. 5 
The chromatic sensitivity //A/ or A/AA. of the eye for spectral colors" 
has two peaks (less than 3:1) at yellow and cyan above a substantially 
uniform level from red (650 millimicrons) to blue (430 millimicrons). 
From the system analysis point of view, it does not appear likely that 
the wavelength functions of the eye's color discriminating system are 





Fig. 22 — Relative intensity of yellow, magenta and cyan 
step-function images. 

as simply related to the wavelength function of its luminance channel 
as they are in a photographic or color television system, because the 
eye's luminance system continues to operate at light intensities far 
below the threshold for color vision. 

It has been demonstrated that the eye can see all of the calculated 
defects in the color detail reproduction of the NTSC system, including 
those in low luminance colors, such as blue, for which the eye's "acuity" 

5 J. W. T. Walsh, Photometry, Constable, London, England, 1926, p. 67. 

6 L. A. Jones, "The Fundamental Scale of Pure Hue and Retinal Sensi- 
bility to Hue Differences," Jour. Opt. Soc. Amer., Vol. 1, p. 63, March-May, 
1017; W. D. Wright, Researches on Normal and Defective Color Vision, 
Kimpton, London, England, 1946. 
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Fig. 23 — Relative intensity of green to magenta step-function image. 

is low. It is important to understand clearly that acuity and resolving 
power are not determined by the spatial sine-wave response alone, be- 
cause they are strongly dependent on the noise level as well. The 
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Fig. 24 — Step functions and line images (NTSC) in munsell lightness 
values. Solid curves horizontal, broken line curves vertical coordinate. 
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remarkable visual constancy of color mixtures over a large part of 
the photopic range may be taken as an indication that the transfer 
characteristics, spatial sine-wave response characteristics, and noise 
levels of the eye's color discriminating system are well matched. Sim- 
ilar noise levels in turn indicate that the "discriminators" receive 
similar signal levels, or more likely, similar signal differences from 
the photo receptors in at least three principal regions of the visual spec- 
trum. Their spatial sine-wave response appears to be well matched, 
because there is no evidence to the contrary in their performance (no 
color fringes, transients, or "smears"). 

It was therefore decided to determine the optical passbands of the 
visual color discriminators by sine-wave response measurements, fol- 
lowing a procedure previously described for the luminance channel of 
the eye. 7 

Sine-Wave Response Functions 

There have been general observations to the effect that resolution 
and sharpness of images (in microscopy, for example) do not change 
appreciably with blue, green, or red illumination, and that the reada- 
bility of colored print of equal luminance is fairly independent of color 
with red slightly better than green or blue. These conclusions are not 
convincing, because they relate to an interpretation of signal combina- 
tions (depending on noise level and sine-wave response) from the 
luminance as well as the color mechanism. The same criticism can in 
part be applied to the single-color sine-wave response characteristics, 
Figure 25, measured by the author on an observer having normal 
vision. The response of a luminance system can, however, be eliminated 
by the use of equi-luminous sine-wave objects, which provide a sinu- 
soidal variation of excitation purity at constant luminance. A red 
purity change, for example, is represented in the red color plane (Fig- 
ure 5) by a sinusoidal modulation of the radius length (horizontal 
distance) , which can be centered near the white axis (near-white back- 
ground) or may be moved out towards a higher purity point. The 
generation of such sine-wave patterns on a color kinescope screen is 
relatively simple, as discussed in the Appendix. 

Sine-wave patterns of constant luminance were set up by measuring 
the sine-wave luminance distribution with a luminance meter,* and 

7 Otto H. Schadc, "Optical and Photoelectric Analog cf the Eye," Jour. 
Opt. Soc. Amur., Vol. 46, p. 721, September, 195G. 

* Weston Foot-Lambert Meter, New Candle Model 931 with Viscor 
Filter, corrected for visual response. The meter is designed for kinescope 
measurements and gave correct luminance readings on the phosphor pri- 
maries and their mixtures. The readings agreed with comparative readings 
taken with a Macbeth Illuminometer. 
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adjusting the component intensity for a constant reading. It was ob- 
served immediately that luminance does not have a constant relation 
to visual brightness. The adjustment for constant luminance (by 
calculation or meter) differed from the incremental mixture values 
required for minimum flicker and constant visual brightness, particu- 
larly for equiluminous blue purity sine-waves in a near white back- 
ground where the increment-component ratio was in error by several 
hundred per cent. The flicker test for equal brightness was made by 
reducing the temporal (electrical) sine-wave modulation frequency to 
20 or 25 cycles and adjusting, for example, the blue component (see 
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Fig. 25 — Response of visual system to single-color sine-wave patterns. 
Appendix) for minimum flicker. Because of the possibility of errors 
arising from the difference in phosphor decay time for different 
kinescope primaries, the brightness equality of the positive and nega- 
tive purity half waves was judged visually at a higher modulation 
frequency producing a wavelength of about 2 centimeters on the 
kinescope screen. It was found to agree with the flicker test. When 
such a stationary wave pattern is observed while one component is 
varied, a sudden phase shift seems to occur in the sine-wave position 
when constant brightness is reached. 

The contant-brightness condition near the white point requires 
more nearly* equal stimulus energy than equiluminance of the com- 



This might be expected from a discriminator near its neutral point. 
A more extended quantitative study of incremental luminance and bright- 
ness ratios in sine-wave objects would be of value. The authcr is aware 
of many publications (in the Journal of the Optical Society of America 
and elsewhere) discussing discrepancies in the total luminance of a color 
from standard sums. 
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ponents as most strikingly evident from blue purity modulation tests. 
It approaches equiluminance when the operating point (background 
color) is moved towards the spectrum locus, which is to be expected 
because the luminosity curve is determined from spectrum colors by 
brightness equality tests. This observation may well explain observed 
discrepancies in color noise visibility with the NTSC color system, 
because equal luminance increments do not guarantee equal increments 
in visual brightness. 

At the time of this writing only one set of sine-wave response- 
measurements taken near the white axis (B w = 7 foot-lamberts) is 
available. The designations red, green, blue, in Figure 26, refer to 
the color plane in which the purity modulation occurs, although the 
color in the areas of reduced "negative" purity (negative half cycles 
of sine wave) appears to the eye as a complementary color in com- 
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Pig. 26 — Response of visual color discriminator system to constant 
brightness sine-wave patterns (purity modulation). 

parison with the color of adjacent areas of higher "positive" purity, 
even at a larger distance from the white point. (The negative blue bar 
appears yellow, negative green appears magenta, etc.) Lacking a 
luminance change, the bars appear to have a constant purity as if 
produced by a square-wave rather than a sine-wave modulation. 

The measured sine-wave response functions of the color discrim- 
inator elements in the visual system differ from one another by much 
smaller factors than the "standard" luminance values. All show a loss 
of low-frequency and d-c response as in the luminance channel. An 
analog of the visual system 7 contains, therefore, a high-pass filter or 
a negative frequency-limited feedback loop which superimposes a nega- 
tive image of low definition on the normal image, "inhibiting" the 
normal low-pass filter response by addition of a negative response as 
reconstructed in Figure 27 (see Appendix). The feedback is least for 
red, and largest for green light, with blue intermediate. More data are 
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required to support a more quantitative analysis of the details of these 
results. 



Equivalent Passbands and Number of Receptors 
in Effective Retinal Sampling Areas 

The equivalent passbands, N e( + ) , of the main low-pass character- 
istics give information on the effective discriminator areas in the retina 
and appear to be in general agreement with Polyak's anatomically well- 
founded deductions 8 that the red end of the color spectrum is signaled 
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Pig. 27 — Positive and negative sine-wave spectrum components of visual 
color discriminator system. 

over direct "private line" connections between optic nerve fibers and 
cones and that it has the best detail response, while green and blue 
signals involve more diffuse matrix networks in the order stated. 

The effective sampling area a (area of equivalent constant intensity 
point-image) can be computed from the equivalent passband by 
a = 1.16 >'N e - (see Reference (3)). It has been shown in Reference (7) 
that at a luminance of 7 foot-lamberts the equivalent passband of the 
eye's lens is 2V B(M = 50 television lines per millimeter on the retina. 
The effective sampling area a L of the optical point image produced by 
the lens of the eye or the retina includes therefore 88 2.2-micron cones 
at this illumination as illustrated in Figure 28. The normal positive 



111 1941 L ' P °' yak ' M ' D '' The Ret " M > University of Chicago Press, Chicago, 



QUALITY OF COLOR-TELEVISION IMAGES 



527 



sine- wave response function (Figure 27) of the visual system is the 
product of the sine-wave response function of the lens and the sine- 
wave response function of the color discriminator neuron system. 
Since the equivalent area of the final point-image in a cascaded system 
is simply the sum of the equivalent areas of its cascaded components, 
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Fig. 28 — Diameter S, (microns) of equivalent round sampling area of 
optical point-image and color discriminators of visual system relative to 
the diameter (2.2/i) of the photo-receptors (cones). 

the equivalent areas, a,,, of the color discriminator system are obtained 
as the difference 

«(i = »m — o-l, 



where a d = effective area of discriminator neuron system, a m = 1.16/ 
(N ei + )i /4.25) 2 = equivalent area of entire system,* and a L = 1.1Q/N ch 2 
= 1.16/2500 square millimeters = equivalent sampling area of lens. 
The diameters obtained for equivalent circular sampling areas a d of 



* The equivalent passband Nh'h of the over-all sine-wave function 
(Figure 27) is measured at an object distance equal to four length units 
(indicated by the subscript 4) ; since the effective distance to the retinal 
image is 17 millimeters (effective focal length of lens) this length unit cor- 
responds to 17/4 = 4.25 millimeters on the retina, requiring division by 
4.25 to obtain lines per millimeter cn the retina. 
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the color discriminator neuron systems are listed in Figure 28.* The 
red discriminator receives signal contributions from approximately 12 
cones, while the more diffuse green and blue discriminators collect 
signals from equivalent areas containing about 39 and 59 cones re- 
spectively. The intensity distribution in the actual point-image of the 
lens and the chromatic sampling areas is approximately gaussian as 
illustrated in Figure 28, and the actual areas include a larger number 
of cones with partial signal contributions. The relatively large number 
of cones associated with a single point-image area does give a logical 
explanation to the known fact that the retina contains a much smaller 
number of matrix units (bipolar and ganglion cells) than cones, and 
that this ratio is adequate to derive color signals from the receptors, 
luith only minor effects on the over-all passbands. 

System Concepts Unkesolved 

The complexity of the visual system, much of it unexplained, leaves 
room for much variety with regards to its possible functioning as a 
color system, and many different theories can be found in the litera- 
ture.** The number of possible systems, however, is reduced by the 
high "quantum efficiency" of the eye 7 ' 11 which does not seem to agree 
with the systems of area-sharing receptors, where each covers only a 
narrow region of the visible spectrum. Given supersensitive receptors, 
even slight dissimilarities in three spectral response characteristics 
similar to a luminosity function are quite sufficient to produce excellent 
color discriminators by signal subtraction and to retain full luminance 
sensitivity by addition. Like the eye, such a system will continue to 
"see" images in monochrome when the scene illumination becomes 
insufficient for operation of its discriminators. t 

The nonlinearity and the matrix system of the retina allow for 



* The signal integration computed from the measured sine-wave spectra 
agrees in order of magnitude with the number of synaptic connections from 
optic nerves over specialized types of ganglion and bipolar cells to the 
cones in Reference (8). 

** Excellent discussions are given in References (9) and (10). 

"f At still lower illumination (scctopic range), cone vision ceases and 
only the rod system remains in operation. 

0 Gunnar Svaetichin and Rune Jonasson, I. A Technique for Oscillo- 
graphic Recording of Spectral Response Curves; II. Spectral Response 
Curves from Single Cones; III. Receptcr Mechanisms for Flicker and 
Fusion; IV. Notes on the ERG Analysis; V. The Cone Function Related 
to the Activity of Retinal Neurons; VI. Aspects on Human Photoreceptor 
Mechanisms," Acta Physiologica Scandinavica, Vol. 39, Supplementum 134, 
Stockholm, 1956. 

10 S. A. Talbot, "Recent Concepts of Retinal Color Mechanism," Jour. 
Opt. Soc. Amer., Vol. 41, p. 895, December, 1951. 

11 A. Rose. "Quantum Effects in Human Vision," Advances in Biological 
and Medical Physics, Vol. V, Academic Press Inc., New York, 1957. 
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variation of amplitudes and weighting of respective gains (negative 
feedback) in combination with the luminance signal, so that the "com- 
puter" (the brain) assigns less importance to signals from dark- 
colored objects in a bright surrounding than it does to small-component 
signals in a color mixture, or when the general "lightness" of all 
colors in the viewing field is more balanced. A similar mechanism is 
effective for pure neutral-luminance signals, where lower-intensity 
signals are suppressed near strong ones, to enhance contrast by cre- 
ating a subjective black-level. In such cases one could replace the low- 
intensity surrounding of an object by a black surrounding without 
introduction of error, but this does obviously not permit the omission 
of all low-intensity signals. For the same reason one cannot replace 
blue, or for that matter any dark-appearing color, by a dark neutral 
or black except under specific conditions. 

Choice of Color Axes in External Systems 
Having Unequal Passbands 

It follows from the discussion of the visual system that, given three 
unequal passbands for an external color-reproducing system, a pre- 
ferred assignment to particular colors cannot be made on the basis of 
the eye's spatial frequency response, since there are apparently no 
major inequalities between color and luminance passbands in the visual 
system. One must look towards other characteristics, such as the 
chromatic aberration* of its lens, the lightness weighting, and the 
subjective black-level (feedback) mechanism, which, in a comparison, 
cause it to attach less importance to "dark" colors, thereby permitting 
errors in their reproduction to be noticed less frequently. 

In the NTSC system the translation of the camera primaries R, G, B 
into the transmission primaries Y, I, Q containing a luminance channel 
(Y) was essential to establish compatibility with monochrome re- 
ceivers. Matrix I also assigns a percentage of the full video passband 
to each color. In view of the above analysis it appears that a somewhat 
more panchromatic distribution of these percentages (higher blue con- 
tent in Y) could be of some advantage. 

Having assigned the normal wide passband to the luminance channel 
of the NTSC color system, there will be no inequality of optical pass- 
bands in the reproduction of black-and-white objects. The axes for the 
remaining two theoretically nonluminous" color signals (/, Q) must 

* The axial chrcmatic aberration ( Jour. Opt. Soc. Amer., Vol. 47, No. 
6, 1947) of the human eye is approximately 2 diopters from 400 to 550 
millimicrons and less than 1 diopter from 550-700 millimicrons. 

** Because of the nonlinearity in an actual system, the / and Q signals 
cause small luminance components, 



530 RCA REVIEW December 195S 

pass through the white point, and thus affect pairs of complementary- 
colors, such as red and blue-green, orange and cyan, yellow and blue, 
yellow-green and purple, green and magenta. A balanced psychological 
weighting condition might be indicated by letting the Q component be 
larger for high-luminance colors (which have a large wide-band Y com- 
ponent), but except for the red and blue-green color pair, all pairs 
contain one higher and one lower luminance color. The choice of I- and 
Q-axis directions is thus not very critical. The two color axes should 
include in general a large angle and the Q axis should not come too 
close to the low-luminance blue primary on one side, and point towards 
the high-luminance green region on the other side of the white point. 
These considerations place the I axis in the orange-cyan region, helping 
the medium-luminance red. The NTSC system axes have these locations 
(See Figure 4) . 

Conclusions 

The above analysis as well as visual tests indicate that the band- 
width restriction in the / and Q channels causes small red-colored 
objects to lose some intensity and blur slightly (relative to white) ; 
green- and yellow-colored objects are sharp but lose some saturation, 
while small blue-colored objects have good saturation but lose intensity 
and sharpness. Because these effects are smaller in high-luminance 
colors, disappear towards white, and occur in one (horizontal) direction 
only, the degradation in a normal two-dimensional color picture trans- 
mitted by the NTSC system is relatively small in comparison with a 
color transmission over three equal independent channels as illustrated 
by plates VIII and IX. The relatively uncritical nature of the band- 
width assignments can be appreciated by deliberately disturbing the 
luminance weighting by an interchange of red, blue, or green signal 
connections at the input to Matrix I, making the corresponding change 
after Matrix II, and observing the results visually at the kinescope 
as illustrated by the photographs plates IX to XI. The reversal of 
luminance weightings for green and blue by interchange of matrix 
connections (plate X) can be detected only by the appearance of a 
transient (following the right side edge of the hat). The reversal of 
luminance weightings for red and green by interchange of matrix 
connections (plate XI) give no detectable effect for this subject as 
compared to the normal reproduction in plate IX. Critical observation 
of test pattern slides reveals that the defects caused by unequal pass- 
bands are readily visible and change to different colors when the above 
interchanges are made, and that there is a preference towards the 
NTSC choice. In a direct demonstration on a kinescope, most observers 
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(engineers) were unable to recognize the equal-bandwidth condition 
or the interchange of matrix connections when shown a variety of 
outdoor and indoor color pictures (SMPTE test series and many 
others) . 

The analysis has shown that commercial color kinescopes have a 
color contrast range resulting in a color space which is larger than 
that of commercial color motion pictures, that color reproduction 
errors caused by the unequal passbands of the system in fine detail 
are relatively small, and that they are considerably reduced by rectifi- 
cation in the nonlinear kinescope. It has also been shown that the 
spurious signals generated in the practical separation of chrominance 
and luminance signals from the composite (NTSC) signal can be 
minimized by proper choice of amplitude response in the various band- 
limiting filters at the transmitter and receiver, and result in a rela- 
tively small loss of bandwidth and picture sharpness as illustrated by 
plates XII and XIII. 

The relative sharpness is well reproduced in the color plates. (Note 
the slight horizontal edge transient in the NTSC reproduction, plate 
XII.) The color saturation, however, is considerably reduced by the 
printing process. 

Regarding the performance of commercially available color re- 
ceivers, it can be stated that the contrast of the color tube is as good 
as that of the tube used in these tests. Pictures received from commer- 
cial color broadcasts can be, and on many occasions have been, as good 
as those observed in these tests. However, with color even more than 
with black-and-white reception, performance depends on proper adjust- 
ments of the receiver controls and, of course, the ambient light level. 
The fine detail monochrome performance should be close to the values 
given (again assuming correct adjustments) , but the horizontal color 
detail is not quite as good, because it is present practice to use a com- 
mon intermediate bandwidth for both / and Q channels in the receiver. 

Appendix 

Generation of Equi-Luminous and Constant-Brightness Sine Waves 

A white background is set up first by adjustment of d-c components. 
The color gun (or guns) producing the desired hue, red for example, is 
then modulated with a sine-wave signal to generate a vertical sine-wave 
field, and the guns producing a complementary signal (white minus 
red = blue plus green) are simultaneously modulated by negative sine- 
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Plates I-IV — Ektachrome performance in recording constant lumi- 
nance images. Plate I, normal chroma, normal exposure; Plate II, 2X 
chroma, normal exposure; Plate III, normal chroma at V2 exposure; 
Plate IV, 2X chroma at Y% exposure. 

Plates V, VI — Monochrome test pattern reproductions by NTSC 
system. Plate V, left, 4-megacycle flat Y-channel input, demodulators 
off (no crosstalk) ; Plate V, right, 4-megacycle flat Y-channel input, 
demodulators on; Plate VI, left, 4-megacycle gradual roll-off Y channel, 
demodulators on, subsequent re-emphasis ; Plate VI, right, 3-megacycle 
flat Y channel, demodulators on, normal Y response in receiver. 

Plate VII — Sample of three line-object reproductions, 4-megacycle 

Y channel and standard NTSC / and Q channels; pure color signals 
obtained from a white pattern by turning off unwanted color channels 
in camera. 

Plates VIII-XI — Effects of unequal passbands and luminance- 
weighting on kinescope image. Plate VIII, three independent 4-mega- 
cycle channels (R, G, B) ; Plate IX, Y + I + Q signals, standard NTSC 
matrix and filters as in Figure 1 ; Plate X, Y + I + Q signals and 
filters, but G and B reversed to .Matrix I and from Matrix II ; Plate XI, 

Y + I + Q signals and filters, but R and G reversed to Matrix I and 
from Matrix II. 

Plate XII — Complete NTSC signal with standard matrix and filters, 
including modulators, multiplexing, and demodulator circuits. 

Plate XIII — Three independent 4-megacycle channels (R, G, B) . 
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wave signals (opposite phase) of the same frequency.* To adjust the 
proper mixture, all guns are first modulated in phase to produce a 
"white" sine-wave, after which the red gun signal is reversed in phase 
and adjusted in amplitude for constant luminance or constant visual 
brightness of the pattern. It is observed that this adjustment varies 
with the background color (horizontal distance from the white-axis) 
which is, therefore, adjusted first by change of d-c components. 

The adjustment for constant brightness (red sine-w^ave amplitude) 
must be made for the actual observer to make his luminance mechanism 
inoperative, because the values (Y) computed from the CIE "standard 
observer" curve are in error near the white axis as mentioned. This 
error is surprisingly large for blue (7 different observers) . 

Measurement of Sine-Wave Response Functions and Interpretation of 
Data 

The observer views, from a fixed distance, a uniformly illuminated 
field on which an extended sine-wave bar pattern is faded in slowly by 
a master control ahead of the potential dividers for the component 
""" : i ' . The pattern becomes visible when the brain receives a just- 
ptible signal. A plot of the required optical input signal modula- 
tion is shown in Figure 29. The operating point in the color plane 
'iear the white point and had a total luminance indicated in Figure 
B w = 7 foot-lamberts. Since the modulation is given in relative 
of the color under test and relative stimulus values are close to 
lues w R : iv G : w n = 1:0.92:1.57 of the white point, the minimum 
tible stimulus increments of the three colors differ by less than 
t of two from one another at any spatial frequency. The normal- 
■( iprocal values of the sine-wave input signals plotted as a func- 
' spatial frequency, f„ or line number, N = 2f s , furnish the sine- 
implitude response of the visual system (Figure 26) provided 
oshold criterion is independent of frequency. Otherwise a cor- 
1 is needed. Upon examination one finds that the visual threshold 
limited and that detection of any object area requires a certain 
signal-to-noise ratio. Because the area of sine-wave bars of 
length and one half cycle width decreases in inverse proportion 
line number, signal-to-noise ratio, R 0 , and sensitivity should 
ortional to l/\/N, as illustrated by curve 1 in Figure 30. It is 
nowever that the eye must vibrate (ocular tremor) in order to 

noninterlaced raster of 500 lines is used to obtain a di"ect propor- 
of optical and electrical frequencies (see Reference (12)). 
tto H. Schade, "A Method of Measuring the Optical Sine-Wave 
Spectrum of Television Image Display Devices," Jour. S.M.P.T.E., 
1 , p. 561, September, 1958. 
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COMPLEMENTARY COLOR 




Fig. 29 — Optical modulation intensities required for threshold vision 
(corrected for kinescope response). 
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Fig. 30 — Effect of signal-to-noise ratio, ocular tremor, and transient 
motion on threshold sensitivity. 
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see detail; i.e., it must generate a-c signals for transmission to the 
brain (its analog contains a blocking capacitor). The deflection ampli- 
tude of the tremor is quite small, 6 microns on the average, 13 which 
does not degrade the image, but is sufficient to generate strong a-c 
signals at line numbers where the half cycle length, 1/N, of the retinal 
sine-wave image is in the order of the deflection amplitude. The a-c 
signal, however, decreases towards lower line numbers and goes to 
zero at N — 0. The average a-c signal developed at different points 
on the retina due to ocular tremor is, therefore, a function of line 
number, and calculations indicate the function shown by curve 2 in 
Figure 30. The deflection amplitude is specified by its reciprocal value, 
expressed as a television line number N TW , for an object distance equal 
to four length units.* 

The action of the tremor counteracts the sensitivity change due to 
the varying signal-to-noise ratio as shown by the product curve 3, but 
overcompensates at low frequencies. This deficiency is decreased by 
a partial d-c response 7 and by the continuous jerky motion of the 
eyeball when a larger image is observed. One becomes conscious of this 
motion and its effect in generating a-c signals when a barely visible 
large object or pattern is viewed, because the pattern fades out when 
fixed steadily. The signal generation in this case is a series of detached 
transients (a continuous large deflection would blur the image). A 
transient image appears after the sudden motion stops, just as in 
tests where ocular tremor and image motion are prevented artificially 
and transient vision is obtained after a change of image content. The 
dashed curves in Figure 30 show that a constant threshold requires only 
a relatively small increase in a-c signal at the low-frequency end of 
the response curve; this is obtainable by intermittent transient ex- 
citation. 

In view of these facts, it appears justified to assume a constant 
threshold criterion in the upper portion of the sine-wave response 
characteristic of the eye which is obtained directly from the measured 
data. The lower portion is less well defined, but indicates that an 
analog system contains a high-pass filter in cascade with the normal 
low-pass filter representing the optic and retinal structure of the eye. 
The low-frequency section of the low-pass filter response can be extra- 
polated by matching the high-frequency section to a possible lens plus 
grain structure response characteristic (see Figure 27). 



13 Lorrin A. Riggs, "The Measurement of Normal Ocular Tremor bv 
Corneal Reflection," Jour. Opt. Soc. Amer., Vol. 42, p. 287 April 1952- 
G. C. Higgins and K. F. Stults, "The Frequency and Amplitude of Ocular 
Tremor, Jour. Opt. Soc. Amer., Vol. 42, p. 872, November, 1952 

* The reciprocal distance on the retina in television lines/mm is eiven 
by N/mm = N (1) /4.25. 



THE EFFECT OF RADIATION ON SILICON SOLAR- 
ENERGY CONVERTERS*t 

By 

j. j. loferski and p. rappaport 

RCA Laboratories, 
Princeton, N. J. 

Summary — The performance of silicon solar cells under simultaneous 
illumination and irradiation by various ionizing radiations has been ob- 
served for the purpose of estimating the useful life of such power sources 
in the environment of the l.G.Y. (International Geophysical Year) earth 
satellite. Experiments were conducted both in air and in vacuum for 
2 mev electrons, 20 mev protons and 40 mev alpha particles. The maximum 
power, open-circuit voltage and short-circuit current experienced a decay. 
These results have been extrapolated to yield an estimated minimum of 
about 10 s years until the output of such cells on a satellite drops to 75 
per cent of the initial value. 

Possible modifications of these estimates based on tentative data ac- 
quired from satellite' flights have been considered. 

The decay of the cells is associated with changes in minority-carrier 
lifetime and deviations from expected behavior are discussed. No definite 
decay was observed during the limited exposure time of this study for the 
ultraviolet radiation between 2200a and 3400a or for x-rays from 
machines operated at 50-2,000 kilovolts. 



Introduction 



V j|"-^HE PURPOSE of this study was to investigate the effects of 
I cosmic and solar radiations likely to be encountered in the upper 
atmosphere on the performance of silicon p-n junction solar 
energy converters. Specifically, experimental study was made of the 
effects on such cells of ultraviolet light, x-rays, gamma rays, electrons, 
protons and alpha particles both in air and in a vacuum. Radiation 
densities comparable to, or greater than, those expected outside the 
earth's atmosphere were used to ascertain at what levels they would 
produce detrimental effects on the electrical characteristics of the 
cells, such as short-circuit current (i s ), maximum voltage (V mai ), 
maximum power (P niax ) , and maximum solar energy conversion effi- 
ciency (?? m . ix ). Experimental results are compared with theory. 

Throughout these experiments the cells were simultaneously illu- 
minated by a light source, whose power density approximated sunlight 



4 Manuscript received September 4, 1958. 

t This work supported by the United States Army Signal Research and 
Development Laboratory Contract #DA36-039-SC-64643. 
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(100 milliwatts per square centimeter), and irradiated by the poten- 
tially detrimental radiation. The irradiated cells were kept within the 
temperature range of 25°C to 30°C because it is anticipated that satel- 
lite temperatures will have an equilibrium value within 15 C C of the 
launching temperature. 



TABLE I 

COMPONENTS OF RADIATION IN THE UPPER ATMOSPHERE AND 
THE EXPERIMENTAL RADIATIONS SELECTED 



UPPER ATMOSPHERE 


EXPER 


1MENTAL 


RADIATION 


ENERGY OR 
WAVELENGTH 


AMOUNT 


ENERGY OR 
WAVELENGTH 


AMOUNT AND 
SOURCE 
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-4 2 

6 x 10 mw/cm 
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10 mw/cm 


50 kev 10.24 A) 
TO 250 kev 
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(AURORAL REG- 
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10- 500 Kev 


1% 

80 7> 
19 7o 
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FLUX 
'l PARTICLE/ 
cm 2 / sec 


2 mev 


UP TO 100 /ia /cm* 
VAN DE GRAAFF 


PROTONS 


2 1 B ev, 
AVERAGE IN 


10 - 20 mev 


UP TO 100 ,ia /cm 2 
CYCLOTRON 


ALPHAS 


NUCLEAR BURSTS 
50M ev 


40 mev 


UP TO 100 fia/cn? 
CYCLOTRON 



Types of Radiation 

Selection of the radiation used in this study was based upon the 
best available information on the nature of radiation components in 
the upper atmosphere. 1 It was, however, necessary to make com- 
promises because it was not possible to duplicate these radiation con- 
ditions exactly. Table I summarizes the information given in Reference 
(1) and also lists the actual radiations available for the study. 

As result of data gathered by the I.G.Y. Satellites 1958 Alpha 
and Gamma, it may be necessary to revise Table I. 2 The number of 
counts reported by the Geiger counters on the satellites was larger by 
a factor of about 3 X 10 l than the values found by extrapolation from 



1 "The Environment of an Earth Satellite," SCEL Technical Memo- 
randum, M-1747, March, 1956. 

2 State University of Iowa Report, SUI-58-5, May, 1958, to be published. 
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the high-altitude-balloon data on which Table I is based. The nature 
of the radiation is unknown, but it appears that it may be either high- 
energy electrons or protons of moderate energy, i.e., much lower in 
energy than the 10 bev protons of cosmic rays. The implication of 
these new data will be discussed. 

In general, as broad a spectrum of electromagnetic radiation as 
possible was used so that predictions could be made concerning those 
areas where experimental sources were not available. For example, 
it is known that solar x-rays and ultraviolet radiation in the wavelength 
range between 3 A and 200 A exist in the upper atmosphere. Labora- 
tory sources in that portion of the spectrum are not available, but 
radiation sources of wavelengths shorter than 0.24 A and longer than 
2200 A are available so that the region in question can be straddled. 

As for the particle radiations, cosmic ray energies are in excess of 
10 bev. However, the highest energies available for this study were 
40 mev for alpha particles and 17 mev for protons. The results of 
irradiations at these lower energies were extrapolated to predict the 
effect of higher energy particles. Sources of radiation used during 
this study are as follows. 

Ultraviolet Radiation 

Three different commercial ultraviolet lamps were used, namely the 
Hanovia Type S-H burner, GE Type BH-6 burner and a GE Type S-4 
burner with glass envelope removed. The intensities of radiation in 
the spectral region 2200 A to 3400 A available from these sources, 
together with the amount of radiation from the sun in the same spectral 
region, are given in Table II. The high intensity available from the 
BH-6 permits the equivalent of one year's irradiation to be delivered 
to a cell in 20 hours. 

X-Rays 

The source of x-rays in the range between 50 and 250 kilovolts was 
the Westinghouse Quadrocondex X-ray machine of the Evans Signal 
Laboratories, while a 2 mev Van de Graaff Electron Accelerator (High 
Voltage Engineering Corp.) of those same laboratories was used for 
the 750 kv and 2,000 kv x-ray spectra. No niters were used to modify 
the continuous spectrum emitted by the target. Periods of irradiation 
were limited to a few hours at each machine energy. 

High-Energy Electrons 

The 2 mev Van de Graaff machine referred to above supplied elec- 
trons at an energy of 1.7 mev. The electrons were admitted into the 
air of the room after passing through a 0.005-inch aluminum window. 
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The electron current to the samples was determined by measuring the 
current to a metal plate of the same area as the samples, but thick 
enough to stop the electrons. 

Protons and Alpha Particles 

Cyclotrons at Princeton University and at the Brookhaven Na- 
tional Laboratories were used to supply 17 mev and 10 mev protons, 
respectively, while 40 mev alpha particles were supplied by the Brook- 
haven Cyclotron. These particles entered the air of the room by pass- 
ing through a 0.001-inch Dural" window. The particle fluxes used were 



Table //—Comparison between Different Criteria for Estimating the Effect 
of X-Radiation on Surfaces 



Type of 
Radiation 

Atmospheric 
X-rays 

200 ev-1200 ev 

665-1500 
1000-1500 
1200-1700 

Machine Pro- 
duced X-rays 
50 kv 
(550 r) 250 kv 
(12,000 r) 750 kv 
(8000 r) 2000 kv 
(5X10* r) 



Incident 
Power/ cm 2 
(mw) 



0.1 
0.06 

3 X 10-" 
10—1-10 — 



4.6 

532.0 
384.0 
3850. 



Power 
Absorbed/ cm 2 
(mw) 



0.1 

0.0G 

3 X 10 — 

10— »-io- 



1.95 X 10- 3 
0.20 
8.03 
53.6 



No. of 

Photons/cm 2 
sec. (A r „k) 

(Approximate) 
3 X 10» 
6 X 10 T 
1 X 10* 
5 X 10-1-10- 



5.8 X 10» 
1.3 X 10 10 
3.2 X 10° 
4.2 X 10" 



in the range of 10 s ' particles per square centimeter per second. Flux 
measurements were made by observing the beam back-scattered from 
a foil according to the following scheme: The machine current was 
increased to values a few orders of magnitude above those to be used 
for the actual bombardment. This higher current, measured by means 
of a Faraday cup and a conventional microammeter, was compared to 
the number of counts registered by a photomultiplier arranged to 
receive some of the scattered particles. It was established that the 
photomultiplier counts were proportional to the current in the Faraday 
cup. Thereafter the number of counts registered by the photomultiplier 
were used to determine the actual flux of particles. To establish how 
the particles are distributed over the irradiated area, copper foils were 
substituted for the cells and activated by the particles. The foils were 
then cut into small pieces and the activity in the regions occupied by 



* Registered trade mark. 
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the solar cells was determined as a fraction of the total activity con- 
tained in the copper. 

The Solar Cells 

The silicon solar cells were manufactured by the Hoffman Semi- 
conductor Division, Evanston, Illinois and have been described in the 
literature. 3 The particular units used were 0.5 X 2 centimeters and 
about 0.05 centimeter thick. Figure 1 shows a sketch of the unit. The 
surface of the originally n-type silicon wafer had been converted to 
p-type by diffusion doping with boron. The diffused layer, which is 
about 1 micron thick, is removed from all surfaces except the top by 
etching. The active volume of the cell from which generated carriers 
can be collected is determined by the thickness of the diffused layer 




Fig. 1 — Diagram of a silicon solar cell. 



and the exposed area. Ohmic contacts are made to a thin top strip on 
the surface and the entire bottom surface by nickel plating and 
soldering. 

At least six cells were used for each radiation type, three in air 
and three in a vacuum. The cells used had the highest values of j? mnx 
from among about 100 cells whose current-voltage characteristics had 
been measured before irradiation. For these cells rj max varied between 
6 and 10 per cent. (rj max is defined as the ratio of maximum electrical 
power output, P max , to solar power input as measured by an Epply 
pyrheliometer.) 

Experimental Arrangement and Procedure 

The samples were mounted on a water-cooled copper block, so 
oriented that the light beam (from a 150-watt projector spotlight) 
was normal to the beam of radiation, and both radiations were incident 
on the cells at an angle of 45° as shown in Figure 2. The intensity of 

3 M. B. Prince, "Silicon Solar Energv Converters," Jour. Appl. Phys., 
Vol. 2G, p. 534, May, 1955. 
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illumination was adjusted to deliver about 100 milliwatts per square 
centimeter of equivalent solar radiation to the samples. Whenever 
possible, constancy of the light source was assured by observing a 
signal from another silicon solar cell which was located well outside 
the radiation field. For the vacuum runs, a similar water-cooled block 
was suspended inside a Pyrex* vacuum chamber fitted with windows 
to admit both light and particles. The pressure in the vessel was main- 
tained at about 10 ~ 5 mm Hg by means of a conventional oil-diffusion 
pump and liquid-air traps. 



Fig. 2 — Experimental arrangement and circuit for measurement. 

Figure 2 also shows how electrical measurements were made on the 
cells. A current-voltage characteristic like that shown in Figure 3 
was obtained by shunting in succession each of fifteen different pre- 
cision resistors across the illuminated cell. The voltage developed 
across the shunt resistor was recorded on a Brown Electronik Recorder. 
From the resistance and voltage, V max , i s , P m „ x and 7? lnax were computed. 
Each run took three minutes and was repeated periodically during 
irradiation. During some runs only three precision resistors were 
monitored to give V, nax , i, and P max , thus allowing for more rapid 
observations during faster decays. 



Ultraviolet Light and X-Rays 

Figures 4 and 5 show i s , V max and P max as a function of irradiation 
time with the S-H Hanovia Lamp for three solar cells in air. It is 



* Registered trade mark. 
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Fig. 3 — Comparison of current-voltage (i-V) characteristics in air and in 
a lO- 5 mm of Hg vacuum at a flux of 110 mw/cm 2 . 



evident from these figures that none of the three parameters changed 
during irradiation for a time nearly equivalent to one year in the 
upper atmosphere. Similar results were obtained using the S-4 burner 
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on cells in vacuum; no changes were observed for periods equivalent 
to about 2.5 years in the satellite environment. The samples were 
irradiated by the BH-6 burner for a period equivalent to 50 years in 
the satellite environment for the spectral range where the burner over- 
laps the upper atmosphere ultraviolet spectrum. Figure 6 shows P max 
versus time for this run, and indicates that a decrease equivalent to 
no more than 0.3 per cent per year may have occurred. The "decay" 
may, however, have been an experimental artifact traceable to the 
difficulty of maintaining constant illumination over the 1000-hour 
laboratory irradiation period. With this reservation in mind, assuming 
that the decay is real, one finds by extrapolation that 85 years would 
elapse before P max would decay 25 per cent. 

Results for irradiation at 50 and 250 kilovolts with x-rays are shown 
in Figures 7 and 8. These runs were made on successive days, some- 

EFFECTIVE YEARS — — 




O 100 200 300 400 500 600 700 600 900 1000 

HOURS—*- 

Fig. 6 — Pmas versus time for high intensity ultraviolet source. 
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times in vacuum and sometimes in air. As is evident from the figures, 
there was no monotonia change in any of the parameters; the scatter 
of the points in these figures is an indication of the limitations of the 
experimental method rather than of any changes occurring in the 
samples. Similar results were obtained with 750 and 2,000 kv x-rays. 
Thus, it is concluded that these x-rays do not affect the operation of 
solar cells. Table II shows both the total dose in roentgens (where 
applicable) and an estimated total number of incident photons, N ph , 
received by the samples during the x-ray and ultraviolet irradiations. 
The quantity, N ph , for the x-rays was computed by assuming that the 
average energy of a photon in the x-ray spectrum was % of the 
maximum photon energy. 

There remains, however, the problem of inferring the life of these 
cells in the x-ray environment described in Table I in which the x-ray 
photon energy is much lower than that used for the irradiations de- 
scribed. Since even the 2,000 kv x-rays did not produce any discernible 
bulk damage it can be concluded that if the low-energy x-rays actually 
present in the upper atmosphere are going to produce a decay of the 
current-voltage characteristics of these cells, it will necessarily be 
the result of a surface effect. Such a surface effect could be propor- 
tional to one of the following parameters of the radiation: (a) the 
integrated number of photons, (b) the incident energy per square 
centimeter, (c) the energy absorbed near the surface per square centi- 
meter. This latter quantity is determined by calculating the energy 
absorbed in the first micron of cell thickness. 

Table II compares the values of these parameters for the x-radia- 
tions likely to be encountered above the atmosphere and the x-rays used 
in these experiments. All values are listed per unit time so that from 
the laboratory irradiation time one can determine directly the equiva- 
lent irradiation time in the satellite environment. As we have indicated 
above, P mnx was constant during the laboratory x-irradiations within 
an experimental error of ± 5 per cent. To arrive at a lower limit on 
the life of the cells in the satellite environment, assume that they had 
experienced a decay of 5 per cent, which was not detected because of 
the limited precision of the experiments. Table III shows how much 
time will be required for the cells to lose 25 per cent of their initial 
P, mx in the satellite environment based on these assumptions. 

High-Energy Electrons 

Three units were run individually at atmospheric pressure and 
three were run simultaneously in vacuum. The distribution of elec- 
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Table III — An Estimate of the Lower Limit on- the Life of Cells Above the 
Atmosphere Based on Three Different Criteria (X-rays) 

Type of Radiation Incident Power Power Absorbed N„ h 

50 kv and 250 kv 2.4 X 10 4 hrs. 9.0 hrs. 39 hrs 

750 kv and 2000 kv 3.8 X 10 5 hrs. 4.9 X 10 3 hrs. 70 hrs 

50 kv and 250 kv 8.1 X 10 1 hrs. 27 hrs. 103 hrs. 

trons on the three samples in vacuum was not known accurately so 
that some of the difference in decay rates of these three units should 
be attributed to this nonuniform electron distribution. The units in 
air were bombarded individually hence no distribution problem existed. 
Figures 9 and 10 show V nwx , i s , and P max plotted against the estimated 
total number of electrons incident on the sample. Table IV summarizes 
the results of the electron irradiation; the symbol <j> 0 represents the 
average flux at which the P max falls to 75 per cent of its initial value. 
It is evident from the table that the value of <j> c is different for the 
samples irradiated in vacuum and those irradiated in air. The reason 
for this difference is uncertain. However, this does not affect estimates 
of the life of cells in the upper atmosphere, since the best available 
data on the electron flux indicates a value less than 10~ 2 per square 
centimeter per second. On this basis, P mnK will decrease 25 per cent 
from its initial value after an irradiation time whose value lies be- 
tween 10 15 and 10 17 seconds, i.e., between about 10 s and 10 10 years. 

Protons 

Two groups of three units were irradiated, the first with 10 mev 
protons in vacuum at Brookhaven, and the second with 17.6 mev protons 
in air at Princeton University. The results are shown in Figures 11 
and 12 where i t , F mnx and P max are plotted against the integrated 
flux, <f>. Differences in the value of <j> 0 are partially due to the uncer- 
tainty about distribution of the flux. As with the electrons, these 
variations in the value of <j> c do not seriously affect an estimate of the 
elapsed time for P tmt to be reduced by a predetermined amount in the 
upper atmosphere. Using the data from the run in vacuum at Prince- 
ton University, we find that on the average, about 3.5 X 10 10 particles 



Table IV — Flux of Particles Required to Reduce Solar Cell Power Output 
to 75% of Its Initial Value. 

Particle Type 0 C (Particles/cm 2 ) 



1.7 mev electrons Between 5 X 10 13 (vacuum) and 

1.5 X 1015 ( ah .) 
17.6 mev protons 3.5 x 10 10 

40 mev a particles 4.4 x 10 n 
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per square centimeter are needed to produce a 25 per cent decrease in 
P mnx , as recorded in Table IV. 

Because of the difference in energy between the protons available 




INTEGRATED FLUX (PflRTICLES/cm ) * 10 — »■ 

Fig. 12 — P nl!lx versus integrated flux of 17.6 mev protons in air. 
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from the cyclotrons and the high-energy cosmic ray protons, the rela- 
tive numbers of defects produced by particles of different energies 
must be calculated. Following Seitz and Koehler, 4 note that the density 
of defects, N D , produced by an integrated flux, 0, of high-energy par- 
ticles penetrating a solid whose thickness is less than the range, i.e., 
a thin target, is given by 

N D = m <j>, (1) 



where n is the number of primary displacements per particle and v is 
the average number of secondary displacements per primary displace- 
ment. If two particles produce equal numbers of defects, it follows that 

"l n l 4>1 = v 2 n 2 $2- 

For a given species of particle of energy E, 

n = k/E so that 

K K 

vi — <j>i = vo — (3) 
E 1 ' E, 

and thus 

VI En <j>, 

— — = — . (4) 

E\ Vo <f>i 

Substitution into this equation shows that approximately 330 ten-bev 
protons must pass through a cell to produce as many defects as a single 
twenty-mev proton does. Since it is estimated that the proton flux in 
the satellite environment is about 1 particle per square centimeter per 
second, it follows that about 3.6 X 10 5 years are required to decrease 
P max by 25 per cent as a result of bombardment by cosmic ray protons 
of the upper atmosphere. 

Alpha Particles 

Six cells were irradiated, three in air, three in vacuum. As before, 
some of the variation in values of ^> ( . from sample to sample is caused 
by uncertainties in determinations of the flux received by the sample. 
The average value of </>,. is shown in Table IV. To estimate the life of 



4 F. Seitz and J. S. Kcehler, "Displacement of Atoms During Irradia- 
tion," Solid State Physics, Vol. 2, Academic Press, New York, N. Y., 1956, 
pp. .107-448. 
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such units in the alpha particle flux of the upper atmosphere, we com- 
bine the flux rate estimated in Table I, i.e., 0.2 alpha particles per 
square centimeter per second, with the average value of <£ c deduced 
from the experiments and take into account the difference in energy 
between the cosmic-ray alpha particles and those available during this 
study. When this is done, one finds that these cells should last 1.4 X 10 5 
years before 25 per cent decay in P majs is reached as a result of the 
alpha-particle bombardment in the upper atmosphere. 

Modified Estimates of Life Based On I.G.Y. Data 

If the level of radiation in the satellite environment is as high as 
given by Van Allen based on data from the first satellites, 2 our estimates 
of the useful life of silicon solar cells must be revised downward. As- 
sume, first of all, that the increased radiation is given by the lower 
limit deduced from the satellite, i.e., greater than previously estimated 
by a factor of 3 X 10 4 . Then if the radiation consists of protons 
of about 20 mev, P 1]mx from the cells would be reduced by 25 per cent 
after about 10 years in a field of this intensity. If the particles are 
electrons, then the high count rate is very likely to have been caused 
by bremstrahlung produced by these electrons bombarding the satellite 
shell and Geiger counter case. As has been shown, x-rays have essen- 
tially no effect on the cells. If we assume that electron current is given 
by the count rate, then such an arrival rate would cause a 25 per cent 
decrease in P maI after a minimum of 3 X 10 3 years, if nothing were 
interposed between the cells and the radiation field. 

Correlation with Theory 

Nuclear radiation can change three parameters which affect semi- 
conductor devices — the conductivity, <r, the minority carrier lifetime, 
t, and the surface recombination velocity, s. Changes in a and t have 
been attributed to a common cause, namely the production of excess 
vacancy-interstitial pairs associated with the displacement of atoms in 
the lattice by the bombarding particles. Such displacement can occur 
only if the energy transferred to a lattice atom exceeds a threshold 
energy, E,, which has been found experimentally to be 13 ev for 
silicon. 6 To transfer such an amount of energy in a collision with a 
silicon atom, an electron must have an energy of 145 kev, while the 
heavier proton needs 98 ev, and the « particle 30 ev. By means of 
Compton scattering, x-ray photons can impart all of their energy to 

5 J. J. Loferski and P. Rappapoit, "Radiation Damage in Go and Si 
Detected by Carrier Lifetime Changes: Damage Thresholds," Phys. Rev., 
Vol. Ill, p. 432, July, 1958. 
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a single electron during an encounter, so that photons whose energy 
exceeds 175 kev can potentially displace silicon atoms. Earlier studies 
of radiation damage in electron-voltaic cells showed that changes in r 
destroy the usefulness of a junction long before changes in a become 
important. It was therefore expected the same would be true for these 
solar cells. 

To provide a framework for discussion, the manner in which 
changes in r would be reflected in the parameters of the photovoltaic 
effect will be reviewed. Figure 13 shows a cross-section of a p-n junc- 
tion where the p-type region to the left of the junction, has a thick- 
ness x. Assume that there is no contribution to the photovoltaic effect 
from minority carrier generation in the n-type region to the right of 
the junction in the figure. This assumption is valid in this case, since 
the light which produces the effect is absorbed principally between the 
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Fig. 13 — Electron energy versus distance from front of cell. 



surface of the material and the junction. A mathematical treatment 
of this model has been presented elsewhere. 0 Here it is sufficient to 
point out that the fraction of injected carriers which survive long 
enough to pass over the junction is proportional to t 1/2 . Now t is 
inversely proportional to the number of recombination centers, N r , in 
the material. Such recombination centers are introduced into the for- 
bidden energy gap by vacancy-interstitial pairs. For low irradiation 
fluxes it is reasonable to assume that the number of defects is propor- 
tional to the integrated flux, <j>, so that 



oc A r r oc = N B t, 



(5) 



where N B is the number of bombarding particles per square centimeter 



6 P. Rappaport, J. J. Loferski, and G. E. Linder, "The Electron-Voltaic 
Effect in Germanium and Silicon P-N Junctions," RCA Review, Vol. XVII, 
p. 100, March, 1956. 
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per second and t is the irradiation time. Thus, plots of t -1 against t 
are linear for a constant flux. This relation has been verified by earlier 
work, 5 according to which the significant parameters of the photovoltaic 
effect depend on t and therefore on <f> in the following way : 



i s oc T i/2 oc 0-1/2 

L OC T-l/2 OC <f,l/2 



(6) 
(7) 



<^max = m oc In r cc In <f>~ 



(8) 



where A = e/kT and i„ is the reverse current of the junction 
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Pmax ~ h V mas °c 0-1/2 l n 0-1. 



(9) 



These relations will be true provided that only t changes during irradi- 
ation. Change of <r will not affect i s and will not produce much of a 
change in the other parameters. 

The validity of the relations shown above for high-energy (> 145 
kev) electron bombardment of silicon alloy diodes has been partially 
confirmed in previous experiments. The data from the present experi- 
ments deviates considerably from the predictions of Equations (6) 
through (9). For example, i s ~ 2 was found to be proportional to <j> n , 
where 1 > n > y 2 instead of the expected value of n = 1. This is 
shown in Figure 14 where ln i s plotted versus ln 0 for four units 
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irradiated by 2 mev electrons. Further experiments are needed to 
explain these deviations. 

Changes in surface recombination velocity could also affect V m .„ 
and />„,„, although it is not expected that surface changes would follow 
the relatively simple form of Equations (G) through (9). Radiation 
which is composed of quanta not sufficiently energetic to produce bulk 
displacements could affect the efficiency of a photovoltaic call by chang- 
ing the surface recombination velocity. Such changes might be caused 
by removing silicon or foreign atoms from the surface, or perhaps, by 
ionizing gas atoms in the vicinity of the surface and the subsequent 
adsorption of such ions onto the surface. If the latter were true, one 
would expect a different behavior of the cells during irradiation in 
vacuo and in air which has not been observed unequivocally. 

Conclusions 

On the basis of these experiments, it is concluded that silicon solar 
cells will not experience any serious deterioration for periods of the 
order of 10 5 years as a result of the radiations to which they will be 
exposed in the upper atmosphere assuming that the actual amount of 
radiation encountered corresponds to the best current estimates. Spe- 
cifically, the only radiations which produced readily observable damage 
were 20 mev protons, 40 mev alpha particles and 2 mev electrons. 
This damage has been identified as a change in the bulk properties of 
silicon. From our data, it is possible to extrapolate to the time which 
shall have elapsed before the cells lose 25 per cent of their initial 
power output in the alpha particle and proton fluxes of the upper 
atmosphere. Such extrapolations indicate a minimum life of 3.6 X 10 5 
years with protons, 1.4 X 10 3 years with alpha particles and 10 s years 
with electrons. If minimum life is calculated on the basis of recent 
satellite data, these figures would be reduced by four orders of mag- 
nitude, or even more, provided that the radiation involved is corpus- 
cular. Such a reduction in life is in the realm of speculation until 
further experiments with satellites yield more complete data on the 
nature of upper atmosphere radiation. 

The electromagnetic radiations used did not produce any changes 
which could be labeled unequivocally as a decay in the power output. 
For ultraviolet radiation, it could be argued that P max may decrease 
by 25 per cent after 85 years. As for x-rays, the best estimates of the 
type of electromagnetic radiation above the atmosphere indicate no 
x-rays of energy higher than 1,200 ev so that none of the electromag- 
netic radiation to be encountered could produce such bulk damage. The 
most pessimistic possible interpretation of the data, according to 
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which a 5 per cent decay during the irradiation was masked by experi- 
mental error, would lead one to estimate that a 25 per cent decrease in 
P mil jc could occur after 6.7 years in the x-irradiation of the upper atmo- 
sphere. It would have been desirable to have extended the x-ray radia- 
tion times so that these extrapolations could have approached the 
minimum life based on the particle irradiation, but this was impossible 
within the time limitations of this study. It should also be pointed out 
that none of these estimates of minimum life have included the effects 
of annealing which are likely to occur and which would increase the 
minimum life. Operation at a somewhat elevated temperature should 
also appreciably increase the minimum life. No significant differences 
between irradiation in vacuum and in air were noted. 

It is important to point out that the decay of P max is approximately 
inversely proportional to the flux so that a second decrease of 25 per 
cent would require 4/3 as much flux, i.e., time, as this first decrease, 
with correspondingly longer periods for further decay. 

Acknowledgments 

Ronald Hand of these Laboratories gave invaluable assistance in 
many phases of this work. B. Markow, L. McSherry, and II. Rast of the 
Evans Signal Laboratories gave generous assistance in operating their 
Van de Graaff accelerator and in supplying information concerning 
the total dosage received by the units. Thanks are also due to Professor 
R. Sherr and Dr. C. Rrockman of Princeton University, and Dr. C. P. 
Baker of Brookhaven National Laboratories whose advice and help 
were indispensable for the Cyclotron experiments. 



THE INFLUENCE OF DEFECT LEVELS ON 
PHOTOEMISSION* 1 

By 

W. E. Spicer 

KCA Laboratories, 
Princeton, N, J. 

Summary — The influence of defect levels on the characteristics of 
photoemitters is considered. These levels may affect a photoemitter through: 

(1) The creation of photoemissive centers, 

(2) The production of a bending of the bands at the semiconductor- 
vacuum interface and the semiconductor-backing interface, and 

(3) The determination of the magnitude of thermionic emission from 
and the conductivity of the photoemitter. 

It is shown that, if photoejnission from the valence band is to be used, 
the photoemitter performance ivill be optimized by acceptor levels near 
the top of the valence band. If photoemission from defect levels is desired, 
these should lie near the bottom of the conduction band. However, quantum 
efficiency in excess of a few percent cannot be expected from defect levels. 

Introduction 

a LL PHOTOEMITTERS of practical interest such as Cs 3 Sb, 
J=\ [Cs] (NaK) 3 Sb, and Ag-O-Cs are semiconductors. In most 
-4- cases 2 the photoemission of interest is due to the excitation 

of electrons from the valence band (intrinsic photoemission), whereas 
the position of the Fermi level and thus the conductivity and thermionic 
emission of the photoemitter are determined by defects. 36 It is the 



* Manuscript received August 26, 1958. 

1 The term "defect" level covers any energy level introduced into a 
semiconductor by an imperfection in the lattice whether by a foreign atom 
(impurity level), a vacancy, or any other type of defect. 

2 The infrared response of Ag-O-Cs is probably not due to valence 
band excitation. 

3 A. H. Sommer, "N-Type and P-Type Conduction in Alkali Antimonide 
Photoemitters," Jour. Appl. Phys., Vol. 29, p. 1568, November, 1958. 

4 J. Sakata, "Some Experimental Studies of the Conductivity and 
Thermcelectromotive Foi-ce of Cs ; .Sb Photo-Cathodes," Jour. Phys. Soc. 
Japan, Vol. 8, p. 125, January-February, 1953; "Hall Effect in Cesium 
Antimonide With Audio-Frequency Currents," Jour. Phys. Soc. Japan, 
Vol. 9, p. 1030, November-December, 1954. 

s J. E. Davey, "Thermionic and Semiconducting Properties of (Ag)- 
Cs 2 0, Ag, Cs," Jour. Appl. Phys., Vol. 28, p. 1031, September, 1957. 

0 W. E. Spicer, "Photoemissive, Photoconductivity, and Optical Absorp- 
tion Studies of the Alkali-Antimony Compounds," Phys. Rev., Vol. 112, 
p. 114, October, 1958. 
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of the order of 2 X 10 4 per centimeter are obtained. Thus, only about 
10 per cent of the light is likely to be absorbed near enough to the 
S-V interface to produce photoelectrons. However, not all of the 
absorbed photons will produce photoelectrons. The maximum quantum 
efficiency obtained from photoemitters where all of the light was 
absorbed within 500 A of the surface (intrinsic photoemitters) is 
about 30 per cent, 2 and there is no reason to assume that defect centers 
are any more effective in producing photoelectrons. Therefore, at best, 
about 30 per cent of the electrons excited from defect levels can be 
expected to escape from the photoemitter. Thus, under optimum con- 
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Fig. 2 — Photoemission from K 3 Sb. Defect photoemission at about 2.1 and 
2.8 ev is produced by potassium addition. The response of S-708 is not 
followed to higher hi' since it was a thick layer in a tube of such a geometry 
that only illumination from the glass substrate was possible. 

ditions, the maximum efficiencies to be expected from defect centers 
would seem to be of the order of a few per cent. It is possible that 
the infrared peak of the Ag-O-Cs photocathode whose yield is of this 
order is due to defect levels. 5 However, this has not been definitely 
established. 

Band Bending 



Because of the presence of surface states at the S-V interface and 
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of a second solid at the S-B interface, it is to be expected that internal 
electric fields will be set up in the photoemitter. This will result in a 
bending of the conduction and valence bands with respect to the Fermi 
level. Such bending may help or hinder photoemission. 

First, let us examine the S-V interface (see Figure 3). Here either 
donor or acceptor surface states may be present. With n-type surface 
states which originally lay above the Fermi level, the surface states 
must be emptied of their electrons until all of the states which remain 
filled lie below the Fermi level. This will result in a positive charge 
layer at the surface. Thus the bands will bend downward at the sur- 
face. Conversely, it can be shown that acceptor states which lay below 
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Fig. 3 — The type of band bending which may be expected for surface states 
near the center of the forbidden gap. The Fermi level is indicated by the 

dashed line. 



the Fermi level before the surface charge is redistributed will cause 
the bands to bend upward at the surface. These possibilities are illus- 
trated in Figure 3 as are the cases where no bending will take place. 
Here we have assumed surface states near the center of the band and 
recognized that the Fermi level will lie in the upper half of the for- 
bidden gap for n-type material and in the lower half for p-type material. 

From these considerations one can make a simple generalization 
about the bending of the bands. With n-type bulk material, the bands 
will tend to bend upward, whereas with p-type material, the bands will 
tend to bend downward. 

From Figure 3, the effect of the band bending on photoemission 
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may be deduced. 12 If the bands bend upward, an electron excited in 
the interior of the material must overcome a retarding potential in 
moving to the surface. Whereas, with bands bending downward, the 
electron will be effectively accelerated toward the surface. Thus bands 
which bend upward will hinder photoemission and those which bend 
downward will aid it. 

It has been found 3 - 0 in the family of photoemitters, [Cs] (NaK) 3 Sb, 
Cs 3 Sb, Rb 3 Sb, K :! Sb, and Na 3 Sb, that the best emitters are p-type in 
bulk conductivity. How much this has to do with band bending and 
how much with the inherent properties of these materials has not been 
established. It is possible that the effect of oxygen on the Cs 3 Sb 
cathode, i.e., the shift of the threshold toward longer wavelengths, is 
due to the introduction of new surface states which produce a favor- 
able band bending. 
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SEMI-CONDUCTOR 




SEMI-CONDUCTOR 

(a) (b) 

Fig. 4 — Band bending at the interface between a metal and the photoemitter. 
is the work junction of the metal and <Psc is that of the semiconductor. 

Thus, there is, at present, little direct evidence of the effect of band 
bending on photoemission, probably because the effect is difficult to 
detect and because a concerted attempt has not been made to detect 
it. However, band bending has been so well established for semi- 
conductors in general 13 that it is certain to take place also in the 
photoemitters. 

In discussing the situation at the interface between the semicon- 
ductor and its supporting contact, we will assume, for simplicity, that 
the supporting material is a metal (similar arguments may be made 
for semiconducting substrates). If the semiconductor work function 
is greater than that of the metal, the situation shown in Figure 4(a) 
is obtained (an ohmic contact for electrons). Whereas if the semi- 



12 For the band bending to be of practical importance, there must be 
an appreciable potential drop within 500 A of the surface. Due to the high 
impurity concentrations present in photoemitters 46 , this condition is likelv 
to be fulfilled. 

13 For example, see the review articles: R. H. Kingston, "Review of 
Germanium Surface Phenomena," Jour. Appl. Phys., Vol. 27, p. 101, Feb- 
ruary, 1950. 
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conductor work function is smaller than that of the metal, the situation 
shown in Figure 4(b) is obtained (a rectifying contact for electrons). 
Remembering that the thickness of photoemitters may be comparable 
with the depth from which photoemission can occur, it is evident that 
such bending of the bands at the supporting surface-photoemitter 
interface can affect the photoemission. In particular, the case shown 
in Figure 4(b) will tend to accelerate electrons away from the support- 
ing surface toward the vacuum and thus enhance the photoemissive 
efficiency, whereas that shown in Figure 4(a) may hinder photo- 
emission. 

Thermionic Emission and Electrical Conductivity 

The thermionic emission from a photocathode at operating temper- 
ature determines the dark current in the tube and, thus, the ultimate 
noise level. For practical purposes, it is therefore desirable to reduce 
the thermionic emission as much as possible. Another important prac- 
tical consideration is that the photocurrent which can be drawn from 
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Fig. 5 — Energy level diagram for a semiconductor without band bending. 
X is the electron affinity, p. is the position of the Fermi level measured from 
the bottom of the conduction band and Eg is the band gap. 

a cathode should not be limited by high cathode resistance. (Cathode 
resistivity is usually not a problem at room temperature but often 
becomes excessive at lower temperatures.) It will be shown here that 
both low thermionic emission and low resistance may be achieved 
simultaneously. 

Thermionic emission, as well as conductivity, of a semiconductor 
is determined by the position of the Fermi level. The expression for 
the thermionic emission is 

J = A 0 T^e-<x+n)/kT ! (1) 

where / is the thermionic current density, T is the absolute tempera- 
ture, k is Boltzmann's constant, x is the electron affinity, ft is the 
position of the Fermi level, and A 0 is a constant, x + f- — $ i s the 
thermionic work function. x> /-<•> an d <fr are shown in Figure 5. /i is 
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measured downward from the bottom of the conduction band and is 
a positive number. 

The expression for the conductivity is 

o- = N c cv e-i'/w (2) 

where N c is the effective density of states in the band in which the 
conduction takes place, t is the electronic charge, v is the carrier 
mobility, and p.' is the position of the Fermi level measured from the 
nearest band. Thus, for an n-type semiconductor, // = /i ; for a p-type 
material, // is measured from the valence band and /t = E G — pf. For 
a simple system with only a single "defect" level, the value of /<. will 
be given by 

E d kT N c 

« d = + In (S) 

2 2 N D 



for a donor, and 



E a kT N c 

V-a = Eo in (4) 

2 2 N A 



for an acceptor level. Here E d is the energy of the donor level meas- 
ured from the bottom of the conduction band, N D is the density of 
donors, E G is the band gap, E A is the acceptor activation energy 
measured from the top of the valence band, and N A is the density of 
acceptor levels. For compensated materials, the Fermi level will 
roughly coincide with the dominant defect level. Thus, it is apparent, 
depending on the defects present, that the activation energy for 
thermionic emission from a given photoemitter may vary from x to 
E G + x) whereas the activation energy for conductivity may vary from 
zero to E G /2. Both minimum thermionic emission and maximum con- 
ductivity may be obtained with an acceptor level located near the top 
of the valence band. 

In experimental confirmation of this type of behavior, it is inter- 
esting to note that the multi-alkali [Cs](NaK) 3 Sb cathode has a 
threshold of response extending to lower energy than the oxygen- 
treated Cs n Sb (S-4 or S-ll surface) 11 while the dark current (ther- 
mionic emission) is roughly identical for these two cathodes. 1445 Thus 



14 A. H. Sommer, "Multi-Alkali Photo Cathodes," Trans. I.R.E. PGNS, 
p. 8, November, 1956. 

15 R. W. Engstrom and W. Widmaier, private communication. 
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it appears that for the multi-alkali cathode, the Fermi level lies closer 
to the valence band than for the Ca 3 Sb surface. This is confirmed by 
the fact that the Cs 3 Sb surface becomes insulating at 77°K, whereas 
the multi-alkali surface does not. 

Schaetti 10 has attempted to add a number of different impurities 
to Cs 3 Sb and has investigated their effect on thermionic emission. 
He found a decrease by about a factor of ten in dark current with the 
addition of Tl, Mg, Ce, and Ba. Thus it would seem that the addition 
of these elements lowers the Fermi level in Cs 3 Sb. It is not possible 
to go beyond this qualitative statement because the experimental con- 
ditions in Schaetti's work were too complex to allow more detailed 
conclusions. 

Conclusions 

It has been shown that for the minimum thermionic emission (and 
thus noise), maximum conductivity, and optimum effect of band bend- 
ing, p-type defects near the top of the valence band are needed. How- 
ever, if photoemission from defect levels (rather than the valence 
band) is desired, n-type levels near the bottom of the conduction band 
are desirable. It is estimated that the maximum quantum efficiency 
which can be expected from defect levels is about one per cent. 
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Summary — A new and more accurate expression for the beam currcnt- 
drive voltage relationship of conventional kinescope guns is presented. It 
is shown that, to a good approximation, the space-charge-free electric field 
at the cathode can be expressed as a quadratic function of radial distance 
from the cathode center. The expression is written in terms of the potential 
applied to grids 1 and 2 and four constants which depend only on gun 
geometry. The actual beam current is then related to the space-charge-free 
electric field. Drive factors [(maximum beam current)/ (cutoff voltage) 3 ' 2 ] 
are obtained for both the grid and cathode drive conditions. 

The values of gamma and drive factor obtained are in better agreement 
with experimental results than previous expressions. Considerable interest 
presently exists in guns which exhibit an increased cathode drive factor 
(and decreased grid drive factor) because of close grid 1-grid 2 spacing 
and low grid 2 voltage. The expressions presented here clearly indicate 
this behavior. 

Introduction 



V ]f HE problem of reduction in kinescope drive requirement is 
I presently of interest for several reasons. 1 " 4 A decrease in the 
required drive voltage can lead to a significant cost reduction, 
since lower power video output stages can be used, or one stage may 
even be eliminated. A small drive voltage is important in the develop- 
ment of a-c/d-c type television sets, where relatively low B+ voltages 
make it inconvenient to obtain a large video output. In addition, low- 
drive kinescopes are particularly advantageous in transistorized tele- 
vision sets, where only small signal levels are normally available. 



* Manuscript received July 24, 1958. 
t RCA Laboratories, Princeton, N. J. 

i Formerly, RCA Laboratories, Princeton, N. J.; now with Kaiser Air- 
craft and Electronics Corporation, Oakland, Cal. 

1 J. M. Shcehe, "Electron-Gun Design Permits Low Ens," Electronic 
Equipment, p. 34, April, 1957. 

2 J. W. Schwartz, "The Annular Geometry Electron Gun: A New 
Device," I.R.E. National Convention Record, Vol. 0, Part 3, p. 13, 1958. 

3 J. W. Schwartz, "A New High Transconductance Electron Gun for 
Kinescopes," RCA Review, Vol. XIX, p. 232, June, 1958. 

4 D. Sillman, "Design Considerations in Transformerless Single Recti- 
fier TV Receivers," I.R.E. National Convention Record, Vol. 6, Part 7, 
p. 133, 1958. 
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This paper discusses the drive and gamma characteristics of con- 
ventional triode-type kinescope guns. A theoretical analysis 5 is pre- 
sented which leads to an expression for beam current in terms of the 
applied grid potentials and four geometric constants. Experimental 
results are in good agreement with the derived expressions. 

Basic Equations 

The beam current in a conventional kinescope is determined pri- 
marily by the geometry and applied potentials in the triode section. 
The masking apertures in succeeding electrodes normally intercept a 
small percentage of the transmitted triode current. Although the 
percentage may increase with increasing beam current, thus reducing 
the over-all gamma, this effect will be ignored. 

The Langmuir-Child law may be used to express the current 
density, /, from a small segment of the cathode in" terms of the space- 
charge-free electric field, E, normal to the segment. Thus, 

j = kE*'"-, (1) 

where fc is assumed constant. Under the Langmuir-Child assumptions 
for plane parallel electrodes and complete space-charge limitation, fc 
contains a factor 1 \/d, where d is the electrode spacing. However, 
there are theoretical and experimental reasons for assuming fc to be 
constant. 

Kinescope control characteristics do not follow the familiar three- 
halves power law. This is because the active cathode emission area 
varies with drive voltage. Because of the axial symmetry of the elec- 
trode structure, the grid apertures, and the applied potentials (grid 1 
negative, grid 2 positive) , the field at the cathode decreases with radial 
distance, r, from the center of the cathode. At some critical radius, 
r,., depending on the applied voltages, the field becomes zero. At greater 
distances (r > r,.) the field is negative. Therefore, essentially no 
emission takes place outside a circle of radius r c . In order to find the 
cathode current, an expression must first be found for the functional 
dependence of field on radial distance from the cathode center. 

From electrostatic theory, the potential at any spatial point can be 
expressed in terms of the electrode potentials and geometric coefficients 
of proportionality. Therefore, it follows that the electric field at the 
cathode can also be expressed in the same terms. Thus, 



5 See also K. Schlesinger, "Transfer Characteristics and Mu Factor of 
Picture Tubes," Proc. I.R.E., Vol. 41, p. 528, April, 1953. 
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E (0, r) = K, (r) V 1 + K 2 (r) V 2 , (2) 



where F x and 7, are the grid 1 and grid 2 potentials, respectively. 
The coefficients K 1 (r) and 7v 2 (r) can each be expressed as a power 
series in r. Equation (2) is, therefore, rewritten as 

E(0,r) = (A+Br2 + ...)F 1 + (C - Dr- + ■ ■ •) V 2 , (3) 

where A, B, C, and 73 are arbitrary geometric constants. For practical 
purposes, the higher order terms can be neglected. This is justified 
as follows : 

In axially symmetric fields, the potential V(z,r) at any point can 
be expressed in terms of the axial potential <j> and its derivatives. 6 

Thus, 

* y+i y»6 

V(z,r)=4> <j!."H </> (1) <£< 0) + • • • (4) 

4 64 2304 

Differentiating Equation (4) with respect to z, 



r- v r 



E (0, r) = <f>' H ^ <j>™ + ■■• (5) 

4 64 2304 

The terms in Equation (5) diminish rapidly at sufficiently small values 
of r. For the purposes of this derivation, only the first two terms will 
be retained. Figure 1 indicates that errors resulting from omission of 
the higher order terms are significant only at values of r > r c . 
Thus, Equation (3) can be rewritten as 

E (0, r) = (A + Br n -) \\ + (C - Dr?) V 2 . (6) 

The value for the critical radius, r c , is obtained by setting Equation 
(6) equal to zero. This gives 



AV 1 + CV 2 

(7) 



DV 2 — BV l 
Combining Equations (6) and (7), 

E(0,r) = (r/-r*)(DV 2 -BV 1 ). (8) 



6 I. G. Maloff and D. W. Epstein, Electron Optics in Television, McGraw- 
Hill Book Co., Inc., New York, N. Y., 1938, p. 77. 
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Equation (1) can be expressed in integral form to give the total 
cathode current, /; 



= 2irk I 
J o 



E(0, r)W-r dr. 



Upon integration, this becomes 

(AVj + cy 2 ) 5 / 2 



/ = M 



{DV 2 -BV 1 ) 



2tt 

where M = k. 

5 



(9) 



(10) 




S o 2 4 G 8 10 IS 14 16 18 20 

RADIAL DISTANCE FROM CATHODE CENTER, r, MILS 

Fig. 1 — Comparison of measured and calculated electric field, from 
resistor-board data. 

In deriving this equation, it should be noted that space-charge 
effects were not rigorously taken into account. Although the Langmuir- 
Child law does include space-charge effects, Equation (2) was obtained 
from purely electrostatic considerations. It ignores the effects of the 
charge distribution on the potential near the cathode. Space-charge 
effects are accounted for in the selection of the constant M. As shown 
in a later section, the expressions derived using these approximations 
are in good agreement with experimental measurements. 

Equation (10) is used to obtain expressions for a drive factor and 
gamma for both grid drive and cathode drive conditions. These ex- 
pressions are in terms of the constants A, B, C, and D. Equation (6) 
can be used to obtain the values of these constants from data taken on 
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an electrostatic analog such as a resistor board. 1 



Drive Factor 

For both grid and cathode drive, the drive factor is defined as the 
ratio of maximum beam current / niax (current at V 1 = 0) to the three- 
halves power of cutoff voltage V co . It has the same units as perveance, 
is a function only of the tube geometry and the cathode emission prop- 
erties," and is independent of applied voltages. The drive factor can 
be expressed in terms of the constants previously defined, as follows: 
From Equation (10) at V x = 0, 

C5/2 

I™* = M ¥ ! . (11) 

D 

For grid drive, / = 0 at V 1 = V co . Thus, 

C 

V co = - -V. 2 . (grid drive) (12) 

A 

Combining Equations (11) and (12), 

C 

/ ma *=vl/ A=/2 (_y co) 3/2. (grid drive) (13) 

D 

The grid drive factor is therefore given by 

/ c 

* max w 

I'a„ - = M ■ A-v-\ (grid drive) (14) 

(-7 C0 )"/= D 

For cathode drive, the relationships 

Vi = -V k , (15a) 

V 2 = V 2 '-V k (15b) 

must hold, where V 2 ' is the grid 2-ground voltage and V k is the cathode- 
ground voltage. Substituting into Equation (10) and proceeding as 
for grid drive, 

7 G. Liebmann, "Solution of Partial Differential Equations with a Re- 
sistance Network Analogue," British Journal of Applied Physics, Vol. 1, 
p. 92, April, 1950. 

* In practice, even under complete space-charge-limited conditions, 
cathode currents are not strictly independent of the cathode properties. 
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C 

IV, (cathode drive) (16) 



A+C 



and 

/„« MC 

p = — : = (A + C) 3 / 2 . (cathode drive) (17) 

(+y co )3/ 2 ^ 

For experimental evaluation purposes, the effects of variations in M 
for different tubes can be eliminated by considering the ratio P KD /r an . 

Thus, 

P KD /A + C 



(18) 



- 1 GI> 

Expression for Gamma 



The beam current in a kinescope is generally considered to be a 
power-law function of the signal or drive voltage, V D . The exponent 
in this power law is called the tube gamma. It is usually determined 
experimentally as the slope of the logarithmic plot of current versus 
drive voltage; 

3 log / V D 3 / 

y == = . (19) 

3 log V D I 3 V D 

Equation (10) can be rewritten to express the beam current in terms 
of drive voltage, from which an expression for gamma is found. In 
doing this, grid drive operation and cathode drive operation are treated 
as separate cases. 

Grid Drive 

The grid drive voltage is given by 

V GD = V x - V co . (20) 

Substituting Equations (12) and (20) into (10), 

I = MA*r- . (21) 

BC\ 

D+ JV 2 -BV GD 
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Thus, 



7nn = 2.5 + 



D C \ Vn 
B A J Van 



(grid drive) (22) 



Cathode Drive 

The cathode drive voltage is given by 

V K o = V C0 -V K 



(23) 



Since the grid is at a-c ground potential for cathode drive, from 
Equation (15a), 



V KD = V co + V,. 
Substituting Equation (16), 

C 

V 1 = V KD - — y/. 



A + C 

Combining Equations (15a), (15b), and (25), 

A 



(24) 



(25) 



V 2 = V KD + 



V,' . 



A+C 

Inserting Equations (25) and (26) into (10), 



(26) 



/ = M (A + C) ■ 



AD + BC 
A + C 



(27) 



V/+ (D B) V Kn 



Thus, 

y KD = 2.5 + ■ 



{AD + BC) TV 
L (A+C)(B-D) V KU 



. (cathode drive) (28) 



Equations (22) and (28) are of special interest because they show 
that gamma is not constant, but depends on tube geometry and the 
applied voltages. In the absence of drive voltage, y = 2.5. As the 
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gfid-cathode voltage is increased from cutoff to zero bias, gamma 
increases to a value on the order of 3 or 4. It is also noted that the 
denominator of the second term in these equations vanishes as the 
drive voltage is still further increased. It should be realized, however, 
that the results are valid only for V 1 < 0. 

The range of values for gamma reported in the literature varies 
from 1.4 to over 3. This analysis provides a theoretical basis for the 
higher values of gamma. Measurements of gamma less than 2.5 may 
be due to an erroneously low measurement of cutoff voltage. The slope 
of the logarithmic current-drive-voltage curve (which is gamma) will 
then also be too low. Figure 2 demonstrates this qualitatively. Incor- 
rect values for drive factor also result from errors in V c „. An accurate 
method for measuring V co is described in the Appendix. 
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Fig. 2— Illustration of how errors in measuring V„ affect the apparent 

gamma. 



Geometrical Dependence of Drive Factor 

The constants A, B, C, and D have been obtained for various gun 
geometries with the aid of a resistor-board potential analog. Voltage 
measurements with two linearly independent sets of electrode voltages, 
at two values of r, are required for each geometry. The field at the 
cathode is obtained by taking the difference in potential between the 
cathode and the first measurable point in an axial direction from the 
cathode at the two different values of r. It has been found that the 
best quadratic approximation for the cathode field is obtained by 
choosing the radial positions r ■- = 0 and r - 2/3 R tl , where R gX is the 
radius of the grid 1 aperture. A schematic representation of the gun 
triode section is shown in Figure 3. 

Figure 4a shows the values of A for variable grid 1-grid 2 spacing 
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and various values of grid 2 radius, E„ 2 . Figures 4b, c, and d show 
the values of B, C, and D, respectively, for the same conditions. Figures 
5a, b, c, and d show measured values of the constants for variable 
cathode-grid 1 spacing with grid 1 thickness as parameter, for a 
grid 1-grid 2 spacing of 6 mils. Figures 6 and 7 show the values for 
the same conditions, but for grid 1-grid 2 spacing of 36 mils and 
100 mils, respectively. Using Equations (14) and (17), numerical 
values of the normalized drive factors P Qn /M and P KD /M are obtained 
from the constants A, B, C, and D for various geometrical configura- 
tions. The accuracy of the constants is in doubt for the 2-mil spacings 
because of the discrete nature of the resistor-board analog. 

Figure 8 shows curves of P Kn /M and P GD /M as functions of 
cathode-grid 1 spacing for a grid 1-grid 2 spacing of 36 mils and 
with grid 1 thickness, t u as a parameter, It is seen that the drive 
factor increases as the cathode-grid 1 spacing is made smaller. Also, 




2"gi " 2Rg2 



CATHODE GRID I GRID 2 

Fig. 3 — Schematic of gun triode section. 

the drive factor is relatively insensitive to grid 1 thickness for large 
spacings, but increases with thinner grid 1 electrodes for small spac- 
ings. Figure 9 shows that the cathode drive factor can be greatly 
increased by decreasing the grid 1-grid 2 spacing. The grid drive 
factor, however, decreases in this case. The diameter of the grid 2 
aperture has little effect on the drive characteristic except at close 
grid 1-grid 2 spacings. 

Experimental Measurements of Drive Factor 

A series of electron-gun triode sections were made with different 
spacings, so that a comparison of the analytical predictions and 
experimental results could be made. Figure 10 shows the essential 
features of a typical gun. The cap at the end of G 2 insures that 
essentially all of the beam is collected and measured at that electrode. 
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Fig. 4 — Geometric constants versus grid 1-grid 2 spacing. 
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(0) GEOMETRIC CONSTANT A vs. SPACING d k 



dk(MILS) 

(b) GEOMETRIC CONSTANT B vs SPACING d|, 
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The cathode-grid 1 spacing was varied from 2 to 10 mils (cold spac- 
ing), and the grid 1-grid 2 spacing from 2 to 140 mils. The beam 
current at zero bias and the cutoff voltage for both grid and cathode 
drive operation were measured, and the drive factor calculated. 

Figure 11 shows measured values of P KD for various grid 1-grid 2 
spacings. The solid curve, obtained from Equation (17), is theoretical. 
In this equation, a value of M = 0.75 X 10-° ampere/volt^'mil 1 / 2 was 





d k (MILS) 

(a) GEOMETRIC CONSTANT A vs. SPACING d k 



d k (MILS) 

lb) GEOMETRIC CONSTANT B vs. SPACING d k 




d k (MILS) 

(C) GEOMETRIC CONSTANT C vs. SPACING dk 



36 MILS 
IS MILS 
' 18 MILS 
i 6MILS 




d k (MILS) 

(d) GEOMETRIC CONSTANT D vs SPACING d k 



Fig. 0 — Geometric constants versus cathode-grid 1 spacing; d, = 36 mils. 

chosen because it gave the best fit for the experimental points. The 
spread in the data is attributed to variations in cathode emission con- 
ditions. 

Equation (18) indicates that the ratio P KD /P GD (which is always 
greater than unity) is independent of the emission constant, M. 
Figures 12 and 13 show the close correlation of experimental data 
with the analytical expressions for P K d/Pgd- 
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Gamma 

Figure 14 shows a logarithmic plot of beam current versus grid 
drive voltage for three different tubes with a grid 2 voltage of 200 
volts. Figure 15 shows a similar plot for one of these tubes at V 2 ' — 50 





(c) GEOMETRIC CONSTANT C vs SPACING d k (d) GEOMETRIC CONSTANT D vs SPACING d k 
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R g! = I8MILS 
Rgi ■ IBMILS 
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T 
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t_ 



Fig. 7 — Geometric constants versus cathode-grid 1 spacing; d„ = 100 mils. 
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volts and V 2 ' = 400 volts ; grid drive and cathode drive curves are 
both shown. In each case gamma is greater than 2.50, and is seen 
to increase rapidly near zero bias. Figure 1G compares the experi- 
mental curves for grid drive with the theoretical curves obtained 




from Equation (21). In this equation, a value of M = 0.62 X lO-" 
ampere/volt 3 / 2 mil 1/2 was chosen in order to permit a direct compari- 
son for both values of V 2 - 
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Fig. 9 — Drive factor versus grid 1-grid 2 spacing. 
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Fig. 10 — Electron-gun triode section. 
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Fig. 11 Theoretical and experimental values of cathode drive factor versus 

grid 1-grid 2 spacing. 
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Fig. 12 — Theoretical and experimental values of Psd/Pod versus cathode- 
grid 1 spacing. 
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LEGEND: 

o DATA FOR TUBES WITH d k » 6MILS(C0LD),t, = t 2 = 5MILS, 
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THEORETICAL CURVES FOR d k , t, AS SHOWN, t 2 = 6 MILS , 

Rg, = Rg 2 = 18 MILS 
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Fig. 13 — Theoretical and experimental values of Pkd/Pod versus grid 1- 

grid 2 spacing. 




GRID DRIVE VOLTAGE V CD ( VOLTS) 



Fig. 14 — Gamma for electron-gun triode section (three tubes, grid drive, 

Vi = 200 volts). 
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Conclusions 

Equations have been derived which are in good agreement with 
experimental observations made by the authors and others. Unlike 
previously available expressions for gamma, experimentally observed 
values greater than three are accounted for. The cathode drive factor 
is shown to be always greater than the grid drive factor and, in some 
ranges, both are sensitively dependent on triode geometry. For close 
cathode-grid 1 and grid 1-grid 2 spacings, and for thin grid 1 discs, 
the cathode drive factor increases markedly. More than a two to one 
variation can be obtained from changes in triode geometry. For grid 
drive, the effect of a thin grid 1 disc and close cathode-grid 1 spacing 
increases the drive factor, but close grid 1-grid 2 spacing reduces the 
drive factor. 
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Fig. 16 — Gamma comparison under grid drive conditions, theoretical and 
for gun triode section. 
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Fig. 17 — Measuring circuit (top-grid drive, bottom-cathode drive). 



DRIVE FACTOR ANT) GAMMA OF ELECTRON GUNS 58:5 



Appendix 

A circuit designed to determine V co to within a few tenths of a 
volt is shown in Figure 17. A fixed a-c voltage is added to a variable 
d-c bias voltage obtained from a ten-turn precision potentiometer. The 
output waveform is viewed on a sensitive oscilloscope, and the potenti- 
ometer is adjusted so that the peak of the superimposed a-c voltage 
just starts to form a small pip in the output waveform (with switch 
S open). The magnitude of the cutoff voltage then equals the d-c 
voltage minus the peak a-c voltage. With switch S closed and both 
grid and cathode grounded, the maximum or zero bias current is 
measured. For the gamma measurements, switch S is closed and the 
a-c voltage source removed. 



TWO BACKWARD-WAVE OSCILLATOR TUBES FOR 
THE 29,000-TO-74,000 MEGACYCLE 
FREQUENCY RANGE*t 

By 

Donald J. Blattner and Fred Sterzer 

RCA Electron Tube Division, 
Princeton. X. J. 

Summary— Two experimental tubes were built which demonstrated the 
feasibility of obtaining voltage-tunable oscillators having useful power 
output over very wide frequency ranges at millimeter wavelengths. The 
design and operating characteristics of these tubes are discussed. Their 
wide frequency range and milliwatt-level power outputs make them well 
suited for signal-generator and local-oscillator applications in the millimeter- 
wave spectrum. 

Extrapolation of the design indicates the possibility of tape-helix 
backward-ivave oscillators delivering useful power at frequencies as high 
as 150,000 megacycles (2-millimeter wavelength). 

Introduction 



T 



K If N HIS paper describes the characteristics and performance of two 
experimental backward-wave oscillator tubes which were built 
to explore the possibilities of obtaining voltage-tunable oscilla- 
tions over very wide frequency ranges in the millimeter-wavelength 
region. These tubes have been given developmental designations 
A-1073 and A-1074. The A-1073 oscillates between 29,000 and 49,000 
megacycles, and the A-1074 oscillates between 48,000 and 74,000 mega- 
cycles. The frequency is varied by a single voltage adjustment. This 
frequency range is considerably larger than that obtained in previously 
described backward-wave oscillator tubes 1 " 6 operating at millimeter 
wavelengths. 



* Manuscript received October 22, 1958. 

t Much of the work described in this paper was performed under a 
contract sponsored by Wright Air Development Center. 

1 R. Kompfner, "Backward-Wave Oscillator," Bell Lab. Record, Vol. 31, 
p. 281, August, 1953. 

2 w - v - Christensen and D. A. Watkins, "Helix Millimeter-Wave Tube," 
Proc. I.R.E., Vol. 43, p. 93, January, 1955. 

3 A. Karp, "Traveling-Wave Experiments at Millimeter Wavelength 
with a New, Easily Ruilt, Spaceharmonic Circuit," Proc. I.R.E. Vol 43 
p. 41, January, 1955. 

4 J. A. Noland and R. E. Lepic, "Backward-Wave Oscillators for the 
17 to 41 KMC Band," The Sylvania Technologist, Vol. X, January, 1957. 
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Table I — Electrical and Mechanical Characteristics 



A-1073 
29,000-49,000 mega- 

cycles 
050-3,200 volts 
3 milliamperes 
5-13 milliwatts 

0.031 inch 
73 

0.040 inch 

1 ampere/centimeter 2 
1,200 gauss 

RG 97/U 



A-1074 
48,000-74,000 mega- 
cycles 
750-3,100 volts 
2 milliamperes 
1.2-2.5 milliwatts 

0.021 inch 
109 

0.030 inch 

1 ampere /centimeter 2 
1,200 gauss 

RG 98/U 



Frequency 

Beam voltage 
Beam current 
Power output (calcu- 
lated) 
Helix inner diameter 
Helix turns per inch 
Cathode diameter 
Cathode current density 
Magnetic focusing field 

strength 
Output waveguide 

Both tubes use a tape-helix backward-wave circuit through the 
center of which a solid electron beam from a parallel-flow gun is 
focused by means of a uniform axial magnetic, field. The r-f output 
signals are taken out of the tubes by means of standard waveguides 
(RG 97/U for the A-1073 and RG 98/U for the A-1074). A summary 
of the electrical and mechanical characteristics of both tubes is giver, 
in Table I. 

Figure 1 is a photograph of the A-1074. Externally it is like the 
A-1073 except for the different size of the output waveguide and 
flange. 

Description of Components 

The major difficulty in designing tubes for the millimeter-wave- 
length range is that physical dimensions must be small, and therefore, 
the power-handling capacity and beam current of the tube are limited. 
Also, efficiency is low due to high r-f losses. The approach taken to 
minimize the problem of small size, with the resultant requirements 
on precision and tolerances in the mechanical design, was to make the 




Fig. 1 — A-1074 backward-wave oscillator tube. 



5 A. Karp, "Backward-Wave Oscillator Experiments at 100 to 200 
Kilomegacycles," Proc. I.R.E., Vol. 45, p. 496, April, 1957. 

8 R. W. Grow, D. A. Dunn, J. W. McLaughlin and R. P. Lagerstrom, 
"A 20 to 40 KMC Backward-Wave Oscillator," Trans. I.R.E. PGED, p. 152, 
July, 1958. 
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design of the tube as simple as possible. Figure 2 is a schematic 
diagram of the A-1073 showing the electron gun, the r-f output 
coupler, the helix, and the collector. 

Tape-Helix Circuit 

A tape helix was chosen as the slow-wave structure for the tube for 
both electrical and mechanical reasons. Electrically, the tape helix 
offers high backward-wave impedance over an octave frequency range, 
larger than that of any other known circuit, and coupling between r-f 
output and the helix over this frequency range is feasible. Moreover, 
the helix presents a comparatively large cross-sectional area over which 
the r-f field can interact with the electron beam. 

Mechanically, the helix is a good circuit to use because conventional 




Fig. 2 — Schematic diagram of tube assembly, 
grid-winding techniques can be used to hold the very tight tolerances 
required in structures which are to be used at such extremely high 
frequencies. The helices are wound of tungsten tape 0.002 inch thick 
by 0.006 inch wide for the A-1073 and 0.002 inch thick by 0.0045 inch 
wide for the A-1074. The mean diameters of the helices in the two 
tubes are 0.0315 inch and 0.021 inch, respectively. These dimensions 
were chosen to assure that the ratio of helix circumference to free- 
space wavelength is less than 0.5 over the entire range of oscillation 
frequency. The helices are made long enough to be self -terminating. 
The tape-helix circuit is supported mechanically by three round quartz 
rods spaced equally around its outside diameter. The rods are held in 
contact with the helix by a fine nickel wire spiraling around their 
entire length, as shown in Figure 2. 

Electron Gun and Collector 

The Pierce-type parallel-flow electron gun 7 shown in Figure 2 was 



7 J. R. Pierce, Theory and Design of Electron Beams, D. Van Nostrand 
Co., Princeton, N. J., 1949. 
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designed for operation with the tape-helix circuit. Confined parallel 
flow was chosen as a simple and effective method of maintaining the 
electron beam close to the circuit. Electrons are emitted from the" 
cathode in the presence of a confining axial magnetic field, and are 
accelerated to helix voltage by several electrodes whose spacing and 
voltage are so arranged as to keep the electrons moving in straight- 
line paths. The elements of the gun are: (1) an impregnated tungsten 
matrix cathode which produces a current density of 1 ampere per 
square centimeter, (2) a beam-forming electrode held at cathode 
potential, (3) an accelerating electrode, and (4) a drift tube, held 
at helix potential, which serves to collimate the electron beam before 

OUTPUT 
WAVEGUIDE 




Fig. 3 — Schematic diagram of helix-to-waveguide matching section, 
it enters the helix. The inside diameter of the drift tube is 0.032 inch 
in the A-1073 and 0.021 inch in the A-1074. 

After passing down through the tape-helix circuit, the electron 
beam is collected in the hollow cylinder of oxygen-free high-con- 
ductivity copper. This shape was chosen to reduce the number of 
secondary electrons reaching the helix. 

R-F Output Circuit 

The electromagnetic wave, which is excited in the tape helix by the 
electron beam, is radiated from an antenna at the gun end of the 
helix into the thin section of waveguide shown in Figure 3. The antenna 
connects to a post which is also the drift tube of the electron gun. 
Choke collars on both sides of the waveguide prevent radiation from 
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leaking out through the coupling apertures in the waveguide. The 
height of the section of waveguide at the antenna is 0.045 inch in 
the A-1073, and 0.030 inch in the A-1074. This dimension is then 
tapered up to standard waveguide size to that the actual output con- 
nections for the tubes are standard RG 97, U and RG 98 /U waveguides 
and flanges. In most of the experimental tubes the output waveguide 
was sealed vacuum-tight by a "gasket-type" mica window. 8 

The transfer of energy from helix to waveguide is efficient over 
the entire frequency range of the tube. The position and penetration 
of the post to which the antenna is connected as well as the antenna 
length, were selected, on the basis of experiments on an X-band scale 
model, to provide optimum matching over the desired frequency range. 
The voltage standing-wave ratio of the scale model is shown in Figure 
4 as a function of frequency. Similar data measured on an A-1073 

legend: 

■ x - sand scale model 

o experimental data on a-1073 
• experimental data on a-1074 




1 1 ' 1 1 1 1 1 1 1 1 ' 1 I I I , i . . I.ii i | 

50 55 60 65 70 75 A-1074 DATA 

FREOUENCY - MEGACYCLES x I0 3 

Fig. 4 — Voltage standing-wave ratio versus frequency for the helix-to- 
waveguide matching section. 



and an A-1074 are also shown. The voltage standing-wave ratio of the 
X-band model was below 2 throughout the frequency range of interest, 
and for the actual tubes the values were all below 3. 

Tube Characteristics 

Oscillation Frequency 

The oscillation frequency of the A-1073 and A-1074 backward-wave 
oscillator tubes is varied by changing the helix potential. In Figure 5, 
the tuning curve for the A-1073, the frequency of oscillation is plotted 
as a function of helix voltage. As the helix voltage is varied from 
650 to 3,200 volts, the frequency of the r-f output signal varies from 
29,000 to 49,000 megacycles. The data given in Figure 5 for three 
different A-1073's indicates that the tuning curve is reproduced closely 



8 F. Sterzer, "Simple High-Temperature Vacuum Tight Mica Window," 
Rev. Sri. histr.. Vol. 28, p. 208, March, 1957. 
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Fig. 5 Frequency of oscillation versus helix voltage for the A-1073. 

from tube to tube. Calculated tuning curves, plotted as dashed lines in 
Figure 5,* are based on dielectric loading factors (DLF) of 0.8 and 
0.9. A comparison between the experimental data and the calculated 
curves indicates that the DLF due to the quartz rods supporting the 
tape helix varies between the values of 0.8 and 0.9 and approaches the 
higher value at the upper end of the frequency scale. 

Figure 6 shows the tuning curve of the A-1074 when the helix 
voltage is plotted against the corresponding oscillation frequency. The 
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Fig, 6 — Frequency of oscillation versus helix voltage for the A-1074. 



* The Appendix outlines the calculations of the tuning curve, start- 
oscillation current, and power output. 
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dashed lines indicate the performance calculated on the basis of 
dielectric loading factors of 0.8 and 0.9. 

Each of the calculated curves has two parts. The lower section, 
which has a positive slope and increases with frequency, is the voltage- 
frequency characteristic of the backward space harmonic wave on the 
helix with which interaction is desired. The upper branch, which has 
a nearly constant negative slope, is the voltage-frequency charac- 
teristic of the fundamental forward wave on the helix. The two curves 
intersect at that frequency at which the helix circumference is equal 
to one half of a wavelength in the presence of the dielectric. As long 
as the helix voltage is below the intersection value, the electron beam 
can interact only with the backward wave, but when the helix voltage 
is slightly greater than the intersection value, the electron beam can 
interact with either the backward wave or the forward wave. In the 
latter case the feedback necessary for oscillation is furnished by reflec- 
tion from the collector end of the helix. 

The experimental data in Figure 6 show that as the helix voltage 
is varied from 750 to 3,100 volts, all other parameters being held con- 
stant, the frequency of the r-f output signal varies from 48,000 to 
74,000 megacycles. When the helix voltage is increased above 3,100 
volts, the forward mode of oscillation takes over and the frequency 
jumps back to 48,000 megacycles. A comparison between the observed 
and calculated backward-wave oscillation curves indicates that the ex- 
perimental DLF is 0.9; in the forward wave mode the DLF is closer 
to 0.8. 

Start-Oscillation Current 

The electron beam current required for the onset of oscillation in 
a backward-wave oscillator depends on the geometry of the beam and 
the length, loss, DLF of the r-f circuit,' and the frequency. Figure 7 
shows the start-oscillation current (measured on the collector) for 
the A-1073 plotted as a function of frequency. The dashed curve gives 
the calculated values of starting current, based on the experimentally 
measured DLF. Data taken on three different tubes are plotted. Im- 
provement in the construction of successive tubes is reflected in their 
lower starting currents, and is probably due to better helix uniformity 
and beam alignment. In each case the starting currents follow the 
shape of the calculated curve, starting below 1 milliampere and rising 
slightly at the higher frequencies. 

Figure 8 shows the start-oscillation current for the A-1074 plotted 
as a function of frequency. The dashed curve gives values of starting 
current calculated on the basis of the measured DLF, while the points 
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Fig. 7 — Start-oscillation current versus frequency of oscillation 
for the A-1073. 

show measured values. The agreement between calculation and experi- 
ment is not very good at the lower frequencies, but in every case the 
starting current lies well below 1 milliampere. 

R-F Poiver Output 

The r-f power output from a backward-wave oscillator tube is, of 
course, zero for electron beam currents below the start-oscillation 
value. As the beam current is increased beyond the start-oscillation 
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Fig. 8 — Start-oscillation current versus frequency of oscillation 
for the A-1074. 
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value, the power output rises sharply at first and then levels off for 
currents above two or three times the starting current (see Figure 10, 
which is discussed below). The output power depends upon beam cur- 
rent and voltage, beam and helix geometry, and helix loss and dielectric 
loading. The output power to be expected from the A-1073 and A-1074 
was calculated for beam currents of 3 milliamperes for the A-1073, and 
2 milliamperes for the A-1074. The curves in Figure 9 show that the 
calculated power output of the A-1073 increases from 5 milliwatts at 
30,000 megacycles to 13 milliwatts at 50,000 megacycles, and the calcu- 
lated power output of the A-1074 rises from 1.2 milliwatts at 50,000 
megacycles to 2.6 milliwatts at 74,000 megacycles. 

Measurements of the power output of these tubes were made with 
a commercial bolometer. This instrument showed a power output of 
several milliwatts for the A-1073, and power of the order of a milli- 
watt for the A-1074. Although the power was by no means monotonic 
over the band, it was observed to be continuous. However, accurate 
power measurements at millimeter wavelengths are extremely difficult, 5 
and commercially available bolometers are inefficient and are useful 
only to indicate relative power values. 

Variation of R-F Power with Beam Current 

Figure 10 shows the experimental values for the variation of r-f 
power output with beam current when all other parameters are held 
constant. Relative r-f power as measured with a bolometer is plotted 
as a function of the ratio of beam current to start-oscillation current. 
The power rises sharply to a point at which the ratio of beam current 
to starting current is between 2 and 3, and then levels off. 

Variation of R-F Power with Focusing Magnetic Field 

The output power from the A-1073 and A-1074 is sensitive to the 
intensity of the axial magnetic field applied to focus the electron beam. 
In Figure 11, the relative power output of the A-1073 is plotted against 
magnetic field strength. Relative beam transmission is shown in the 
same figure to make it clear that the loss of output power is not due 
to interception, but rather to the more subtle variation in the position 
of the electron beam. Evidently, the electrons flow closest to the helix 
over the greatest length for a magnetic field of 1,200 gauss. These 
data illustrate once again how essential it is that the electron beam 
be close to the circuit in millimeter wavelength tubes. Similar perform- 
ance was observed in the A-1074, where, again, the optimum magnetic 
field was found to be 1,200 gauss. This value is remarkably low for 



9 W. M. Sharpless, "A Calorimeter for Power Measurements at Milli- 
meter Wavelengths," Trans. I.R.E. PGMTT, p. 45, September, 1954. 
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Pig. 10 — Relative r-f power output versus the ratio of beam current to 
starting current for the A-1073 as measured on commercial bolometer. 
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tubes operating at millimeter wavelengths and is probably due to the 
simplicity of the mechanical features of the design. 

Extrapolation to Higher Frequencies 

Because the calculated and measured performance of the tubes 
were in fair agreement, the performance of tape-helix backward-wave 
oscillator tubes at higher frequencies can be predicted with some degree 
of confidence. 

An extrapolation of the present design to a backward-wave tube 
which would operate in the frequency range between 90,000 and 150,000 
megacycles results in a tube which would have a helix 0.6 inch long, 
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MAGNETIC FLUX DENSITY - GAUSS 

Fig. 11 — Relative r-f power output and beam transmission through helix 
versus magnetic focusing field strength for the A-1073 as measured on 

commercial bolometer. 

made of 0.001 by 0.0022 inch molybdenum tape wound at 218 turns 
per inch on a 0.0105 inch inside diameter. 

The simplicity of the design makes this structure mechanically 
feasible. The tube should oscillate from 90,000 to 150,000 megacycles 
as the helix voltage is varied from GOO volts to 3,300 volts. Starting 
current is calculated to be 0.5 to 1.9 milliamperes over the frequency 
range. For a current density of 5 amperes per square centimeter 
(which is not impractical with impregnated cathodes) , the collector 
current is 2.8 milliamperes and the expected power output is in the 
order of one milliwatt. It should be pointed out that while theory 
predicts certain power outputs, obtaining them over the entire operat- 
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ing range of an actual tube is another matter, and a great deal of 
developmental work must still be done before reliable performance and 
uniformity can be expected from this type of tube. 
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Appendix — Calculation op Tape-Helix Backward-Wave 
Oscillator-Tube Performance 

Calculation of helix voltage, starting current, and r-f power output 
as functions of frequency are described below. These quantities* can 
be determined reliably from the formulas if the physical properties of 
the tube, especially beam diameter and helix loss, are accurately known. 

Helix Voltage versus Frequency 

For oscillation, the velocity of the electron beam must be equal to 
(^-i)n> the phase velocity of the first backward space-harmonic wave 
on the tape helix. 

The subscript D denotes the presence of dielectric loading. There- 
fore, the beam voltage is given by 



where tj is the charge-to-mass ratio for the electron. The value of 
(V-i)i> is found from the frequency, /, helix pitch, p, the dielectric 
loading factor, DLF, and the phase velocity of the fundamental wave 
on the unloaded helix, v n , by the relation 




(1) 



2tt/ 



2tt/ 



(2) 



(DLF)v 0 
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(DLF)v 0 



(3) 



(DLF)v l 
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pf 



Therefore 



* MKS units are used. 
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The value of v 0 is quite accurately given by 

v 0 = , (5) 

2-n-a 

where c is the velocity of light in free space, p is the helix pitch, and 
a is the mean radius of the helix. A correction to Equation (5) is given 
by Pierce. 10 

The helix voltage required for synchronism between the electron 
beam and the forward wave is 

1 

forward = (DLF)v 0 * 

= 2.84 X 10- 12 (DLF) V- (6) 

Start-Oscillation Current 

The current required to start oscillation depends on the length of 
the helix, the helix impedance, and the distributed loss on the helix. 
The frequency dependence of these parameters makes the starting 
current a function of frequency. Harman 11 gives the value of CN 
required for onset of oscillation as a function of distributed helix loss. 
N is the length, I, of the helix in backward wavelengths 

I v 

N = = , (7) 

'^■-i v_ 1 

and C is the helix gain parameter 




In this expression, / is the beam current, V osu is the helix voltage 
corresponding to the oscillation frequency involved, and K_ 1 is the 



10 J. R. Pierce, Traveling Wave Tubes, D. Van Nostrand Co., Princeton, 
N. 3., 1950; Fig. 3.3, p. 25. 

11 W. A. Harman, "Backward Wave Interaction in Helix-type Tubes," 
TR-13, Electronic Research Lab., Stanford University, p. 3, April 26, 1954. 
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interaction impedance which the helix presents to the beam. Watkins 
and Ash 12 plot K_ l as a function of frequency, and show how to correct 
its value for dielectric loading effects. 



UNLOADED 



It: fa 



(DLF) 2 — 



27r/a 



(9) 



A further correction to J£_j to allow for the relationship between 
beam and helix diameters, is given by Harman 13 . Then 



START OSG 



4 [(CAT) n ' 



(10) 



Poiver Output 

The power output is computed from the efficiency expression of 
Grow and Watkins 14 



Efficiency = (4/3) RC, 

where R is an efficiency reduction factor (plotted as a function of helix 
loss by Grow and Watkins) and C is the helix gain parameter used 
above. R-F power output is therefore 

P rf = (4/3) RCIV OS0 . (11) 



This expression is valid at twice the start-oscillation current. 



12 D. A. Watkins and E. A. Ash, "The Helix as a Backward-Wave 
Circuit Structure," Jour. Appl. Phys., Vol. 25, p. 782, June, 1954. 

13 W. A. Harman, op. ext., p. 60. 

11 R. W. Grow and D. A. Watkins, "Backward-Wave Oscillator Effi- 
ciency," Proc. I.R.E., Vol. 43, p. 848, July, 1955. 



MEASUREMENT OF TRANSISTOR 
CHARACTERISTICS IN THE 3-250 MEGACYCLE 
FREQUENCY RANGE* 

By 

John IT. O'Connell and Thomas M. Scott 

RCA Laboratories, 
Princeton, N T . J. 

Summary — The transistor characteristics which- are of especial con- 
cern to the engineer applying transistors in HF-VHF circuits arc: input, 
output, and unilatcralization impedances; transconductance; and current 
gain. 

Techniques and special apparatus which have been developed for the 
measurement of these characteristics, at frequencies up to 250 megacycles, 
are described in this paper. Commonly available commercial test instru- 
ments arc used to the greatest practical extent. 

Introduction 



^ HE particular transistor characteristics which are of interest 
the design of HF-VHF circuits, such as the r-f and i-f por- 
tions of a television or FM receiver, are defined as follows: 



(1) Input impedance components: the shunt resistive component, 
R iu , and the shunt reactive component, X in , of incremental input 
impedance, for either the common-emitter or common-base connection, 
with the collector-to-common voltage held constant (collector shorted 
to common at the signal frequency) . 

(2) Output impedance components: the shunt resistive component, 
R ont , and the shunt capacitive component, C oat , of incremental collector- 
to-common impedance, with the base-to-emitter voltage held constant. 
This impedance is the same for both the common-emitter and common- 
base connections. 

(3) Unilateralization impedance components: the series resistive 
and capacitive components of the common-emitter unilateralization 
impedance. This impedance, when placed between the base terminal 
and a source of signal voltage equal and opposite to that at the col- 
lector, produces an input impedance at the base terminal which is 



* Manuscript received July 16, 1958. 
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independent, over a restricted frequency range, of collector load im- 
pedance. Unilateralization also renders the output impedance inde- 
pendent of input source impedance. The input and output impedances 
become those discussed in (1) and (2) above, except for the incidental 
shunting effect of whatever practical unilateralization circuit compo- 
nents are used. 

(4) Transconductance: the magnitude of the incremental ratio 
of collector current to base-to-emitter voltage, with either the collector- 
to-base or the collector-to-emitter voltage held constant. In the meas- 
urement procedures to be described the former connection is used 
(base common). While theoretically some difference could exist be- 
tween measurements made in the two connections, for practical pur- 
poses it may be assumed that the common-emitter transconductance 
will be identical to the common-base transconductance. 

(5) Current gain: the magnitude of the incremental ratio of 
collector current to emitter current with the collector-to-base voltage 
held constant. A knowledge of this ratio, or more particularly the 
manner in which it varies with frequency, is a prerequisite to analysis 
of the transistor as a device. 

Items (1), (2), (3), and (4) bear directly on circuit design. The 
available unilateralized gain, in terms of these quantities, is given by : 

(Transconductance) 2 R in R out 

Power Gain = . (1) 

4 

The usable gain of a tuned amplifier is usually determined by consid- 
erations relating to stability and transistor interchangeability and, 
depending on the tolerance allowed in unilateralization, may be sig- 
nificantly lower than that obtained from the above expression. 1 

A figure-of -merit 2 for the transistor itself is the frequency at which 
the available common-emitter power gain falls to unity. This is, of 
course, the maximum frequency of oscillation in a loss-free circuit, 
and is given by 




1 T. O. Stanley and D. D. Holmes, "Stability Considerations in Tran- 
sistor IP Amplifiers," Transistors I, RCA Laboratories, Princeton, N. J.. 
1956, p. 403. 

2 L. J. Giacoletto, "Study of PNP Alloy Junction Transistors from D-C 
Through Medium Frequencies," RCA Review, Vol. XV, p. 506, December, 
1954. 



coo 



RCA REVIEW 



December 1958 



The factor /„, in the foregoing expression, is termed the alpha-cutoff 
frequency, and is the frequency at which the common-base current 
gain has declined to 0.707 of its low-frequency value (3 decibels down). 
The base lead resistance, R w , which is not directly measurable, can 
be inferred from 

fa 

R b f= . (3) 

dmfg 

where g m is the intrinsic transconductance (38 millimhos per milliam- 
pere of collector current at 25 °C in an ideal transistor) and f g is the 
extrinsic transconductance-cutoff frequency, the frequency at which 
the transconductance, as measured at the transistor terminals, has 
declined to 0.707 of its low-frequency value. Alternatively, R bt > can 
be assumed to be the input resistance measured at a frequency so high 
that the intrinsic base-to-emitter capacitive reactance is low compared 
to the base-lead resistance; that is, at frequencies well above the 
intrinsic transconductance cutoff frequency. The collector-to-internal- 
base capacitance, C Vc , appearing in Equation (2) can be assumed to 
be the difference between the transistor common-base output capaci- 
tance at zero emitter current, and the incidental shunt capacitance in 
the transistor mount. Typical values for the incidental capacitance 
range from 0.5 to 1.0 micromicrofarad. 

Input and Output Impedance 

Transistor input and output impedance components can be measured 
conveniently and accurately at frequencies up to 250 megacycles using 
a Boonton type 250-A RX Meter. Both the resistive and reactive 
components of output impedance are generally within the range of 
this meter and may be measured with full meter sensitivity (from 
0.1 to 0.5 volt r-m-s across the transistor) , provided the collector bias 
voltage is greater than about three volts. 

For input impedance measurements, the voltage across the RX 
meter terminals is reduced so as not to overdrive the transistor. A 
simple modification of the RX meter, as recommended by the manu 
facturer, allows for adjustment of the signal voltage appearing at the 
meter terminals. To find the small-signal input impedance, the signal 
drive applied to the transistor is reduced by means of the added 
adjustment knob until the meter balance stops changing. 

Adaptors (Figures 1 and 2) are used with the RX meter for input 
and output impedance measurements. These adaptors are arranged 
for common-emitter test. Short leads are used in signal-carrying 
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Fig. 1 — (a) Output impedance adaptor; (b) input impedance adaptor. 



circuits. Two different values of by-pass capacitance are used in the 
adaptors, as shown in the schematics. The higher value capacitors 
are suitable at frequencies from 2 to 20 megacycles, the lower value 
capacitors from 20 to 250 megacycles. 

The range of reactive balance of the RX meter is +20 to —100 
micromicrofarads. A transistor input capacitance in excess of 20 
micromicrofarads necessitates an auxiliary reactance of opposite sign 
across the meter terminals, shown as L201 in Figure 1, to enable off- 
setting the meter reference balance from zero on the capacitance scale. 



Unilateralization Impedance 

The series resistive and capacitive components of the unilateraliza- 
tion impedance are measured using a specially designed test set in 
conjunction with a signal generator, a receiver, and an oscilloscope. 
The test set, shown in Figures 3 and 4, is also used in measurement 
of transconductance and current gain. 




Fig. 2 — Input and output impedance adaptors. 
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Fig. 3 — Front view of transistor test set. 



In the unilateralization measurement, a CW signal generator, set 
to about 100 millivolts output, is coupled via a relatively high impedance 
to the base of the transistor under test, as shown in Figure 5. The 
emitter is common. The collector load consists of a tightly coupled 
phase-reversing transformer tuned by the capacitance of a reverse- 
biased junction diode (the collector junction of a 2N139 transistor). 
The reversed-phase section of the transformer is coupled back to the 
transistor base through a calibrated, adjustable, series E-C unilateral- 
izing impedance. The tuning of the collector circuit is varied back and 
forth past the fixed frequency of the input signal by varying the 
magnitude of the bias applied to the diode, and thus its capacitance, 3 




Fig. 4 — Rear view of transistor test set. 



3 J. H. O'Connell and L. J. Giacoletto, "A Variable Capacitance Ger- 
manium Junction Diode for UHF," RCA Review, Vol. XVII, p. 68, March, 
1956. 
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at a GO-cycle sweep rate. The resulting variation of collector load 
impedance produces a simultaneous variation in base input impedance 
if the circuit is imperfectly unilateralized. Such impedance variation 
results in a variation in signal voltage at the base. The base voltage 
is monitored with a communications-type receiver. The wave-form of 
the receiver a-f output is a replica of the detected signal envelope, and 
is observed on an oscilloscope, using a 60-cycle, variable phase, sine- 
wave sweep. Figure 6 shows waveforms corresponding to various 
degrees of unilateralization. In these oscillograms, the forward and 
reverse traces are not quite superimposed because of low frequency 
phase shift. 
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Fig. 5 — Unilateralization test circuit. 



Transconductance 

Transconductance may be measured at a low signal frequency using 
the simple circuit shown in Figure 7. An adaptor incorporating this 
circuit is shown in Figure 8. In Figure 7, the input and output 
resistance values indicated have been chosen so that the transconduct- 
ance, in mhos, is simply the ratio of v 2 to %. The input voltage, v x , 
should not be so high as to overdrive the transistor; a maximum value 
of 0.1 volt is suggested. An audio signal generator set to about 2,000 
cycles is considered a suitable signal source. Input and output voltages, 
v x and v 2 , may be measured successively, using any of several commer- 
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Fig. 6 — Unilateralization test oscillograms. 

cially available multi-range electronic voltmeters. Alternatively, the 
voltage at v 2 can be set to a suitable indicated value such as 0.04 volt, 
the voltmeter transferred to v u and the v 2 reading reproduced at v x 
by adjustment of a calibrated attenuator on the signal generator. The 
transconductance is then the voltage ratio indicated by the difference 
in attenuator settings, both of which should be at least 20 decibels 
below maximum attenuator output to avoid possibility of differential 
loading error. 
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Fig. 7 — Audio frequency transconductance test circuit. 
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Pig. 8 — Audio frequency transconductance test adaptor. 

At high frequencies, transconductance is measured on a relative 
basis as a function of frequency. The circuit used is shown in Figure 9. 
This circuit is incorporated in the test set of Figure 3. 

A 30-millivolt, 90 per cent amplitude-modulated signal is connected 
to the input jack. The input is terminated by the bolometer, RG, and 
the voltage divider R7, R8, R9. The bolometer is an envelope detector, 
used to monitor the amplitude of the signal at the input to the voltage 
divider. The assumption is made that the input impedance of the 
transistor is high at all frequencies compared to the G-ohm shunt 
section of the voltage divider, and that the voltage input to the tran- 
sistor is, therefore, independent of frequency. (This is generally the 
case for high-frequency transistors operated at emitter currents less 
than a few milliamperes. Input impedance measurements, as previously 
described, will disclose any instance in which correction factors might 
be needed, and will provide the data for estimating such correction 
factors.) 

A low-impedance path for collector current is provided by bolo- 
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Fig. 9 — HF-VHF transconductance test circuit. 
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meters R2 and R3, which are in parallel for signal frequencies, and 
which serve as an envelope detector for the collector current. The ratio 
of the signal modulation components available from the input and 
output bolometers is the measure of the relative transconductance at 
the frequency of measurement. Such ratios, obtained at frequency 
intervals of several megacycles, and plotted as decibels versus frequency, 
give a complete picture of the transistor transconductance character- 
istic. Such a curve is shown in Figure 10. The low-frequency, flat 
portion of the curve has the absolute transconductance value obtainable 
by measurement at a-f, as previously described. 

Bias is applied sequentially to the input (R6) and output (R2, R3) 
bolometers from a Hewlett Packard 415B Standing-Wave Indicator, 
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Pig. 10 — Relative transconductance versus frequency plot. 

which also provides narrow-band amplification and decibel-level indica- 
tion of the detected signal. Input bolometer R6, a 1/200 ampere Littel- 
fuse, is biased at approximately 4 milliamperes (the "high" current 
position on the HP 415B), at which current it exhibits approximately 
500 ohms resistance. Output bolometers R2 and R3 are 1/100 ampere 
Buss fuses, and are series-biased at about 7 milliamperes (on the 
"high" position), at which current their parallel resistance is about 
56 ohms. Collector signal current is coupled to the bolometers via 
capacitor CI, which is of sufficiently high audio impedance to divert 
any incidental a-f component of collector current from the bolometers 
and output indicator. 

The signal generator used in the transconductance test should 
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preferably be one capable of being modulated at 80-90 per cent, so 
that operation well above the noise level is possible without using 
excessive signal drive on the transistor. Increasing the modulation 
from 30 to 90 per cent improves the signal-to-noise ratio by 4.7 decibels. 
A further improvement of about 2 decibels is obtainable by substituting 
a 400-cycle filter (available from the manufacturer) for the 1,000-cycle 
filter normally supplied in the standing-wave indicator, and using 
400-cycle signal modulation. This latter improvement relates to the 
thermal time-constants of the bolometers. 

Both the magnitude of transconductance and its frequency depend- 
ence change with emitter current, so care must be taken to maintain 
the emitter current at the desired value. 

DC BIASING NOT SHOWN 
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Pig. 11 — HF-VHF current gain test circuit. 



Current Gain 

The method of measurement of relative common-base current gain 
as a function of frequency is quite similar to that described for trans- 
conductance measurement, and requires the same apparatus. The 
current gain test circuit in the test set is shown in Figure 11. This 
differs from the relative transconductance test circuit only in that a 
constant-current signal source is substituted for the constant-voltage 
source. In Figure 11, the 2000-ohm resistance, R10, is assumed to be 
invariant with frequency, and high compared to the transistor input 
impedance, so that with constant signal voltage applied, the emitter 
current is likewise constant. About 0.1 volt of 90 per cent modulated 
signal is used. As in the case of transconductance measurement, the 
transistor input is held constant by monitoring the signal detected by 
the input bolometer, R6. The required precision of measurement 
generally is not obtainable using the signal generator output indication 
for setting the transistor input at the successive frequencies. 
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The alpha-cutoff frequency can be found by a successive approxima- 
tion approach, or by plotting the collector current as a function of 
frequency. Since the current gain is measured on a relative basis, 
factors such as the shunting effect of the resistances through which 
biases are supplied are of no concern as long as such factors are not 
frequency dependent. To obtain the current gain on an absolute basis, 
it may be assumed that the low-frequency current gain is unity. The 
plot of collector current versus frequency then transforms directly to 
a plot of common-base current gain versus frequency. Such a plot, 
obtained by measurement of a 2N384 drift transistor, is shown in 
Figure 12. 



2N384 
V C =-I2V, I E =l.5mo 
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Fig. 12 — Relative current gain versus frequency plot. 



Correction Factor 

When the measurement frequency, the collector load resistance, 
and the capacitances associated with the collector, have values such 
that an appreciable portion of the collector current is by-passed or 
fed back, a correction factor should be added to the indicated current 
gain and transconductance. It is shown in the appendix that in the 
case of grounded-base current gain measurement, the effect of feedback 
can be neglected. The numerical value of the first-order-approximation 
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current-gain correction factor, F, which accounts only for by-passing, 
is 

L 



= V'l + (^CJ-, 



(4) 



l bol 



where i c is the collector signal current, i bol is the signal current deliv- 
ered to the output bolometer, a is the angular frequency of measure- 
ment, R L is the collector load resistance, and C is the collector-plus- 
incidental capacitance shunting the transistor output. This correction 
factor, expressed in decibels, is plotted in Figure 13 for the particular 
load resistance value, 56 ohms, which occurs in the circuit here dis- 
cussed. The indicated number of decibels is added to the decibel-level 













\6 


















\ 0.5db 






























1 — 













MEASUREMENT FREQUENCY - !.IC 

Fig. 13— Correction factor for excess collector capacitance. 

indication of the standing-wave indicator. In entering the plot, the 
test-set incidental capacitance, 3.5 micromicrofarads, is added to the 
transistor collector capacitance. 

It is shown in the appendix that in transconductance measurements, 
feedback to the transistor base lead resistance alters the apparent 
value of collector capacitance as a function of frequency. The equiva- 
lent shunt value of the feedback capacitance, C Ve , is given by 

9m Rib' 
1 +" - 

u 



(5) 
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in which / is the measurement frequency. C,, auiv is added to the inci- 
dental transistor mount capacitance (typically 0.5 to 1.0 micromicro- 
farads) and the test set incidental capacitance (3.5 micromicrofarads) 
to find the ordinate of the correction factor in Figure IS. The factor 
f g in Equation (5) is the true transconductance-cutoff frequency which 
cannot be known until after the correction factor is applied. However, 
the measured transconductance cutoff frequency can be substituted 
with small loss in precision. 

Test Set 

The test set consists of two principal sections: (1) an r-f 
compartment containing the unilateralization, transconductance, and 
current-gain test circuits, and (2) a power supply which provides 




Fig. 15 — Close-up view of r-f box. 

biases for the transistor under test. The complete circuit diagram 
of the test set is shown in Figure 14. 

The r-f compartment is shown in Figure 15. A separate transistor 
socket is used for each type of test. Partitions in the box bridge all 
sockets for input-output isolation, and provide direct low-impedance 
grounds for transistor shield terminals, and for base terminals of 
the transconductance and current-gain sockets. 

The phase-reversing transformer (shown in Figure 15) is wound 
on a ferrite toroid core. For operation at the television intermediate 
frequency, 44 megacycles, the transformer is tuned by the collector, 
sweep diode, and incidental circuit capacitance. By plugging in a shunt 
capacitor (C23), this transformer can be tuned to 24 megacycles. 
Similar components provide for operation at the FM and television 
intercarrier i-f frequencies, 10.7 and 4.5 megacycles, respectively. 
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Fig. 16 — Frequency response characteristic of the transconductance 

test circuit. 

The phase reversing transformers have a collector-to-feedback sec- 
tion voltage ratio of seven to one. Correspondingly, the unilaterali- 
zation adjusting controls, R u and C„, have scales indicating, respectively, 
seven times and one seventh the values at which they are set, so that 
when adjusted for the unilateralized condition, they indicate the values 
of the series unilateralization impedance components as defined earlier. 

The frequency-response characteristic of the transconductance test 
circuit can be checked by plugging a minimum-lead-length 1,000-ohm 
high-frequency type resistor into the socket, between emitter and 
collector terminals. A typical characteristic so obtained is shown in 
Figure 16. 

The current-gain test circuit was checked using a "short" between 
emitter and collector terminals of the socket. The resulting frequency- 
response curve is shown in Figure 17. The increased response at 
high frequencies results from a decrease in impedance of the emitter- 
current-limiting resistor, RIO. This amount of deviation from the 
ideal flat response is not considered significant enough to warrant use 
of a correction factor. 

Power Supply 

The power supply operates from the a-c line. It provides two 
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Fig. 17 — Frequency response characteristic of the current gain 

test circuit. 
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independently adjustable outputs. The collector voltage supply is 
adjustable from 2 to 25 volts and is regulated. The emitter current 
supply, which is unregulated, is adjustable from 0.7 to 5.0 milliamperes. 
A PNP-NPN switch (SI) provides for polarity reversal. With switch 
SI set at either PNP or NPN, the power supply output is available, 
at binding posts on the panel, for biasing a transistor on which input 
or output impedance measurements are being made. In the OFF 
position of this switch, an external bias source can be connected to 
the binding posts to take care of requirements beyond the ranges of 
the internal supplies. 

The collector voltmeter, V c , if returned to common (ground) , would 
read in error by the amount of the IR drop in the 1,000-ohm collector 
r-f isolating resistor (Rl or R17). Since the emitter current closely 
approximates the collector current, a correction is obtainable by re- 
turning the voltmeter to a point similarly "above" common potential 
in the emitter current return path. R2G and R27 provide this com- 
pensation. In PNP tests, the voltage developed across R26 is placed 
in series with the regulator reference voltage. This produces a varia- 




Pig. 18— Transistor common base equivalent circuit. 



tion in output voltage of the collector bias supply, as a function of 
transistor current, of such direction and magnitude as to maintain 
the bias voltage constant at the collector. This convenient feature is 
not available in NPN tests. 

Approximately 9 volts d-c reverse bias for the variable capacitance 
diode, V3, in the unilateralization test circuit, is obtained from the 
input of the collector bias supply at the tap between R23 and R24. The 
a-c sweep voltage for the diode is coupled in through C17 from a 10,000- 
ohm potentiometer (R25) across the 18-volt transformer secondary. 
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Appendix 

The equivalent shunt value of the feedback capacitance which is 
to be used in the correction factor of Equation (4) is obtained from 
an analysis of the transistor common-base equivalent circuit shown in 
Figure 18. The circuit elements of the equivalent circuit are denned 
as follows: C Ve is the sum of the emitter diffusion and transition 
capacitances, C Vo is the collector transition capacitance, C s is the inci- 
dental capacitance which shunts the transistor output and includes 
transistor mount and socket capacitances and circuit capacitances, 
Oc — (<l/kT)I e is the emitter junction conductance, r bb . is the base-lead 
resistance, a is the transistor common-base current gain, and R L is the 
parallel resistance of the two bolometers used to measure collector 
current. 

The two nodal equations below which are used to calculate transistor 
current gain and transconductance are obtained by equating the cur- 
rents entering and leaving nodes b' and c. 

i* = V v (g w + ji>C„. c ) + ai e — ja>C b . c V c , (6) 

ai e = —jmC VB V v + V c UmC b . c + Y L ) . (7) 

Transistor current gain, which is denned as the ratio of collector 
current, i c , to emitter current, i e , can be found by simultaneous solution 
of Equations (6) and (7). If the substitution V c = i c Z L is made, then 
the ratio of i c to i e is 

% (a + 3v>C b . e r w ) 

— = . (8) 

'u 1 + h"C b - c Z L + jo>C,,. l: r hh , 

For most transistors, faC b . e r w « 1 up to the alpha-cutoff fre- 
quency and can be ignored. Equation (8) can now be written 

4 a a 

— = = - -. (9) 

i e 1 + juC }> . c Z L i«>c b . c R L 

1 + 

1 + j*C t R,, 

The relationship between the measured current, i bol , and the collec- 
tor current, i c , is 



% = hoi (1 + faCJt,,). 



(10) 
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If Equation (10) is substituted into Equation (9), 



hoi 



l+j u >R r .(C b . c + C s ) 



(11) 



Vi + ^-Rj/HCs + cv,); 



(12) 



It may be seen from Equation (12) that the ratio of bolometer 
current to input current multiplied by the correction factor, F , gives 
the magnitude of transistor current gain. The equivalent value of 
shunt capacitance to be used in the correction factor is the sum of 
collector capacitance plus incidental capacitance which shunts the tran- 
sistor output. 

Equations (6) and (7) can be simultaneously solved to give the 
ratio of output current, i c , to the input voltage, V in , which is denned 
as transistor transconductance. To eliminate the input current from 
Equations (6) and (7) , the substitutions i e = Y e (v tn — v v ) and ai e = g m 
{v in — v b .) were made. The simultaneous solution of Equations (6) 
and (7) for transistor transconductance is 



ju>C Ve (g m — jo>C h . c ) 



1 + ju>C yc Z L + 



jaC Vc Z L (g m — j<oC b . c ) 



(13) 



Equation (13) can be simplified by assuming that G Vc « C Ve and 
g m » jaC V( . up to the transconductance cutoff frequency, /„. Equation 
(13) simplified is then 



(I, 



1 + faC b r e r bi . 



(14) 



i> { „ 



1 + /«A 



R,. 



1 + ]<»C s R h 



1 + 



.966' + jmCyc 



The numerator of Equation (14) is the intrinsic transconductance 
of the transistor and the denominator represents a measurement error 
caused by a finite value of load resistance which is to be corrected by 
Equation (4) . The ratio of measured current, i boh to input voltage is 
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obtained by substituting the relationship of Equation (10) into Equa- 
tion (14) : 



1 + i»>C Vl .r,, 



1 + jo>R L 



C s + C, yc 1 + 



1 + josC Ve r b 



(15) 



— 



1 + ^Rr 2 



c s + c h . c I 1 




(16) 



u 



From Equation (16) the value of shunt capacitance to be used in 
the correction factor is 



C = C. + C b . c 1 + 



Sm r » 



(17) 



1 + 



/ \ 2 



and Equation (16) multiplied by the correction factor will give the 
ratio of measured output current to input voltage equal to the tran- 
sistor transconductance. 



GENERATION OF SECOND HARMONIC IN A 
VELOCITY-MODULATED ELECTRON BEAM 
OF FINITE DIAMETER* 

By 

F. Paschke 

RCA Laboratories, 
Princeton, N. J. 

Summary— This second-order computation reveals that the amplitudes 
of the second-harmonic current and velocity do not in general vary period- 
ically with distance. The maxima of second-harmonic current are closer 
to the input cavity than the corresponding maxima of the fundamental 
current. At very small beam diameters, the second harmonic exhibits an 
instability and both the current and, velocity amplitudes rise with distance. 
This can be interpreted as a parametric oscillation where the two first-order 
Hahn-Ramo waves serve as low-frequency-pumping sources. The analysis 
agrees, at least qualitatively, with experimental results reported recently 
by Mihran. 

Introduction 

IN klystron frequency multipliers and high-power klystron amplifiers 
a knowledge of harmonic generation is required. One-dimensional 
theories 1 - 2 of velocity-modulated electron beams yield, for small 
signals, the convection-current densities and velocities shown in Figure 
3. The distributions are periodic for both the fundamental frequency 
and the second harmonic. The maximum convection-current density 
of the second harmonic is predicted to appear at the same distance 
from the input cavity as the maximum convection-current density of 
the fundamental. Mihran, 3 however, has shown experimentally and by 
machine computations on a disc-electron model that for beams of 
finite size 

(1) the second-harmonic current does not vary periodically with 
distance; 

(2) the maxima of the second harmonic current are closer to the 



* Manuscript received September 22, 1958. 

1 D. R. Hamilton, J. K. Knipp, and J. B. H. Kuper, Klystrons and 
Microwave Triodes, McGraw-Hill Book Company, Inc., New York, N. Y., 
1948. 

2 F. Paschke, "On the Nonlinear Behavior of Electron-Beam Devices," 
RCA Review, Vol. 18, p. 221, June, 1957. 

3 T. G. Mihran, "Measurement of Harmonic Currents in a Velocity- 
Modulated Electron Beam," 16th Annual Conference on Electron Tube Re- 
search, Quebec, Canada, June, 1958. 
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input cavity than the corresponding maxima of the fundamental; 

(3) for thin beams, the second-harmonic current grows with 
distance. The growth resembles that occurring in a rippled-beam or 
velocity-jump amplifier. 

Mihran attributed these phenomena to the presence of fringe fields 
of the beam. In an earlier paper 2 the nonlinear space-charge-wave 
equation has been developed by successive approximation into three 
simultaneous linear differential equations of second order. These have 
been solved 2 for the one-dimensional case, i.e., for negligible fringe 
fields. In the present paper the first two simultaneous equations will 
be solved with the fringe fields taken into account in order to demon- 
strate analytically the phenomena observed by Mihran. 

Solution of the Nonlinear Space-Charge-Wave Equation 
The following assumptions underlie the analysis: 



Fig. 1 — Current flow through an infinitesimal beam element. 
(1) The electron flow is confined by an infinitely strong magnetic 



(2) The electron and phase velocities are small compared to the 
velocity of light. 

(3) The electric r-f field E in the direction of the stream is 
approximately constant over the cross section of the beam. This is 
tantamount to assuming that the transverse electric field varies linearly 
across the beam. 

(4) The effects of the potential depression caused by space charge 
across the beam and the thermal velocity distribution are negligible. 

Figure 1 shows the current flow through an infinitesimal element 
of the beam; E is the axial electric field strength, i is the convection- 
current density, and S is the area of the beam cross section. The 




field. 
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displacement current per unit length M which flows away from the 
beam accounts for the fringe field. This field is particularly important 
in traveling- wave amplifiers, for it establishes the coupling between 
beam and delay line. According to the continuity equation of Maxwell's 
field theory, the total current flux into a given volume is zero. Thus 
for the infinitesimal beam element, 

di d 2 E M 

+ e 0 + — = 0. (1) 

dx dxdt S 

As a variable the "displacement," a scalar quantity defined by 

y = Cl) E + m (2) 
is used where m represents the fringe field and is given by 

3 2 m M 



(3) 



3a;3£ 



The fields are approximated by first- and second-order perturbations 
which are denoted by the subscripts 1 and 2. The subscript 0 refers 
to unperturbed quantities. The perturbations satisfy the following 
system of equations (taken from Reference 2) : 

3ty, 2 3% 1 3% 

- + + + /V (2/1 — m a ) = 0, (4a) 

3.r 2 v 0 dxdt v 0 2 dt 2 

d-y, 2 d~y n 1 d 2 y, 

+ - + PpHv-.-m,) (4b) 

3.T- v u dxdt v 0 2 dt 2 

2 / Zy x 2> s y x 3 2 2/ t \ 2 / dy x d 2 y x 

+ - 



V'o V 3t dxdt dx dt 2 / i 0 \ dt dx 2 dx 
3r) Zy x 



dxdt I 



( 0 v„ 2 dx 



(Vi — Vh). 



The velocities are given by 



1 / 32/i a 2/i , 

v 1 — + v 0 ) , (5a) 

p 0 \ dt dx 
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1 / dy 2 dy 2 Zy 1 \ 

v 2 = h v 0 ■ h v-l I , (5b) 

p 0 \ Zt dx Zx / 

and the current densities by 

i\~ , (6a) 

dt 



Zy 2 

i 2 = • (6b) 

dt 

ft p is the plasma wave number of the infinitely extended beam, and 
p 0 the d-c space-charge density. The first-order perturbations deter- 
mine the fundamental components and the second-order perturbations 
the second-harmonic components of the modulation. Although the 
theory is applicable only in cases where the transverse electric field 
strength varies approximately linearly across the beam, the analysis 
can be generalized as follows: The first-order exact solutions, the 
Hahn-Ramo waves, are well known. The propagation constants of 
these waves are given by 

P=P.±pM p p (7) 

where /?„ is the beam wave number and 7>(w) the plasma-frequency- 
reduction factor. The latter is a function of the frequency and of the 
geometry- 4 The fringe-field parameter m 1 in Equation (4a) is now 
determined so that the solution of Equation (4a) coincides with the 
exact solution, Equation (7). Thus 

Vi — m i = P 2 M Vi ■ (8) 

Similarly, in Equation (4b) for the second harmonic, 

y 2 — m 2 = p-(2u>) y 2 . (9) 

It is expected that this generalization predicts the dispersion of the 
second harmonic and therefore the field distribution along the axis 
satisfactorily. However, one can not calculate the field distribution 
across the beam by this method. 

To obtain the proper solutions of Equations (4), it is assumed that 



4 G. M. Branch and T. G. Mihran, "Plasma Frequency Reduction Fac- 
tors in Electron Beams," Trans. I.R.E., PGED, Vol. 2, p. 3, April, 1955. 
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the beam is modulated at x — 0 by a double layer with the small voltage 
Y, sin at. If the further assumption is made that 



Up 

— « 1, (10) 



one can neglect the excitation of harmonic velocities by the double 
layer. 2 Then the velocity modulation produced by the double layer at 
x = 0 is (77 Vi/vo) sin *>*• Tne proper solution of Equation (4a) is 

to Vi 

y l= — ■ ■ sin (3 p p (0) x sin (at — /3,.a;) , (11) 

where V 0 is the d-c velocity of the beam in volts. With the aid of 
Equations (5a) and (6a) 

Vi 

v 1 = v 0 cos p p p (<») x sin (mt — /3 e x) (12) 

2V 0 

V 1 

— sin f3 p p(u>)x cos (a>t — fS e x). (13) 

w„p(<o) 2F 0 

With the knowledge of the first-order perturbations Equation (4b) 
can be solved. The boundary condition of no second-harmonic modu- 
lation at x = 0 is met by the solution 



2/3 = — V 



(2<„) \4V„ 



3 

1— - cos 2j8j,p (01) a; 

4^-1 

(14) 



4(1 - I 2 ) 

H ■ cos /3,,?> ( 2<d) .t sin 2 (<«? — /3 c ..c) . 

4£ 2 - 1 

This solution is valid only under the Condition (10). The parameter 
£ is defined by 

p M 

£ = . (15) 

p(2o>) 

£_ has been evaluated from the reduction factors given in Reference 
(4) for a pencil beam with radius b inside of a drift tube of radius a. 
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The result is shown in Figure 2. | approaches unity for extremely 
thick beams where practically no fringe fields are present and ap- 
proaches Y 2 for very thin beams where the fringe fields are strong. 
From Figure 2 it can be seen that for practical beam diameters the 
assumption of one-dimensional fields ($=l) is unrealistic. With the 
aid of Equations (5b) and (6b), the second-harmonic velocity and 
convection-current density are found to be 



v., = v 




sin 2ft p p(u)x 



(16) 



2— Fringe-field parameter, I, versus normalized beam diameter. 



l., = i„ 



1 



+ 



4(1 — | 2 ) 



cos 2/? p p (<u) a; 



4I 2 



■ COS P v p (2a) X 



(17) 



cos 2 (mi — /? t .a;) . 



It can be seen from Equations (15), (16), and (17) that the velocity 
and current-density amplitudes are periodic functions of distance only 
if | happens to be a rational number. Even if the pattern is periodic, 
it is not as simple as that of the one-dimensional theory where | = 1. 



GENERATION OF SECOND HARMONIC 



G23 



For very thin beams p (w) is approximately a linear function of fre- 
quency (unless b/a = 0) and £ = %. Then the inhomogeneous "driving" 
term on the right hand side of Equation (4b) contains a function which 
is a solution of the homogeneous differential equation. Consequently 
there is a resonance, and Equations (14), (16), and (17) become 

/ v x v 

y., - — i H [1 — cos 2(3 p p (o>) x 

+ 3/?„t> U) ^ sin 2^ p p («.) .t] sin 2 M — /3 c x) , (18) 
/ Vi \ 2 

d., = v 0 ■ [sin 2p p p (<o) a; 

o) p p(o.) \ 87 0 / 

+ 6/3 p p (<o) a- cos 2fi p p (o>) a] sin 2 M — /? e a;) , (19) 

2 M 2 / Vi \ 2 
i n = i 0 [1 — cos2^ p p(«) a; 

+ 3/? p p («) a; sin 2/JpP («>) x] cos 2 («t — fox) . (20) 

Thus for very thin beams where | = x /2, the harmonic amplitudes 
increase with distance. The increase is predicted to be linear. How- 
ever, perturbation terms of higher than second order in V x /V 0 which 
have been neglected in the present analysis contain amplitude functions 
of higher than first power in .r. Thus the series whose first two terms 
are given by Equations (11) -(13) and (18) -(20) converge satis- 
factorily only for distances relatively close to the excitation plane, 
x — 0. It may very well be that the actual increase in harmonic ampli- 
tudes is nonlinear. From the following, one can understand this 
instability physically. The second-order nonlinearity of the beam pro- 
duces second-harmonic waves with twice the propagation constants 
of the first-order waves (Equation (7) ) : 

j8 = 2j8,±2p(«.) fl p . (21) 

The waves described by Equations (7) and (21) have equal phase 
velocities and therefore travel in synchronism. In very thin beams 
where | = the second-harmonic waves are as strongly supported 
by the beam as the first-order waves, because the propagation con- 
stants from Equation (21) are identical to those from Equation (7) 
for the frequency 2w. Thus there is a situation where both the funda- 
mental and second-harmonic frequencies are supported by the electron 
beam and where the waves associated with the two harmonic frequen- 
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cies travel in synchronism. Both conditions are required to be present 
if the instability is interpreted as parametric oscillation as suggested 
by Kompfner. 5 In an ordinary traveling-wave-type parametric amplifier 
or oscillator, 6 energy is injected at a pump frequency o> pp and converted 
into energy at a signal frequency «., and an idling frequency w 4 . That 
the process can be reversed with energy injected at u> s and a { and con- 
verted into energy at m pp is compatible with the conservation of energy.* 
The following relations must hold for the frequencies and associated 
propagation constants: 5 

<"iii> — °>i + <»*< (22) 
Ppp = ft + P, (23) 

These conditions are met by the waves described by Equations (7) 
and (21) if <o s = to and o); = (i). Thus the instability of the second 
harmonic can indeed be interpreted as parametric oscillation at the 
frequency 2o> with low-frequency pumping at m. The threshold of this 
oscillation is zero because the resonant structure, the electron beam, is 
lossless. It is expected that instabilities occur also for the higher 
harmonics as soon as 

p(no>) = np(oi) (24) 

becomes a valid approximation. 

Figures 3 to 5 illustrate how the fundamental and harmonic ampli- 
tudes vary along the beam for various values of £, that is, for various 
geometries. Figure 3 shows the familiar pattern of the one-dimensional 
theory, f = 1. Both the velocity and convection-current density are 
periodic functions of distance. However, for £ = 1/V2. which cor- 
responds to a beam with f3 c b = 1.2 and b/a = 2/3, the amplitudes of 
neither the velocity nor the convection-current density are periodic 
functions of distance. This is shown in Figure 4. It is noteworthy 
that the maxima of second-harmonic convection-current density are 
closer to the input cavity than the corresponding maxima of the funda- 
mental convection-current density. Figure 5 shows the instability 
phenomenon which occurs in very thin beams (| = %). It should be 
mentioned that the amplitudes in Figures 3-5 are not drawn in correct 

6 R. Kompfner, Comments to Reference (3), 10th Annual Conference 
on Electron Tube Research, Quebec, Canada, June, 1958. 

0 P. K. Tien and H. Suhl, "A Traveling-Wave Ferromagnetic Ampli- 
fier,'^ Proc. I.R.E., Vol. 46, p. 700, April, 1958. 

* The power relations cannot be given here since a fourth-order theory 
would be required. 
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relation to each other. In Figures 3 and 4 the amplitudes are shown 
in relation to their maximum possible value. The second-harmonic 
velocity amplitude shown in Figure 5a is the function in the paren- 
thesis of Equation (19) divided by 50. The amplitude of the second- 
harmonic current density shown in Figure 5b is the function in the 
parenthesis of Equation (20) divided by 20. 




Fig. 3 — Amplitude functions of the fundamental (dashed) and second- 
harmonic (solid) velocities (a) and convection-current densities (b) versus 
distance from the input cavity measured in plasma-phase angle for infinite 
beam diameter (one-dimensional theory, 1 = 1). 
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Conclusion 

The generation of second harmonic in a velocity-modulated electron 
beam has been analyzed. The fringe field which is neglected in one- 
dimensional theories, has thereby been taken into account. The phe- 
nomena observed experimentally by Mihran 3 have been demonstrated 
analytically. 





Pig. 4 — Amplitude functions of the fundamental (dashed) and second- 
harmonic (solid) velocities (a) and convection-current densities (b) versus 
distance from the input cavity measured in plasma-phase angle for b/a = 2/3 
and /3 c6 = 1.2 (£ = 1/V2). 
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Fig. 5 — Amplitude functions of the fundamental (dashed) and second- 
harmonic (solid) velocities (a) and convection-current densities (b) versus 
distance from the input cavity measured in plasma-phase angle for thin 

beams (£ = 1/2). 
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NOTES ON ERROR-CORRECTING TECHNIQUES 
I. EFFICIENCY . OF SINGLE-ERROR-CORRECTING 
CODES WITH A CONSTANT BIT RATE OF 
TRANSMISSION* 

By 

Jacques Dutka* 

RCA Defense Electronic Products, 
New York, N. Y. 

Summary — Under the stated assumptions, the efficiency of single-error- 
correcting codes for binary data compared with unprotected codes is 
investigated. In particular, the range of values of the bit error probability 
for which the former is advantageous, probability-wise, is determined as 
well as the values for which there is maximum relative efficiency. The 
relative power gain obtained by using single-error-corrccting codes is also 
determined. 

Introduction 

IN EECENT years, methods for constructing error-detecting and 
error-correcting binary codes have been developed by R. W. Ham- 
ming and others, 1 ' 2 but comparatively little has been published on 
the efficiency of such codes compared with unprotected codes under a 
variety of conditions. Such an analysis would have to consider such 
comparisons as the bandwidth, the rate of information transfer, the 
relative power gain (or loss), the feasibility of mechanization, etc. 
This paper considers some limited questions in this area. 

Assumptions 

J. A Bernoulli distribution for bit errors is assumed. This means 
that the occurrence of an error in one bit is independent of the occur- 
rence of an error in any other bit and that the probability of an error 
per bit is a constant, p. 

2. Errors due to fading are not considered since there is little 
analytical information available on their occurrence or on the distri- 
bution of their amplitudes and durations. 

* Manuscript received August 27, 1958. 

t Dr. Dutka is also associated with Columbia University, New York, 
N. Y. 

1 M. J. E. Golay, "Notes on Digital Coding," 1'roc. I.R.E., Vol. 37, p. 657, 
June, 1949. 

2 R. W. Hamming, "Error Detecting and Error Correcting Codes," Bell 
Sys. Tech. Jour., Vol. 29, p. 147, April, 1950. 
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3. The time of transmission of a bit is assumed constant so that the 
time required to send a word is proportional to the length of the word. 
This insures that no additional bandwidth is required in transmitting 
a message in a protected code as compared with transmitting the 
message in an unprotected code. (A protected code is one in which 
redundant bits are added to an m-bit word before transmission so 
that errors which are in the received word may be detected or cor- 
rected.) 

Efficiency of Single-Error-Correcting Codes 

Let p denote the probability that a received bit is in error and let 
7 = 1 — p. The probability that there is at least one error in a word 
m bits in length is P 0 = 1 — q m , and this is the probability that a re- 
ceived word in an unprotected code is in error. In a single-error- 
correcting code, k redundant bits are added to the m information bits 
to produce an w-bit word. Here k is the smallest integer such that 
1 + n ^ 2 k . For an w-bit received word in a single-error-correcting 
code, the probability that there will be at least one error (after single- 
error correction is applied) is P 1 = 1 — q" — npq"— 1 . 

It is now of interest to determine the range of values of the bit 
error probability for which a single-error-correcting code is advan- 
tageous as compared with an unprotected code. More precisely, it is 
desired to determine the values of p for which P 1 P 0 . 

The answer may be summarized as follows : 

P\ — Po f or 0 — V — 1 — a where a depends on n and k, and 
1/2 — a < n (k — l)/k (n — 1). Moreover a — 1/2 if and only if 
2 k — n + 1. Single-error-correction coding is most efficient when 
p = p' = (n — k) /k (n — 1). 

Proof : 

To show that P x — P 0 for 0 — p — 1 — a, it is sufficient to prove 
that the polynomial f(q) = (1 — n) q lc + nq 1 *- 1 — 1 is non-negative 
for a ^ q 1 where a is the smallest positive zero of f(q). For if we 
have 

(1 — n) q k + nq''- 1 — 1^0 

for a — q — 1, then on multiplying this inequality by — q"> and re- 
arranging, 

(n — 1) q" — nq"- 1 — q m 

I _ qn _ npqn-l ^\ — qm 



or 



Pi^Po- 
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It remains to show that /(<?) ^=0 for a ^ q ^ 1. By Descartes' 
Rule of Signs, the polynomial f(q) has k zeros of which two are positive 
(one of the zeros is q = 1), and the remaining zeros are complex ex- 
cept when k is odd, in which case there is one zero of f(q) which is 
negative. An examination of the derivative f'(q) = q k ~" [kq (1 — n) 
+ n (k — 1)] shows that there are critical points when q = 0 and 
q~q' — n {k — l)/k (n — 1). For q s= 0, f (q) is an increasing func- 
tion until q = q' where f(q) has a relative maximum, and f (q) is a 
decreasing function for q > q'. Since /(l) = 0, it follows that f(q') 
is positive. Also /(1/2) = (1 + n — 2 k )2~ h is negative or zero ac- 
cording as 2 k > n + 1 or 2" = n + 1. If /(1/2 < 0, then since 
f(q') > 0 and /(</) is continuous, there exists an a, 1/2 < a < q', such 
that a is a zero of f(q) . If /(l, '2) = 0, then 2* = n + 1 and a = 1/2. 
Since /((/) has a maximum when q = 1 — p = n (k — 1) /k (n — 1), 
it follows that single-error correction is most efficient (in the sense 
that the ratio (1 — Pj)/(1 — P 0 ) is a maximum (when p = (n 
-k)/k(n-l). 

The statement is sometimes made with respect to single-error- 
correcting codes that they are useful only when the bit error probability 
is small and that at medium or high bit error probability rates, they 
offer no advantage over unprotected codes. In fact the claim has been 
made that under the latter conditions it is better to use unprotected 
codes. The validity of this statement will now be investigated. 

It is sufficient to restrict ourselves to the case where p < 1/2. For 
if the probability of error of a received bit is known to exceed 1/2, the 
bit may be inverted (a "0" is replaced by a "1" and vice versa) and 
the bit error probability becomes less than 1/2. (Of course, if p = 1/2, 
neither protected nor unprotected codes will be of any value.) It is 
easily verified that 

As n co, p' = ( n — fc) /k (w — 1) -> 0. The value of the bit error 
probability for which single-error-correction coding is most efficient 
tends to zero as the ivord length becomes infinite. 

For all n sufficiently large, P 1 ^ P 0 , and thus a single-error-correct- 
ing code is advantageous compared with an unprotected code if only 
the word length is large. More precisely, 

For n ^ 13, P x ^ P 0 for 0 ^ p ^ 1 - a where 1 - a > (1/2) + 
/(l/2)//'(l/2) and f(q) = (1 - n)q" + nq'" 1 - 1. As n -> co, 
(1 — «) -» 1/2. Thus the upper bound of the bit error probability for 
which single-error correction is advantageous tends to 1/2 as the ivord 
length becomes infinite. 

Proof : 

To show that a — > 1/2 as n — * » where a is the smallest positive 
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root of /(</) = 0, it suffices to find a number a x such that % > a — 1/2 
and to verify that a x —> 1/2 as n — > oo. A convenient choice for a x is 
the Newton-Raphson approximation 

1 /(1/2) 1 1+71 — 2" 



2 /'(1/2) 2 2 [w (fc — 2) + fc] 

Since log 2 (w + 1) — < 1 + log 2 (« + 1), a t — > 1/2 as n — > to. It 
remains to show that a x > a for w — 13. Subtracting a from <x 1( 



a, — a- 



(1/2 -a) f (1/2) -/(1/2) 
f (1/2) 



(1) 



where /'(1/2) > 0 and it must be proved that the numerator on the 
right is positive. On expanding /(«) by Taylor's formula about the 
point q — 1/2, 

(a -1/2) 2 

/(«) =/(l/2) + (a -1/2) /' (1/2) +- -/"(|) 



2! 



where 1/2 < $ < a. Since /(a) = 0, 



(a -1/2) ^ 

(1/2 -«) /'(1/2) -fa/2) =- — rca 

2! 



(2) 



The left-hand side of Equation (2) is the same as the numerator 
of the right-hand side of Equation (1). The problem thus reduces to 
showing that /"(<?) is positive in the interval [1/2 to a] for n — 13. 
Now q" = n(k — 2)/k(n — 1) where /" (q") = 0 and /" (q) > 0 for 
0 < q < q". To show that a < q" (and thus that /" (q) > 0 in the 
interval [1/2 to «] it will now be proved that f(q") > 0. 

f(q") = (?")* [(1-n) +«( 9 ")- 1 - (<7")- fc J 
2(n— 1) /&(« — !) 



/(<7") =■ 



= (tf") fc 
2 (?")" (n-1) 



n (Zc 



n(fc — 2) 



fc 2 



2 (fc-2) 



1 \ 



1 - 



(3) 



By direct substitution, it is found that the bracket on the right-hand 
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side of Equation (3) is positive for k = 4, with 13 ^ n ^ 15 and also 
for k = 5 with 16 ^ n ^ 19. Now since (1 — l/n)*- 1 < 1, 

/(<?")> n ■ 1+ ) . (4) 

n(&-2) L 2(i-2) fc-2/ J 

By substitution, it is found that the bracket on the right-hand side 
of Inequality (4) is positive for 20 ^ n < 64 (and * denned by 
log 2 (1 + n) =s/c < 1 + log 2 (1 + n)). Finally, since [1 + 2/(fc -2)]*-s 
< e 2 for all & > 2, 

2 (<f) fc (»-l) T ft=e* - 

/(<?")> n 



n(k-2) L 2 (A; — 2) _ 

Now 2''-i < n < 2" for all k, (n cannot be an integral power of 2), and 
it is readily verified that for k ^ 7, 2*-* - A; 2 e 2 /2 (fc - 2) > 0. Hence, 
on summarizing the above, we see that f(q") >/(<*) =0 for n 13' 
and since f(q) is an increasing function for 0 < q < q', it follows that 
q" > a and f"(q) > 0 in the interval [1/2, a]. Thence we see that 
% > a for » ^ 13, Q.E.D. 

Table /—Single-Error-Correcting Code Compared with Unprotected Code 



■m 



p'= (n — k)/k(n—l) 



20 i ok - 1923 - 4845 

40 r 5c - 1667 - 4643 

60 S 4 £ -1482 .4595 

100 ? in? - 1329 • 460G 

lm I 107 .1348 .4827 

r ? -1667 .3596 

4 4 9 -1563 .3664 

5 H - 1500 -3906 

191 o J? - 129 5 .4273 

121 8 129 .1182 4390 



Table I summarizes the results obtained above as applied to som< 
words of convenient length. The first part of the table applies to words 
in which the number of information bits is m = 10, 20, 40, 60 80 anc 
100. The second part of the table applies to words in which the length 
is n — 5, 9, 17, 33, 65 and 129, where n is of the form 2*-» + 1 This 
form represents the cases which are most unfavorable to single-error 
correction. For since 2'*-i < « < 2 * if we choose n = 2*-» + 1, then 
m/n, the ratio of the number of information bits to the total number of 
bits, is a minimum. The column marked p' denotes the value of the bit 
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Table II — Probability of Error for Unprotected Word 



v m 

P x 


5 


10 


20 


40 


OU 






.0001 


.00049 


.00100 


.00200 


.00399 


.00598 


.00797 


.00995 


.0005 


.00250 


.00499 


.00995 


.01981 


.02956 


.03922 


.04878 


.001 


.00499 


.00996 


.01981 


.03923 


.05826 


.07692 


.09521 


.0025 


.01244 


.02472 


.04883 


.09528 


.13945 


.18147 


.22144 


.005 


.02475 


.04889 


.09539 


.18168 


.25974 


.33035 


.39423 


.007 


.03451 


.06784 


.13107 


.24496 


.34392 


.42991 


.50464 


.01 


.04901 


.09562 


.18209 


.33103 


.45284 


.55248 


.63397 



10 




p (Probability of error per bit ) 



Fig. 1 — Probability of undetected error per word versus probability of 
error per bit for unprotected code. 



Table III — Probability of Error for Protected Word 



p N 


5 


10 


20 


40 


60 


80 


100 


.0001 


.00000 


.00000 


.00000 


.00001 


.00002 


.00004 


.00006 


.0005 


.00001 


.00002 


.00008 


.00026 


.00054 


.00091 


.00137 


.001 


.00003 


.00009 


.00030 


.00101 


.00212 


.00354 


.00529 


.0025 


.00022 


.00056 


.00180 


.00601 


.01241 


.02032 


.02981 


.005 


.00088 


.00219 


.00695 


.02237 


.04462 


.07077 


.10063 


.007 


.00171 


.00422 


.01321 


.04139 


.08040 


.12442 


.17269 


.01 


.00344 


.00840 


.02576 


.07753 


.14487 


.21632 


.29010 
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error probability for which single-error correction is most efficient, and 
the column marked 1 — a denotes the value of the bit error probability 
which must not be exceeded for single-error correction to be advan- 
tageously compared with an unprotected code. 

Table II gives the probability P 0 versus p and Table III gives the 
probability P x versus p for words containing m = 5, 10, 20, 40, 60, 80 
and 100 information bits. Figures 1 and 2 denote the corresponding 
graphs. It will be noted that both sets of curves are approximately 
linear in form when drawn on log log graph paper. Table IV gives 

I.Oi 1 1 1 1 — | — , 




p ( Probability of error per bit) 

Fig. 2 — Probability of undetected error per word versus probability of 
error per bit for protected code. 

the number of information bits versus the total number of bits per 
word for 3 ^ 100. Figure 3 shows a plot of Ic versus n and Figure 4 
shows a plot of the ratio m/n versus n for the same range of n. Table 
V gives the probability of undetected error for a word when single- 
error correction is applied to a word of 52 information bits, to two 
subwords of 26 information bits each, and to four subwords of 13 
information bits each. Figure 5 shows the corresponding graphs. It 
appears that comparatively little protection is gained by breaking up 
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log 2 (n+ I ) < k < I * log g ( n + I ) 



20 



40 60 80 

n ( total number of bits per word ) 



100 



Pig. 3 — Number of redundancy bits versus total number of bits per word. 



« .8 



n = m + k 
log 2 (n-H) 


< k < 1 + log 2 


n + 1) 















































E — c 0 20 40 60 80 100 

n { total number of bits ) 

Fig. 4 — Ratio of information bits to total bits per word versus total 

number of bits. 
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Table IV — Information Bits Versus Total Number of Bits Per Word 



k n m/n 



OA I 3 - 3333 

^■f, 3 5-7 .4000 -.5714 

.I'll 4 9-15 .5556 -.7333 

H'tl I 17 - S1 .7059 -.8387 

M S 33 - 63 .8182 -.9048 

58 - 9 ° 7 65-100 .8923 -.9300 

n = to + fc, log, (w + 1) & < 1 + ] og2 ( W + !) 

TViMe V — Probability of Undetected Error for Protected Word 



Two 26-bit words with Four 13-bit words with 

(Jne 5Z-bit word with 5 redundancy bits 5 redundancy bits 

V 6 redundancy bits each each 

.0001 -00002 .00001 .00001 

.0005 .00041 .00023 .00027 

•001 .00159 .00091 00061 

.0025 .00941 .00553 00372 

•005 .03435 .02100 01443 

.007 .06260 .03944 ; 0 2755 

•01 .11467 .07828 



.05390 



the 52-bit word into protected subwords when the number of extra 
redundancy bits required is taken into account. 

Assume that a mark is represented by a d-c pulse of amplitude E 
volts and a space by a pulse of zero volts amplitude. It is also supposed 
that the transmitted signal is corrupted by additive Gaussian noise of 
zero mean and root mean square value v r . An error in a received bit 
occurs when a mark (space) is sent and the noise is of negative (posi- 
tive) polarity of magnitude £7/2. If we suppose that the probability 
of transmitting a mark (space) is 1/2, then the probability of error 
in a received bit is 



-dt. 



1 r-urn, e-'va 1 ^ x e _ (!/2 r<n 2/2 
P = ~ / — — dt + - ■ / —dt= / 

2 J -« V2. 2 J vnTt VStt" A /2 , v 

Let p l denote the effective bit error probability for an m bit word 
denned by the relation 1 - (1 - Pl )« = P} . That is, Pl is what the 
bit error probability would have to be in order to yield the same word 
error probability for an m-bit word as the word error probability for 
an w-bit word which includes single-error correction and in which the 
bit error probability is p. Let q t = 1 - Pl and let E u the effective pulse 
amplitude, be defined by 



■ ■ E x /2r,. V2tT 
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The quantity 20 \og 10 E^'E is a measure, in decibels, of the power 
gain (loss) obtained by using a single-error-correcting code (and 
more generally of a protected code). This quantity versus p and m 
is given in Table VI and the results are shown graphically in Figure 6. 



i.o 



.- .01 



W| = 1 word of 
W 2 =2 words of 
W4 = 4 words of 


52 info 
26 info 
3 info 


motion bi 
motion bi 
motion bi 


s plus 
ts plus 
t s plus 


6 p 
5 p 
5p 

/ 


ori t 
ari t 
arit 

/ 


y bits- 

y bits each. 

y bits each. 




/ 






'/ 








/ 


/ / r 











.0001 5 .001 25 5 7 .01 

p (Probability of error per bit) 



Pig. 5 — Probability of undetected error per word versus probability of 
error per bit for Single Error-Correcting Code applied to subwords. 



It would be of interest to obtain the quantity 20 log 10 EJE as an 
explicit function of p 1 and p, and, as an initial step, E x as an explicit 
function of p t . It does not appear likely that a usable exact form for 
E x as a function of p x can be obtained by inverting the above integral. 
Instead we shall obtain upper and lower bounds for E x as a function 
of p u and for the quantity 20 log 10 EJE. We have 



2 iv 



In 4p 1 q 1 ^ E x 2v r V — 2 In Ap x q x 



(5) 



and 
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Ei \ /ln4p 1 <? 1 
20 log,,, I -101og 10 



E 



In 4pq 



10 1og 10 — = 1.049 db. (6) 



Proof: 

Consider the surface z = (1/2tt) exp {— (x- + y-) /2). The volume 
under this surface bounded by the right cylinder whose base is the 
square —a =s x, y ^ a is 



a *s — a 



z dxdy. 



Table VI — Decibel Gain In Use of Protected Code 



m 


V 


E/2v r 


Ex/tv, 


20 log M Ei/E 


.5 


.0001 


3.7190 


5.2626 


3.0157 




.0005 


3.2905 


4.6339 


2.974 




.001 


3.0902 


4.3387 


2.947 




.005 


2.5758 


3.5759 


2.850 




.01 


2.3263 


3.1991 


2.767 


10 


.0001 


3.7190 


5.2148 


2.936 




.0005 


3.2905 


4.5854 


2.882 ' 




.001 


3.0902 


4.2889 


2.847 




.005 


2.5758 


3.5184 


2.708 




.Ul 




0.1387 


2.602 


20 


.0001 


3.7190 


5.1234 


2.782 




.0005 


3.2905 


4.4662 


2.654 




.001 


3.0902 


4.1736 


2.611 




.005 


2.5758 


3.3922 


2.392 




.01 


2.3263 


3.0092 


2.236 


40 


.0001 


3.7190 


5.0204 


2.606 




.0005 


3.2905 


4.3603 


2.445 




.001 


3.0902 


4.0532 


2.356 




.005 


2.5758 


3.2492 


2.017 




.01 


2.3263 


2.8751 


1.840 


CO 


.0001 


3.7190 , 


4.9529 


2.489 




.0005 


3.2905 


4.2889 


2.302 




.001 


3.0902 


3.9738 


2.184 




.005 


2.5758 


3.1714 


1.807 




.01 


2.3263 


2.7944 


1.592 


80 


.0001 


3.7190 


4.8980 


2.392 




.0005 


3.2905 


4.2360 


2.194 




.001 


3.0902 


3.9197 


2.065 




.005 


2.5758 


3.1152 


1.651 




.01 


2.3263 


2.7435 


1.432 


100 


.0001 


3.7190 


4.8560 


2 317 




.0005 


3.2905 


4.1942 


2.107 




.001 


3.0902 


3.8762 


1.969 




.005 


2.5758 


3.0732 


1.534 




.01 


2.3263 


2.7048 


1.309 
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The volume under the surface of a right cylinder whose base is the 
circle x 2 + y- — a- inscribed in the square is 



i r 2 * r 

2* Jo Jo 



e _rV2 rdrdO. 



and the volume under the surface of a right cylinder whose base is 
the circle x 2 + y 2 - Aa?/ir, the area of the circle being equal to that 
of the square, is 



2n- /*2a/\rt 



-r f 

2» Jo Jo 



c —r*/s rdrdQ, 




.0001 



5 .001 5 

p (probability of error per bit) 

Pig, 6 — Decibel gain in using protected code. 



It may be readily shown that 



. 1 f I e-^rdrdO^ 
2" Jo Jo 

r r 1 r r 

_L,.L, zdxdy ^Jo Jo 



e -rV2 rdrdO 
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and on evaluating the integrals on the left and right, 



I zdxdy ^ 1 — e-w* 

a J — a 

Since 

/ / zdxdij = ( / c-' 2 /2 dt/sJTr 

J —a J —a \J —a 



it follows that 



(The relative error involved in approximating the middle term by the 
term on the right never exceeds 0.75 per cent.) Let a = E\/2v r . Then 
the middle term is equal to 1q x 1, and solving these inequalities for 
E t gives Equation (5). 

From these inequalities we may write a similar set involving E 
and p and from these results we find 

/ ln4p l9l \ 4 E t 

101og 10 — -101og 10 — ^201og 10 

\ In 4pq J it E 

and 

E x f In4p,7, \ 4 
20 1og lo — ^ — — +101og 10 — 
E \ \n4pq / tt 

which is equivalent to Equation (6) . 

The practical usability of these results in any particular application 
depends, of course, upon the magnitude of the error which is acceptable 
in the application. 

Conclusions 

At the cost of a reduction in the rate of information transmission 
by a factor m/n and the requirement of a temporary storage to com- 
pensate for the delay, the use of an error-correcting code has, under 
certain conditions, resulted in efficiency gains. 

In this paper we have established: 

(a) the ranges of bit error probability for which single-error 
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correction is advantageous and the dependence of these values 
on the word lengths; 

(b) the word lengths for which the efficiency of single-error-cor- 
recting codes as compared to unprotected codes is maximum ; 

(c) the fact that for a typical case of a 52 information bit word, 
relatively little gain is obtained in applying single-error 
correction to subwords; 

(d) the relative power gain in applying single-error correction 
for typical values of bit error probability and for typical 
word lengths; 

(e) upper and lower bounds in analytical form for the relative 
power gain as functions of the actual and effective bit error 
probabilities. 



AUTOMATIC CONTROL OF VIDEO TAPE 
EQUIPMENT AT NBC, BURBANK* 

By 

Robert W. Byloff 

National Broadcasting; Company, Inc., 
New York. N. Y. 

Summary — The presentation of television programs at the proper local 
time in the different time zones in the United \States has created unusual 
program recording demands upon the network television broadcast indus- 
try. To solve this problem, the National Broadcasting Company has 
developed an integrated automatic system using video tape recorders and 
associated automatic control equipment. 

_ This paper describes the equipment and physical layout of this new 
facility, modification of equipment and methods employed to provide auto- 
matic cueing of video tape recorders, equipment that automatically main- 
tains the television color signal at the required level, and sequencing 
devices that automatically switch the input and output circuits, auto- 
matically stop, start and rewind the video tape recorder machines all to 
provide one-, two-, three- or four-hour program delays or any combination 
thereof. 



Introduction 

V 



I 



^ HERE are four time zones across the United States — Eastern, 
Central, Mountain, and Pacific. In addition, during the summer 
months "daylight-saving" time is in effect in certain areas 
causing a fifth time zone to exist. Thus, there may be a four hour 
time differential between the east and west coasts. 

Television programming is arranged with the time of day in mind ; 
that is, there are morning, afternoon, and evening programs each 
designed to appeal to a particular audience group. It is important, 
therefore, that programs designed to go on the air at about eight 
o'clock in the evening, go on at that time in each part of the country. 
Unless the program is delayed by some means, a program going on 
the air at eight o'clock in the evening in New York would be on at four 
o'clock in the afternoon in Seattle, Washington, during summer months. 
In this case, the program broadcast in Seattle would not reach the 
audience for which it was intended. 

Kinescope recording has been used for a number of years to delay 
programs in the far west, but until video tape recorders with their 
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ability to play back immediately, became available, it was impractical 
to delay programs „ for as little as one hour. Now it is possible to 
provide a network program service in which programs are aired in all 
sections of the country at the same local times. 

Such a service, was inaugurated by the National Broadcasting 
Company (NBC) in April, 1958. By means of video tape machines in- 
stalled in Burbank, California, programs recorded from the eastern 
NBC network are played back to the central and Pacific coast networks 
one and three hours later, respectively. The map of Figure 1 indicates 
the time differential by the use of shaded areas. Clocks in these areas 
show a representative time in each zone. Thus when it is eight o'clock 
in the northeast, it is seven o'clock in areas on eastern standard and 
central daylight times, six o'clock in areas on central standard time, 
five o'clock in areas on mountain standard and Pacific daylight time, 
and four o'clock in a small area in the Pacific northwest which stays 
on Pacific standard time. The complete solution to this problem would 
be to provide services of one-, two-, three-, and four-hour delay to the 
respective areas. However, by providing only a one-hour and three- 
hour time delay, it is possible to provide the entire network with pro- 
grams which are on, in terms of local times, at most only one hour 
different from the local time of origination in the east. In general 
this is being done with the entire NBC program schedule. Figure 1 
shows the time schedule for a typical evening program across the 
country. In terms of local time; it is on at 8:00 PM in New York, 
8:00 PM in Chicago, 7:00 PM in Omaha, 8:00 PM in Los Angeles, 
and 7:00 PM in Seattle. 

Video tape recorders are presently arranged to record and play back 
reels of tape one hour in length. If a one-hour delay is required, it is 
impossible to record for a full hour, since time must be allowed to 
rewind the recorded reel and ready it for playback. To avoid discon- 
tinuities in program material, it was decided to record for half-hour 
periods. Thus, any switching operations could occur during "station 
break" periods. A recorder involved in the delay operation records for 
a half hour, rewinds and waits during the next half hour, plays back 
the one-hour delay during the third half hour, then rewinds and waits 
during the next one and a half hours so that during the seventh half 
hour, it can play back the three-hour delay. It then rewinds and waits 
during the eight half hour to be ready to record again during the 
ninth half hour. Other machines on a staggered basis fill out the sched- 
ule as shown in Figure 2. Note that it takes eight machines using this 
recording-playback format to do a continuous one- and three-hour 
delay. It is normal practice in network broadcasting to record the 
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same program and play it back on two machines simultaneously to 
afford a "protection" copy recording. Thus, if an element in one chain 
of equipment should fail, the program can continue on the protection 
chain. They are normally termed "A" and "B" copies. For complete 
protection, therefore, sixteen machines are required. Fortunately, the 
program schedule is such that continuous recording and playback is 
not required for the entire day, and so it is possible to perform this 
operation with twelve machines which will allow recording and play- 
back on a one- and three-hour delay basis for three hours. If an hour 
is then skipped, the cycle can be repeated. 



Plan of Facility 

The plan of the tape recording plant at Burbank, California, is 
shown in Figure 3. Since color programs occur at random times during 
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Fig. 2 — Record-playback cycle for video tape recorders for one- and three- 
hour delayed broadcasting. 



the schedule, and will increase in frequency, it was decided to "colorize" 
all the machines. The plant contains, at present, four RCA prototype 
color video tape recorders and eight RCA-colorized Ampex tape re- 
corders. Each RCA recorder consists of six standard equipment racks 
and a monitor housing containing color and black-and-white monitors, 
a cathode-ray oscilloscope, VU meter, and a control panel. Each color- 
ized Ampex machine consists of four video racks, a large console and a 
monitor. All of these equipments are installed in one space to gain 
maximum efficiency. It is possible to observe the operation of each 
machine from any spot in the room. The status panel at one end of 
the room shows the status of the machines in the automatic delayed 
broadcast system, and permits rearrangement of the operational se- 
quence at any time. In the quality control room, the quality of the 
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Fig. 3 — Video tape equipment layout. 




Fig. 4 — Video tape recording room 
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Fig. 5 — RCA color video tape recorders. 



delayed playbacks is monitored in color and in black-and-white, and 
the choice of "A" or "B" copy is made. 

The remaining space is used for maintenance, erasing facilities, 
parts storage, and for future expansion. 

The tape room is shown in Figure 4; the RCA recorder facility is 
shown in Figure 5; the Ampex recorder facility is shown in Figure 6; 
and the quality control room is shown in Figure 7. 




Fig. 6 — RCA-colorized Ampex recorders. 
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Video and Audio System 

The video and audio systems are identical in form; the two systems 
are shown by the block diagram of Figure 8 which also shows the 
control system. The twelve video tape recorders are connected between 
input and output switching systems. 

The input switching system, which is used during the recording 
process, is fed by the incoming network line from the east, all the 
Burbank studio outputs, and test signals. These test signals are re- 
corded at the beginning of each tape to act as a check on the operation 
of the machine before it is used on the air. Automatic-gain-control 
amplifiers, both in audio and video, immediately precede the video 




Pig- V — Tape-quality control room. 



tape recorder input to insure that the recorded signal level is as high 
as permissible without overloading; this gives maximum signal-to-noise 
ratio. The a-g-c amplifiers also eliminate signal level variations. 

Each video tape recorder is provided with a color monitor, a black- 
and-white monitor, a cathode-ray oscilloscope, a VU meter, and a loud- 
speaker. The monitors are switched automatically to the input or out- 
put of the machine according to whether the machine is recording or 
playing back. 

The purpose of the output switching system, used during playback, 
is to connect the outputs of the twelve recorders to outgoing lines for 
delayed broadcasting and for playback to studios when integrating 
with other program material. The four lines provided for normal 
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delayed broadcast playback are nominally, the one-hour delay to the 
midwest, the three-hour delay to the Pacific coast, a two-hour delay, 
which is used for special purposes, and a four-hour delay for future 
use. Monitoring on these four lines is available in the quality control 
room, and control of output switching is provided to permit switching 
between "A" and "R" copy tapes. 



Types of Control 

The control system has several operating modes. Controls at each 
machine are provided to permit selection of the control mode desired 
as follows: 
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PRESET OPERATION AT MACHINE 
D3 SEQUENCE 



Fig. 8 — Video and audio systems. 



A. Manual Control at Machine — In this mode, the recorder oper- 
ator can select the input line, the output line, and the recorder operation 
(i.e., record, playback, rewind, fast forward, or stop). This mode is 
used mainly for testing. 

B. Manual Control in Studios — In this mode, the machine oper- 
ation is controlled by any one of the Burbank studios. The studio to 
have control is selected at the machine. 
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C. Delayed Broadcast (D-B) Sequence — In this mode, all machine 
operations including input and output line selection and switching are 
assigned to the automatic D-B sequencer. The recorder enters the 
cyclic operation required for delay broadcasting as outlined above. 

D. Preset Operation at Machines — In this mode, the operator 
may preset the machine to record or play, preselect the record or play- 
back line, and set the time and duration of operation. The machine 
will automatically start, select the proper lines, record or play back 
for the specified period, and then automatically rewind. This mode is 
used for scheduled operations of a noncyclic type. 

Almost all operating situations can be handled by control modes B, 
C, and D thereby eliminating manual switching operations in the tape 
room which was the goal in the design of the control system. 





MONOCHROME SECTION 






1 




CHROMINANCE SECTION 






SYNC SECTION 


COLOR SENSING SECTION 



Fig. 9 — Color video automatic-gain-control amplifier. 

Automatic Control Devices 

In order to simplify operation, certain devices were added to the 
tape recorders. 

A. A-G-C Amplifiers — Automatic-gain-control amplifiers are used 
in the audio channel of the recorders. In addition, color video a-g-c 
amplifiers are used in the video paths. These amplifiers, which were 
developed recently by NBC, perform many functions as shown in Fig- 
ure 9. The incoming video signal is passed through a filter which 
separates luminance and chrominance components. In the luminance 
channel sync is clipped, the base line cleaned up, and the signal con- 
trolled to a particular level. In the chrominance channel the burst is 
measured, and by this means chrominance is held to a predetermined 
level. A sync channel separates sync from the incoming video signal, 
clips and shapes it, and adds it to the luminance and the chrominance 
information at the output of the amplifier, through suitable delay net- 
works. Thus the amplifier automatically sets and maintains separately 
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the levels of luminance, chrominance, and sync of a color signal. The 
color sensing section in the amplifier detects whether the incoming 
signal is color or monochrome, and if monochrome, automatically shuts 
off the chrominance channel, and widens the bandwidth of the lumi- 
nance channel. 

B. Automatic Cueing Equipment — The recorders are arranged so 
that they "cue up" automatically to a predetermined point, thus pro- 
viding a selected length of playing time before the program material. 
The problem was to insure that after a recording was made for delayed 
broadcast operation, the machine would rewind and "cue up" to the 
proper position for playback. 

Equipment for special cue tracks was not immediately available so 
it was necessary to improvise temporarily. The equipment was ar- 
ranged to detect the absence of control track during rewinding and 
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Fig. 10 — Automatic cueing. 



to detect the presence of 1000-cycle audio tone when playing to come 
up to cue. This tone is recorded at the beginning of each recording 
as a test. A diagram of the functioning of the cue system is shown in 
Figure 10. The rewind signal causes the machine to rewind and also 
to connect the "absence of control track detector." When the tape has 
rewound to the point where no control track is detected, a relay closes 
putting the machine in the. fast forward mode. This is done to bring 
the machine to a stop without using the brakes and thereby reduce 
brake wear. After four seconds in the fast forward mode, the machine 
is given a stop signal, and two seconds later a play signal. The 
machine plays until it detects 1000-cycle audio tone on the audio track, 
then stops. The 1000-cycle signal is recorded automatically during 
the recording cycle exactly ninety seconds ahead of program material. 
When cue heads are available, the system will be changed to use 
specially recorded signals on this track for cueing. 
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Fig. 11 — Delayed-broadcast sequence status panel. 



Delayed-Broadcast Sequencer 

Automatic delayed broadcast operation is available through use 
of the D-B sequencer. This sequencer consists of a status panel and 
racks of control equipment. The status panel is shown in the photo- 
graph of Figure 11. The system has two similar parts: one for the 
"A" and another for the "B" copy tapes. Machines are assigned to 




Fig. 12 — Delayed-broadcast sequence control racks. 
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the "A" sequence at the top section of the panel. There are ten so- 
called "slots" for machine assignment, permitting ten machines to be 
assigned, so that continuous one-, two-, three-, and four-hour delays 
can be handled. Machines may be assigned in any order to these delay 
slots. Assignment is made by selecting the desired machine with the 
selector knob at the right side of the panel and then pressing the 
"enable" button under the appropriate slot. The read-outs indicate the 
machine assignments. A similar setup is provided for "B" copy tapes 
in the next section of the panel. Controls at the lower section of the 
panel provide input and output line assignments. These controls per- 
mit any input line to be used for recording and also permit assignment 
of any delay to any output line. The delay modes can also be selected 
at this location. Thus, the system can be set up for any combination 
of one-, two-, three-, and four-hour delays. 

Figure 12 is a photograph of the control equipment racks asso- 
ciated with the D-B system. 

Figure 13 shows the D-B sequence control system, in a very much 
simplified representation, as it would be if set up to record and play 
one- and three-hour delays, "A" copy only. Eight recorders are required 
to provide the complete delay. The heart of the system is the sequence 
stepping switch. This is a multideck-"rotary" switch mechanically 
ganged and electrically pulsed from clock impulses. The switch moves 
one step every half hour. The points around the switch correspond to 
machine positions. The drawing shows recorders 1-8 around the switch. 
Each deck represents a different machine operation. Thus, the first 
deck shown at the top of the drawing is for record. The drawing 
shows that this switch deck is in position to make VTR 1 record. The 
record impulse also sets up the input switching system to select the 
proper input line and connect it to VTR 1. One-half hour later, this 
switch will rotate one step and will cause VTR 2 to record, and the 
input switching system to switch the input line to VTR 2. This will 
continue on succeeding half hours to make all recorders record in 
sequence. The second deck is a rewind deck. It is shown in position 
so that VTR 2 is rewinding. Deck 3 is the three-hour delay play deck. 
Referring to Figure 2, note that while recorder 1 is recording, recorder 
3 is playing the three-hour delay, and recorder 7 is playing the one- 
hour delay. This is the condition shown in Figure 13. Deck 3 also 
sends impulses to the output switching system to set up the proper 
output line to the recorder playing back the three-hour delay. A 
selector switch permits preselection of the line to be used for three- 
hour delay. The other decks in the sequence switch as shown in 
Figure 13 are: deck 4, rewind; deck 5, play two-hour delay; deck 6, 
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Fig-. 13 — Delayed-broadcast sequence system control. 
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rewind; deck 7, play one-hour delay; deck 8, rewind. Note that while 
deck 5 is calling for recorder 5 to play back the two-hour delay, the 
delay mode switch for two-hour delay is open, and therefore the play- 
back does not occur. For the sake of simplicity, the system which 
allows any recorder to be assigned to any position in the sequence, and 
the complete system which permits four-hour delay are not shown. 

Preset Operation 
Figure 14 shows the time preset system with which each recorder 
is separately equipped. It enables an operator to preset recorders to 
record or play from any of the recording or playback lines at specified 
times and for specified intervals. When there is coincidence between 
the clock and the time set into a recorder, a pulse passes to the recorder 
which starts the recorder and switches the lines. At the same time 




Fig. 14 — Time preset system. 

this pulse starts an interval timer, which stops the recorder and re- 
leases the lines when it has completed the interval. 

This preset system supplements the D-B sequencer to provide essen- 
tially automatic operation of the tape equipment used for delayed 
broadcasting by NBC at Burbank. 

Operating Experience 

This equipment has been in operation since April, 1958. At the 
present time over four hundred machine hours per week are being 
logged for actual recording and playback. During this period the 
system has provided NBC with a reliable delayed broadcast operation 
with a minimum of confusion and errors. It has been very important 
in achieving efficiency in this operation. 



POST-INSTALLATION PERFORMANCE TESTS OF 
UHF TELEVISION BROADCASTING ANTENNAS* 



By 

D. W. Peterson 

RCA Laboratories, 
Princeton, N\ J. 

Summary — Determination of the power gain and radiation pattern of 
a modem television broadcasting antenna is difficult after it is installed. 
However, if ordinary coverage surveys do not give the expected results, 
tests on the installed transmitting antenna may be necessary to measure 
the performance of this part of the system. This paper undertakes a 
summary of the pertinent propagation information and outlines a pro- 
cedure for field checking antennas that requires the minimum effort con- 
sistent with obtaining reliable results. Propagation effects due to terrain 
shadowing, clutter, and the earth's curvature, are avoided as completely 
as possible. 

A check of antenna performance does not constitute even a partial 
coverage survey. Coverage surveying is quite another matter which includes 
the effects of terrain shadowing , clutter, and the earth's curvature. 

Introduction 

A NY post-installation check of a UHF television broadcasting 
i=i\ antenna is generally recognized as much less accurate than 
factory tests made under relatively ideal conditions. Even 
at a properly designed factory test-site, antenna gain and pattern 
measurements unavoidably include earth reflection effects, although 
these effects are now usually negligible 1 . The philosophy adopted for 
post-installation antenna checking avoids a large part of the com- 
plication of propagation theory. There remains a simplified theory 
that must be applied. UHF broadcasting antennas need rarely be 
checked beyond a distance of ten miles. By thus restricting the dis- 
tance and carefully choosing propagation situations, the theoretical 
requirements can be limited to simple free-space and plane-earth the- 
ories. Certain tests are required to know which theory to apply to a 
specific field measuring situation. If, in spite of the obvious superiority 
of factory test-site measurements, field testing of an antenna becomes 
necessary after installation, this may be done with reasonable reliability 

* Manuscript received September 9, 1958. 

1 E. H. Shively, "Pattern-Testing the TF-24B UHF Antenna," Broad- 
cast News, Vol. 69, p. 42, Mav-June, 1952, and E. H. Shively and L. D. 
Wetzel, "Measurements of RCA UHF TV Antennas," Broadcast News, 
Vol. 82, p. 14, February, 1955. 
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provided care is exercised and the limitations imposed by propagation 
effects are fully recognized. 

Problems Associated with Field Measurement 

In order to make a post-installation check of an antenna, one must 
define the size, location, and character of areas, known as first Fresnel 
(earth reflecting) zones, that may cause significant earth reflection 
effects. Also, the first Fresnel (shadowing) zone must be described 
quantitatively to know that significant terrain shadowing effects have 
been avoided. For the exceptionally smooth-earth situation there must 
be quantitative propagation laws for prediction, and a smooth-earth 
criterion is required to determine when to apply smooth-earth theory. 
Also, the near-field region must be determined. Means of establishing 
these factors are provided in this paper in the form of easily used 
curves and formulas. There is also a brief review of the applicable 
theory and recommended procedure for field measurement of antenna 
performance. 

Since the object is to field-test the antenna, unpredictable propaga- 
tion effects must be avoided as far as possible. This means that terrain 
shadowing, clutter effects, and shadowing from the earth's curvature 
should be completely avoided. In other words, all field measurements 
for post-installation antenna checking should be made in unmistakably 
uncluttered, line-of-sight locations. 

However, one propagation effect, namely earth reflection, cannot 
be completely avoided. This effect can cause an increase above the 
free-space value of field strength (electric intensity) up to 2 to 1 and 
in rare cases even more. Also, earth reflection may reduce field strength 
10 to 1 or more below the free-space value. Since the effect cannot be 
avoided, it must be taken into consideration. Examination of the 
field strength versus height characteristic often gives clues as to the 
existence and magnitude of earth reflection. However, these clues 
cannot always be relied upon since the limited height range examined 
may not reveal the entire character of the interference pattern result- 
ing from earth reflection because of the complexity from more than 
one effective reflecting surface and the resulting multiple rays. It 
should be kept in mind that significant reflecting surfaces (portions of 
hills) may occur at any distance or may be situated away from the 
direct path. 

Because of the unavoidable propagation effects mentioned, there is 
no positive assurance that even carefully selected measuring sites will 
yield the predicted field strengths. Nevertheless, field strengths in the 

i 
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order of the free-space value are usually found over rough terrain 
provided the precautions discussed below are observed. 

Some idea of the field strength to expect over smooth terrain can 
also be reached. However, a complication results from the fact that 
exceptionally smooth terrain may cause a lobed two-ray interference 
pattern. When this condition is suspected, the theoretical field strength 
versus height distribution must be calculated and measurements must 
be made to determine whether or not the condition exists. Sufficiently 
smooth earth will limit the distance of usable measuring sites to those 
sites where interference maxima and minima can be reached with 
feasible receiving-antenna heights. 

UHF propagation is subject to severe shadowing effects from ter- 
rain obstacles. These effects can and should be avoided by maintaining 
at least first Fresnel (shadowing) zone clearance above terrain ob- 
stacles. Similarly, the effect of the local clutter of trees and buildings 
at the receiving site is to be avoided. Both the terrain shadowing and 
the clutter effects are highly frequency dependent. Since they can be 
avoided, these effects need not unduly complicate field performance 
checks of a broadcasting antenna. 

Measurements of field intensity can be made too near as well as 
too far with respect to the transmitting antenna. Thus, the near-field 
region should be avoided because within this distance range the pattern 
has not fully formed. On the other hand, when measurements are 
attempted at too great a distance the field intensity will be time 
variable, another condition to avoid. 

In cases where the azimuth pattern as measured at the factory 
test-site is accepted, only one radial line need be examined. The 
distance range of measurements along this one radial will usually 
include radiation from the entire useful part of the elevation pattern, 
with sufficient sampling detail to establish the general character of 
the elevation pattern. Often only the maximum of the elevation pattern 
need be checked, preferably in several receiving sites. When all or 
part of the azimuth pattern must be checked, the process of examining 
the elevation pattern may be repeated along additional radial lines. 
A reasonable compromise to reduce the amount of measurement for 
a horizontal pattern check is to check pattern circularity by measure- 
ments on the elevation pattern maximum only. 

When field strengths are below expected values, all relevant theo- 
retical possibilities must carefully be considered in accounting for the 
condition. It is unfortunate but true that most of the propagation 
effects reduce rather than increase field strength with respect to the 
predictions of simple theory. 
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Free-Space Theory of Propagation 

When propagation is presumed to occur in empty space, attenuation 
results from the divergence of energy from a point source. Antenna 
engineers like to think of an isotropically radiating source at the 
point; consequently, power flow through the space can be thought of 
as radial and with uniform power density flowing across any sphere 
concentric with the point source. An isotropic source has come into 
wide use as the reference for antenna gain. However, because of the 
ready availability of half-wave dipole antennas, this kind of reference 
is used instead for television antennas. 

Transmitting-antenna power gain is commonly expressed as a power 
ratio, specifically the ratio of power radiated in free space on the 
elevation pattern maximum to the power maximum of a reference 
half-wave dipole, assuming that the power inputs are equal. The 
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Fig. 1 — Azimuth pattern, 
azimuth field intensity patterns of most broadcasting antennas are 
not precisely circular although they usually closely approach circularity. 
The quoted antenna power gain, G, is the gain that an antenna with 
the same radiated power and the same elevation pattern would have 
if the azimuth field intensity pattern were precisely circular with 
unity magnitude. The circle is the root-mean-square circle of the 
azimuth pattern as shown in Figure 1. The true gain, G$ , in a direc- 
tion specified by 4> 1 can be related to the quoted gain thus: 

G H = [f(^)r-G (1) 

where G^ is the antenna power gain in the direction specified by 

f(<pi) is the azimuth field intensity pattern evaluated for the 
direction <j> lt 
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G is the r-m-s antenna power gain. 

Because of the use of a half -wave reference dipole for power gain, 
the free-space field intensity on the maximum of the dipole pattern 
will be useful. 

137.6 vp; 

E 0 = (2) 

r 

where E 0 is the free-space field intensity in millivolts/meter for 
a half-wave dipole, 

P d is the power input to the dipole in kilowatts, 

r is the distance in miles between the antenna and a point 
in space. 

In the direction specified by <j> ± the field intensity on the maximum 
of the elevation pattern from an antenna of power gain, G, is 

137.6 VGP 

A*, =" /(<£i) (3) 

r 

where E <t>l is the field intensity in the <j> x direction in 
millivolts/meter, 

P is the power input of the transmitting antenna in 
kilowatts. 

The theoretical free-space field intensity is expected under certain 
conditions when the smooth-earth criterion (discussed later) is not 
satisfied. Before this becomes a useful tool for checking an antenna 
in the field, methods of measurement and the necessary conditions must 
be understood. On the other hand, when the smooth-earth criterion is 
satisfied, plane-earth theory must be considered. 

Ray Theory of Propagation Over Plane Earth 

Propagation from a transmitting antenna located at P 2 to a point 
I\ over smooth, reflecting earth is described in Figure 2. The earth 
is assumed to be a relatively smooth plane. The lobing which results 
from interference of the direct and earth-reflected rays is described 
and formulas are given for the calculation of field intensity over 
smooth, plane earth. 
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Fig. 2 — Plane-earth reflection. 



According to simple plane-earth ray theory, radiation from P 2 can 
reach P x by way of two rays, namely P 2 P 1 (the direct ray) and P- 2 P s Pi 
(the earth-reflected ray) . The incident grazing angle fa of the earth- 
reflected ray must equal the angle fa of the ray upon reflection. With 
horizontal polarization, as used for television broadcasting in the 
United States, there is a 180° phase reversal upon reflection at P 3 , 
and there is complete reflection (magnitude of the reflection coefficient 
is unity) for smooth-earth conditions. 

If the transmitting antenna consists of a half-wave dipole with a 
circular elevation pattern, the field intensity at P t will vary with height 
hi according to a simple periodic law, resulting from the alternately 
constructive and destructive interference of two rays, provided a 
simplifying restriction is made. The restriction is that the grazing 
angle f will be small (less than 10°) so that the reflection coefficient 
will remain close to unity and so that the inverse-distance law of 
free-space propagation will result in only small differences in the 
magnitudes of ray 1 and ray 2 field intensity contributions. The law 
of field intensity variation is, for horizontal polarization and ^ < 10°, 



This formula, shown graphically in Figure 3, gives a limit on the 



E %1 = -E 0 V2-2cosA, 



(4) 



where E hl is the field intensity at height 7ij Oh i s related to A in 
Equation (5)), 

A is the difference in path length in electrical degrees 
between rays 1 and 2. 



E 




A = f(h,,h 2l ^,X) 
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Fig. 3 — Field-intensity variation with height. 
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maximum possible excursion of field intensity variation with height 
for any transmitting antenna as long as the receiving site is below 
the maximum of the elevation pattern. When the antenna has power 
gain by virtue of a narrowed elevation pattern, the variation of field 
intensity with height will be similar as long as both ray 1 and ray 2 
are near the maximum of the elevation pattern. Under most other 
practical conditions, the excursions about the free-space field intensity 
will be less. With high-gain shaped-elevation-pattern antennas, the 
ray 1 and ray 2 contributions must be calculated separately if the 
magnitudes of the maxima and minima to be found with height are 
desired because portions of the pattern exhibit rapid changes with 
elevation angle. However, it usually suffices to approximately locate 
the maxima and minima with Equation (4) and then depend upon 
measurement in the field for the magnitudes. For a well-chosen meas- 
uring site, the average of these two values will be the free-space field. 

The path difference A depends upon the geometry shown in Figure 
2, and the wavelength and can be calculated from the formula 2 

k A = 0.136 (5) 

# dk 

w 

where /i x and lu are the measuring and transmitting antenna 
heights in feet, 

d is the distance between P t and P., in miles, 

980 

A is the wavelength in feet. A = , where / is in 

megacycles/sec. 

Test Equipment for Field Measurement 

A reliable field-intensity meter is, of course, indispensable. When 
reliable measurements are needed, it is advisable to check the instru- 
ment against laboratory generators or against a bolometer. This, of 
course, only checks the calibrating generator in the field-intensity 
meter. The entire instrument can be checked by establishing and 
measuring known field intensities. 

An over-all check of the field intensity meter can be made in a large, 
flat, clear area (say a 10 acre field) by establishing a radiated field 
with a signal generator and a matched A/2 dipole. The field-intensity- 
meter dipole height should be adjusted to pass through both inter- 
ference maxima and minima at a distance great enough so the direct 

York, N^Y^l 953 f^ SeU ' UHF Propagation, John Wiley & Sons, New 
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and earth-reflected rays are essentially the same length. The free-space 
field which has been established (actually the field strength which 
would have existed without the presence of the reflecting earth) can 
be calculated from Equation (2). Of course, transmission-line attenua- 
tion must be subtracted in the calculation of established field strength. 
The free-space field intensity from the measurement will be 

■E'min + E max 

E 0 = , 

2 

where E min is the field intensity on the interference minimum, and 
£7 max is the field intensity on the interference maximum. 

When measurements are made in hilly terrain or in the vicinity 
of possible reflectors (this cannot always be avoided), the standard 
dipole furnished with the field-intensity meter should be replaced with 
a corner reflector or paraboloid antenna. This can be calibrated against 
the dipole in a large, flat, clear area by direct substitution, taking 
care that both antennas are placed successively in the same field 
intensity. The field-intensity distribution in the substitution test area 
should first be checked for constancy by probing with the dipole. 

Procedure for Field Measurement op Performance 

Two kinds of useful measurements associated with two distinct 
types of propagation paths can be made in an effort to check antenna 
performance in the field. 

The first kind consists of field strength versus height measurements 
with conditions such that two-ray (plane-earth) theory clearly applies. 
There will be no question about the existence of first Fresnel (shadow- 
ing) zone clearance since there will be no hills or other shadowing 
obstacles. It will be necessary to estimate the theoretical heights of 
maxima and minima and measure a field strength versus height dis- 
tribution from which to calculate (by averaging a maximum and a 
minimum) the free-space field intensity. The area surrounding Miami, 
Florida, and the Illinois prairie are illustrative of places where the 
terrain is flat enough to apply smooth-earth theory. Paths over bodies 
of water can also be included in this category. 

The second kind consists of field strength versus height measure- 
ments over terrain so rough that the two-ray (plane earth) theory 
as tested by the Rayleigh roughness criterion does not apply. Over 
rough terrain, there must be a test for first Fresnel (shadowing) zone 
clearance. 
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Both kinds of measurements must be made under line-of-sight 
conditions without the clutter of trees or buildings in the vicinity 
of the receiving antenna. They should be made as close as possible 
(outside the near-field region) to the transmitting antenna to avoid 
spherical earth effects. Other kinds of propagation paths will be 
encountered but, for the purpose of antenna checking, are to be 
avoided. 

A straightforward procedure for checking television broadcasting 
antenna performance will now be outlined: 

1. Measure the power fed to the broadcasting antenna after checking 
the VSWR of the antenna. An independent substitution wattmeter 
is preferable to simply relying on built-in power-monitoring devices. 
Subtract transmission line loss to learn antenna input power. 

2. Select a number of measuring sites outside the near-field region 
(Equation (14)), but no further than necessary from the transmit- 
ting antenna, which are completely uncluttered with buildings or 
trees. From these sites, the transmitting antenna should be clearly 
seen to be well above all intervening terrain features. On-the-spot 
observation is necessary. If clutter at the receiving site is unavoid- 
able, the receiving antenna should be kept well above the clutter. 

3. Plot elevation profiles 3 for the intervening terrain for each potential 
test site. 

4. Calculate the size and location of the first Fresnel (earth reflecting) 
zone, Equations (10), (11), and (12), and test this area by the 
smooth-earth criterion (Equation (6)). It must be recognized that 
terrain features outside the first Fresnel (earth reflecting) zone 
may alter the received field intensity. 

5. If the smooth-earth criterion indicates smooth earth, calculate 
the theoretical location of field-intensity interference maxima and 
minima as a function of receiving antenna height (Equation (4)). 
If no theoretical maximum is accessible, the site should not be used. 

6. If the smooth-earth criterion indicates rough earth, test the inter- 
vening terrain for first Fresnel (shadowing) zone clearance and 
abandon the site if the clearance test is not satisfied. The test 
should be applied to hills, buildings, and trees. (Figure 4) 

7. Measure and plot field intensity versus height in the acceptable 
sites and compare data with appropriate theory. 

The procedure outlined has been aimed at avoiding unpredictable 



3 If topographic maps are not available from either Federal or State 
sources, salient features can sometimes be measured with an altimeter. 
See H. Reed and C. Russell, op. ext., p. 52. 
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propagation effects. This is not always possible; however, it will 
usually be possible to find sites which are satisfactory. The weakness 
of a measuring site will not always be apparent. Sometimes the field 
intensity versus height recordings will suggest likely shortcomings 
of a measuring site. 



ANTENNA I 




- ANTENNA 2 

\ .FRESNEL ZONE BOUNDARY 

-TERRAIN ELEVATION PROFILE 



USE WHERE 0 << b 



USE WHERE Q = b 




Pig. 4 — First Fresnel (shadowing) zone clearance. 



Expected Field Intensity 

When the smooth-earth criterion indicates smooth earth in the 
first Fresnel (earth reflecting) zone, values ranging from the free- 
space field intensity to twice this value are to be expected at heights 
where interference maxima occur. Either the reflection coefficient can 
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be lower than unity or ray 1 and ray 2 of Figure 2 may differ signifi- 
cantly in length to account for this range of values. The field intensity 
versus height data should confirm the applicability of smooth-earth 
theory by the presence of interference maxima and minima. Several 
field intensity versus height distributions are desirable at each meas- 
uring site. Recordings at constant height of field strength versus 
distance are also valuable. 

When the smooth-earth criterion indicates rough earth in the 
first Fresnel (earth reflecting) zone, free-space field intensities are 
to be expected. If the surroundings of the receiving site are hilly, it 
will be usual to observe complicated field intensity versus height dis- 
tributions which will be a consequence of the immediate surroundings. 
Under these circumstances field intensities both above and below the 
free-space value are possible. 

Another condition which may be encountered is intervening terrain 
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Fig. 5 — Measured field strength versus height in hilly terrain. 

in the first Fresnel (earth reflecting) zone which satisfies the smooth- 
earth criterion, accompanied by cluttered and/or hilly surroundings at 
the receiving site. Then there is the possibility that there will be an 
average reduction of field intensity below free-space values because of 
destructive interference throughout the measuring site. In other words, 
there may be superimposed on the observable maxima and minima an 
average trend of increasing (or decreasing) field intensity with height 
which is not fully observed because of the limited height range which 
is accessible. When this happens, there is no way (except exploration 
at unreasonable heights) of knowing the extent of the interference 
hidden in the measurements. Thus, in the typical field strength versus 
height distribution illustrated in Figure 5, examination of field 
strength up to 30 feet would fail to reveal that there had been a 
reduction in field strength at all heights up to 30 feet as a result of 
destructive interference. The pattern of destructive interference 
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would only have been apparent upon examination of received field 
strength up to 200 feet in height. The rapid variations with height 
are the result of reflecting surfaces on relatively nearby hills whereas 
the slower variation results from a more remote smooth plane in the 
first Fresnel (earth reflecting) zone. 

"Before-and-After" Comparisons 

There may be a requirement to test the improvement in antenna 
performance upon changing the transmitting antenna while keeping 
the same center of radiation; e.g., there may be an increase in gain, 
a change in pattern shape, or a shift in frequency. Presumably data 
both before and after the antenna change can be obtained. The problem 
here is basically easier than the absolute test of a single antenna 
because with due precaution the propagation effects can be virtually 
eliminated from the comparison. However, there are important ex- 
ceptions. 

When the first Fresnel (earth reflecting) zone is smooth accord- 
ing to the smooth-earth criterion, it will be necessary to proceed with 
care. It should be recognized that the direct and earth-reflected rays 
may come from significantly different parts of the elevation pattern. 
Thus, when there is a change of pattern shape, the contributions from 
the earth-reflected ray and direct ray both change in magnitude but 
not necessarily by the same factor. Unless this is understood there 
may be apparent inconsistencies in before-and-after comparisons 
even though the frequency and the height of the center of radiation 
remain the same and the antennas are functioning normally. 

If the antenna change includes even a small change in location there 
will be drastic changes in the detailed field intensity distribution in 
rough terrain. This can easily be true even in areas of fringe reception 
at great distances. Furthermore, a shadowed fringe-area valley may 
receive median field intensities which are changed even though the 
same transmitting antenna is moved only a few miles transversely with 
respect to the propagation path. 

When an antenna change is accompanied by a shift in frequency 
there will be propagation changes which must be given due regard in 
before-and-after comparisons. This will be true even though the 
elevation pattern and the center of radiation remain unchanged. 

Another kind of precaution relates to the time dependence of 
received field intensity beyond the horizon of smooth earth or in the 
deep shadows of rough terrain. Because there is often substantial time 
variation, it is well in before-and-after comparisons to include only 
line-of-sight receiving sites. When necessary, co-channel interference 
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should be kept to a minimum by use of a unidirectional receiving 
antenna. 

Comparisons can be made by, say, eight radial recordings of field 
strength measured at constant height as for a coverage survey. How- 
ever, if the only object of the test is demonstration of the improvement 
or change of the transmitting antenna performance, it can be done 
much more easily by excluding cluttered surroundings and making 100 
foot recording samples at line-of-sight locations chosen to measure on 
significant parts of the elevation pattern. Of course, the measurements 
will be made in exactly the same places before and after the antenna 
change. 

Accuracy 

The accuracy of field measurements as an indication of broadcasting 
antenna performance is determined by (1) the care exercised in 
choosing measuring sites and (2) the quality and amount of the data. 

The matter of expected accuracy for absolute field-intensity meas- 
urements over commonly encountered kinds of terrain merits discus- 
sion. Field-intensity meters depend, for absolute accuracy, upon a 
bolometer or a crystal diode in the calibrating generator. High-quality 
UHF signal generators are claimed to have absolute voltage accuracies 
of ± 1.5 to 2 decibels. One would hardly expect better from a field- 
intensity meter calibrating generator. Furthermore, there will ordi- 
narily be mismatch error in field intensity measurements unless the 
A/2 dipole effective length factor incorporates mismatch correction. 
Thus, an accuracy of ± 2 decibels for a field-intensity meter would, by 
present practice, be good. More can be expected with respect to relative 
accuracy at a given time and frequency since the instruments depend 
upon accurate and reliable piston attenuators for relative values. Piston 
attenuator accuracies in the order of tenths of a decibel are easily 
attained. 

Other accuracies which are involved in field-intensity measurements 
include the transmitting antenna azimuth and elevation patterns, 
transmitter power, transmission line attenuation, diplexer attenuation, 
and transmitting antenna gain. Any field evaluation of transmitting 
antenna performance should be based upon knowledge of all of the 
quoted factors and their probable accuracies. 

Appendix 

First Fresnel (Shadowing) Zone Clearance 

If free-space propagation theory is to be applicable, there must be 
a certain minimum clearance between the terrain and a ray from the 
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transmitting antenna to the measuring site. The first Fresnel (shad- 
owing) zone is bounded by an ellipsoid with the transmitting antenna 
and the measuring point at the foci. When there is first Fresnel 
(shadowing) zone clearance over terrain which is rough by the 
smooth-earth criterion, free-space propagation is expected. The re- 
quired clearance can be determined from Figure 4. 

Test to Determine Validity of Ran/ Theory Over Plane Earth 

It has been convenient to treat earth reflection as though occurring 
at a point (ray theory) . Actually, there is a significant reflecting area, 
elliptical in shape, known as the First Fresnel (earth reflecting) zone. 4 
This area can be located and defined. The length of the ellipse is 2a?!, 
the width is 2y lt as shown in Figure 6. If the first Fresnel (earth 
reflecting) zone surface is smooth enough to satisfy the Rayleigh 

POINTS P, AND P 2 LIE IN 
VERTICAL EAST-WEST PLANE 




SOUTH DOWN 

Fig. 6 — Reflecting zone. 

smooth-earth criterion, it is probable that the simple ray theory 
applies. 

According to the smooth-earth criterion, a reflector behaves essen- 
tially like a mirror 5 if 



H< (6) 

16 sin f 

where H is the height of a surface irregularity, 

$ is the angle of incidence. 

According to this rule-of-thumb, reflection is specular or mirror-like 



4 K. A. Norton and A. C. Omberg, "The Maximum Range of a Radar 
Set," Proc. I.R.E., Vol. 35, p. 4, January, 1947. 

5 H. Reed and C. Russell, op. ext., pp. 236-238. 



070 RCA REVIEW December loss 

when the height, H, Figure 7, of surface irregularities satisfies the 
smooth-earth criterion. 

When reflection can be assumed to be specular, based on the smooth- 
earth criterion, the reflection coefficient is assumed to be unity Then 
the field-intensity contribution of either the direct or the earth- 
reflected ray can be computed from the free-space formula 

137.6 VF 

®o — ■ millivolts per meter, 

r 

where r is the length, in miles, of the rays. The interference maxima 
and minima will be, respectively, the sum and difference of the two 
contributions. 

First Fresnel (Earth Reflecting) Zone 

The formulas for computing the dimensions and location of the 
first Fresnel (earth reflecting) zone on plane earth follow. Refer to 
figure 6 for the geometry and symbols. 




surface irregularity 



Fig. 7 — Roughness criterion geometry, 
wavelength, A, is in feet 
frequency, /, is in mc/sec. 
heights h x and lu are in feet 
distances x 0 , z t , y u d u d.,, and d are in miles. 
The wavelength is 

980 



A=- 



(7) 



The grazing angle for the earth-reflected ray is 

hi 



4> = tan -1 



5280 <Z, 



(8) 
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If h 2 » h lt A. the distance to the earth-reflected-ray reflection point is 

h 2 

The distance to the center of the first Fresnel (earth reflecting) zone 
ellipse is 



x 0 sa d ± 



— 1 

i 1 sin xp _] 



1+ |. (10) 

2 h. 



The length of the ellipse (major diameter) is 



2 x x = 0.00378 



1 + 

4 hi sin <]/ _ 



(ID 



sin 3 >p 

The width of the ellipse (minor diameter) is 

2y 1 = 2 x x sin (12) 

While these formulas provide a useful guide, it should be recognized 
that significant reflections may occur outside the first Fresnel (earth 
reflecting) zone. 

Transmitting Antenna Pattern and Power Gain 

Since antenna gain is quoted in terms of r-m-s power gain, G, 
antenna azimuth and elevation patterns must be taken into account 
in the estimation of expected field intensity at a specific measuring 
site. Designating the azimuth pattern as /(<£) and the elevation 
pattern as f{6), these patterns may be related to the coordinates of 
Figure 8. If the azimuth direction angle is fa and the elevation angle 
0! for a specific measuring site, the theoretical free-space field intensity 
at a point r x miles from the transmitting antenna will be 

137.6 VGT 

E^6 in = f(fa)f(0 1 )- (13) 

Avoidance of the Near-Field Region 

The near-field region of the radiation pattern should be avoided 
since only far-field patterns are ordinarily published. Thus, measure- 
ments should not be made closer than 
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(b) 

Fig. 8(a) — Azimuth pattern in horizontal NSEW plane; 
(b) — Elevation pattern in east plane. 



where d miB is the minimum measuring distance in feet, 

W is the transmitting antenna aperture (length of 
radiating portion in feet), 



(14) 



\ is the wavelength in feet. 
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